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ABSTRACT  OF  THESIS 


> 

> 

AN  ANALYSIS  OF  SIMULATED  AND  ACTUAL  DMSP  SSM/T-2  BRIGHTNESS 

TEMPERATURES 

I 

Defense  Meteordo^cal  Satellite  Program  (DMSP)  Special  Sensor  Microwave/T-2 
(SSM/T-2  or  T-2)  global  brightness  temperature  maps  are  related  to  some  of  the  known  ob¬ 
served  general  circulation  features  of  the  atmosphere  and  surface  characteristics  of  Earth. 

The  brightness  temperatures  from  two  re^ons  with  different  large  scale  circulations  are 
compared  and  contrasted. 

The  concept  of  the  invariance  of  the  weighting  functions  with  respect  to  water  va¬ 
por  overburden  is  tested.  It  is  demonstrated  that  the  levels  where  the  183  ±  1  and  3  ^ 

GHz  brightness  temperatures  match  the  thermodynamic  temperature  are  levels  of  ap- 
projdmately  equal  overburden.  Water  vapor  overburden  (integrated  water  vapor  above  a 
given  level)  was  retrieved  and  mapped  on  a  constant  pressure  surface  from  T-2  brightness  | 

temperatures  and  atmospheric  temperature-pressure  data.  A  simple  physical  retrieval 
scheme  based  on  the  approximately  equal  overbrnden,  irrespective  of  the  moisture  profile, 
for  the  183  ±  1  GHz  frequency  was  developed  and  demonstrated.  This  simple  retrieval  ^ 

coidd  be  used  to  initialize  more  sophisticated  retrievals. 

Simulation  studies  of  some  factors  which  affect  brightness  temperatures  at  the 
SSM/T-2  frequencies  were  conducted  and  the  residts  discussed.  These  factors  included 
moisture  profiles,  backgrounds  with  different  emissivities,  and  low  level  water  clouds. 

FRANK  A.  LEUTE  IV 
Department  of  Atmospheric  Science 
Colorado  State  University 

Fort  Collins,  Colorado  80523  ^ 

Summer,  1993 
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Chapter  1 


INTRODUCTION 

1.1  BACKGROUND 

Water  vapor  plays  a  vital  role  in  establishing  our  present-day  climate.  It  is  the  ^ 

most  significant  contributor  to  the  observed  greenhouse  effect  of  the  planet  (Stephens  and 
Tjemkes,  1992).  Water  vapor  links  the  processes  of  evaporation,  cloud  formation,  and 
precipitation.  Latent  heat  release  plays  a  large  role  in  the  dynamical  flow  fields  of  the  > 

tropical  atmosphere  and  knowledge  of  the  distribution,  transport,  and  divergence  of  water 
in  all  phases  is  critical  to  work  pertaining  to  the  difficult  problem  of  cloud  feedback  and 
its  role  in  climate  change  (Stephens,  persoxial  communication).  ^ 

The  source  of  water  vapor  is  primarily  at  the  earth’s  smface  whereas  its  sink  is  in 
the  process  of  precipitation.  Combining  these  sources  and  sinks  with  the  atmospheric 
circulation  means  the  distribution  of  water  vapor  in  the  atmosphere  is  highly  variable 
both  in  spatial  and  temporal  scales  (Prabhakara  and  Dalu,  1980).  Thus,  it  is  necessary 
to  establish  a  dense  network  for  observing  water  vapor  aid  its  distribution.  Present  day 
observations  are  largely  based  in  the  use  of  radiosonde  data  and  these  data  lack  in  the 
necessary  global  coverage  needed  for  rmderstanding  the  role  of  water  vapor  in  climate.  ^ 

Water  vapor  also  affects  the  exchange  of  radiation  between  the  atmosphere  and  space 
directly  by  its  influence  in  atmospheric  emission  and  indirectly  by  its  role  in  cloud  for¬ 
mation  processes  (Bates  and  Stephens,  1991).  GEWEX  (the  Global  Energy  and  Water  ► 

Cycle  Experiment,  part  of  the  World  Climate  Research  Program)  has  the  basic  theme  of 
developing  a  quantitative  rmderstanding  of  how  water  shapes  the  energy  budget  of  the  cli¬ 
mate  system  of  earth.  The  connections  and  relationships  between  the  global  hydrological  I 

cycle  and  the  energy  budget  are  critical  to  the  question  (and  potential  problem)  of  climate 
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change  (Stephens  and  Tjemkes,  1992).  The  typical  climate  change  scenario  (2  x  CO2) 
predicts  a  global  warming  of  approximately  twice  (or  treble,  depending  on  the  model)  the 
C02-induced  warming  with  water  vapor  feedback  than  without  this  feedback  (Houghton 
et  aL,  1990).  The  first  part  of  the  water  vapor  feedback  is  due  to  the  direct  greenhouse 
effects  of  doubled  CO2  which  directly  leads  to  increased  sea  surface  temperature  (SST). 
The  increased  SST  increases  evaporation  in  the  boimdary  layer  which  increases  humidity 
in  the  boundary  layer.  The  increase  of  water  vapor  in  the  boundary  layer  through  its 
greenhouse  effect  increases  the  SST.  The  second  part  of  the  feedback  process  was  hy¬ 
pothesized  by  Lindzen  (1990).  Water  vapor  feedback  is  also  thought  to  be  modulated  by 
the  changes  in  the  vertical  distribution  of  water  vapor.  Lindzen  argues  that  the  increased 
warming  near  the  ground  results  in  increased  and  deeper  cumulus  convection.  This  in¬ 
creased  cumulus  convection,  he  clmms,  leads  to  drying  of  the  upper  troposphere  (above 
5  km)  and  elevation  of  the  altitude  at  which  convected  heat  is  deposited.  Since  green¬ 
house  absorption  is  most  important  above  5  km  these  cumulus  effects  are  negative  and 
hence,  should  diminish  the  effect  of  CO2  warming.  In  contrast  to  this  view,  present  day 
climate  models  actually  produce  a  moistening  of  the  upper  troposphere  associated  with 
this  convection.  Both  hypotheses  are  depicted  in  Figure  1.1  (Stephens,  personal  commu¬ 
nication).  Attempts  to  study  quantitatively  the  distribution  and  transport  of  water  vapor 
in  the  atmosphere  have  been  mcide  using  the  global  radiosonde  network  (e.g.  Peixoto  and 
Oort,  1983).  Unfortimately,  these  studies  have  met  limited  success  since  these  radiosonde 
observations  are  confined  msdnly  to  land  areas  since  this  is  where  the  radiosonde  stations 
are  primarily  located.  Lnportant  atmospheric  phenomenon  restricted  to  oceanic  areas  are 
thus  largely  ignored  (e.g.  many  aspects  of  El  Nino,  including  westerly  wind  bursts  and  the 
location  of  moisture  convergence/divergence).  Early  work  on  transports  based  on  infrared 
(IR)  water  vapor  data  have  been  reported  by  Wittmeyer  (1990).  However,  Wittmeyer 
and  Vonder  Haar  (1993)  and  others  report  serious  adverse  cloud  impacts  in  baroclinlc 
zones.  Schwartz  and  Doswell  (1991)  contend  that  more  observations  are  essential  in  both 
time  and  space  to  improve  mesoscale  weather  forecasting  and  further,  that  remote  sensing 
systems  have  yet  to  approach  the  radiosonde’s  capability  to  resolve  the  vertical  thermody¬ 
namic  structiire  of  the  atmosphere.  While  the  later  may  or  may  not  remain  true  in  the  age 


Figure  1.1:  Hypotheses  on  water  vapor  feedbztck  mechanisms  (Stephens,  personal  com¬ 
munication).  The  left  portion  of  the  diagram  is  the  feedback  as  presently  perceived,  with 
increasing  temperature  leading  to  increased  water  vapor  largely  in  the  boundary  layer.  The 
right  portion  of  the  diagram  is  the  additional  modification  to  the  feedback  by  increasing 
the  upper  tropospheric  moisture  through  convection.  > 
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of  new  microwave  (millimeter  wave)  instruments,  the  former  may  be  improved  drastically 
by  these  new  remote  sensing  systems.  Certainly,  the  future  lies,  in  many  respects,  in  the 
advancement  of  the  capabilities  of  remote  sensing  end  the  application  of  these  capabilities 
to  the  measurement  of  atmospheric  phenomena. 

Satellites  are  used  for  many  applications  of  remote  sensing.  The  theory  and  practice 
of  remote  sensing  spans  many  disciplines  and  is  intimately  related  to  inversion  theory  (see 
Twomey,  1977).  The  deployment  of  satellites  expanded  the  opportunities  for  deriving  the 
distribution  of  water  vapor  over  the  entire  globe  (as  opposed  to  the  limitations  with  land- 
and  ship-based  radiosondes).  Moreover,  many  of  the  satellites  launched  in  the  1970s  and 
1980s  were  restricted  to  measurements  of  the  total  precipitable  water.  The  launch  of  the 
Special  Sensor  Microwave/ Water  Vapor  Sounder  (SSM/T-2  or  T-2  which  is  described  in 
more  detsdl  in  the  next  chapter)  heralds  a  new  era  in  water  vapor  studies  with  its  apparent 
capability  of  profiling  atmospheric  water  vapor  although  not  with  the  vertical  resolution 
of  current  radiosondes. 

Observing  system  simulation  experiments  (OSSEs)  were  carried  out  by  Hoffinan  et  al. 
(1990)  to  attempt  to  detemune  the  impact  of  the  Special  Sensor  Microwave/Temperature 
(SSM/T)  sensors,  in  addition  to  other  sensors.  In  these  experiments  the  addition  of  the 
SSM/T-2  and  SSM/T-1  data  improves  moisture  analysis,  especially  in  the  tropics  and  the 
extratropics  of  the  Southern  Hemisphere.  Root  mean  square  (RMS)  errors  for  relative 
humidity  forecasts  were  decreased  from  29%  to  22%  in  the  OSSEs  with  the  SSM/T-2 
data  included.  In  the  Southern  Hemisphere,  the  SSM  data  improved  the  500  hPa  height 
forecasts  by  12  hom-s. 

1.2  MICROWAVE/MILLIMETER  WAVE  SATELLITE  RETRIEVALS 

The  potential  benefits  of  using  passive  microwave  sensors  to  measure  atmospheric 
parameters  was  first  recognized  in  the  1940’s.  Nonetheless,  it  was  only  during  the  early 
1960’s  when  the  first  microwave  radiometer  was  used  on  a  spacecraft  for  the  purpose 
of  atmospheric  remote  sensing.  The  U.S.  Mariner  2  Venus  Probe  in  1962  took  the  first 
satellite  microwave  emission  observations  of  a  planetary  atmosphere.  The  two-channel 


radiometer  determined  an  upper  limit  on  water  vapor  in  Venus’  atmosphere  and  retrieved 
a  planetary  surface  temperature  near  600°K  (Barath  et  al.,  1964).  The  first  U.S.  satellite 
to  make  microwave  observations  of  the  earth’s  atmosphere  was  Nimbus  5  in  1972. 

Microwave  temperature  retrieval  on  satellitss  operated  by  National  Oceanic  and  At¬ 
mospheric  Administration  (NOAA)  began  with  the  Nimbus  5  (Nimbus  E)  Satellite  Mi¬ 
crowave  Spectrometer  and  continued  with  the  Nimbxis  6  Scanning  Microwave  Spectrometer 
(SCAMS).  The  SCAMS  used  three  channels  (52.85  GHz,  53.85  GHz,  and  55.45  GHz)  in 
the  oxygen  absorption  band  between  50  and  60  GHz.  The  follow-on  to  the  SCAMS  is  the 
Microwave  Soimding  Unit  (MSU).  It  serves  as  a  complement  to  the  infrared  temperature 
soimding  imit  called  the  High  Resolution  Infrared  Sounder  2  (ElRS/2).  The  MSU  has  four 
channels  (50.30  GHz,  53.74  GHz,  54.96  GHz,  and  57.95  GHz)  in  the  oxygen  absorption 
band.  The  MSU  is  a  ‘Dicke’  radiometer  where  the  output  of  the  receiver  is  proportional 
to  the  (hfference  between  the  brightness  temperattne  of  the  scene  being  viewed  and  the 
temperature  of  the  intemai  radiation  source.  'Ihe  primary  purpose  of  the  MSU  is  to  make 
temperatme  soimdings  in  the  presence  of  clouds  although  MSU  brightness  temperature 
data  are  also  currently  being  analyzed  to  understand  certain  aspects  of  the  Earth’s  climate 
(Spencer  et  al..,  1990). 

Water  \'apoi-  retrieval  from  satellites  occurs  primarily  in  two  regions  of  the  electro¬ 
magnetic  spectrum,  the  infrared  and  niicrowave/millimeter  wave  regions.  The  HiRS/2 
component  of  the  polar  orbiting  TIROS-N  operational  vertical  sounder  (TOYS)  (Smith  et 
aL,  1979;  Ehllger  and  Vonder  Haar,  1981)  and  the  VISSR  atmospheric  sounder  (VAS)  on 
Geostationary  Operational  Environmental  Satellites  (GOES)  (Smith,  1983)  are  examples 
of  infrared  instruments.  These  instruments  use  the  absorption  in  the  6.7  fxm  infrared  vi¬ 
bration  hand.  They  have  had  limited  success  in  obtmning  operationally  useful  water  vapor 
profiles,  in  part  because  of  the  limitation  of  all  infrared  sensors  that  moisture  sounding  is 
impossible  at  levels  below  cloud  top. 

The  appeal  of  microwave/millimeter  wave  sounding  is  that  temperature  retrievals 
suffer  limited,  if  any,  degradation  due  to  the  presence  of  clouds.  For  years  the  concept 
of  using  microwave  remote  sensing  for  water  vapor  retrieval  was  discussed.  The  first 
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efforts  focused  on  the  use  of  the  weak  H2O  rotational  line  at  22.235  GHz  for  sensing  total 
precipitable  water.  The  Nimbus  5  (Nimbus  E)  Satellite  Microwave  Spectrometer  (NEMS) 
used  channels  at  22.235  and  31.4  GHz  to  measure  atmospheric  water  vapor  and  liquid 
water  over  the  oceans,  even  in  the  presence  of  many  types  of  clouds.  NEMS  also  measured 
the  atmospheric  temperature  profile  (as  mentioned  above)  from  0-20  km  (Staelin  et  aL, 

1976).  ' 

Millimeter  wavelengths  also  possess  the  advantage  of  low  emissivity  values  for  the 
ocean  s\uface  versus  infrared  wavelengths.  The  low  emissivity  values  provide  a  relatively 
cold  background  that  provides  enough  contrast  to  allow  the  retrieval  of  low  level  moisture  * 

fields  (Hoffinan  et  aL,  1990).  The  accuracy  of  these  low  level  retrieved  values  decreases 
with  high  surface  wind  speeds  due  to  the  formation  of  foam  on  the  ocean  surface.  This 
foam  causes  the  emissivity  of  the  ocean  to  increase,  thus  effectively  reducing  the  contrast  I 

of  this  background  surface.  AdditionaUy,  millimeter  wavelengths  are  able  to  sense  through 
clouds  so  we  can  measure  water  vapor  in  the  atmosphere  below  the  clouds. 

The  Special  Sensor  Microwave/Imaga*  (SSM/I)  which  has  been  carried  on  DMSP  ^ 

satellites  since  1987  is  also  used  to  retrieve  precipitable  water.  Jackson  (1992)  presents  a 
review  of  several  of  the  methods  used  for  retrieval  of  precipitable  water  with  the  SSM/I. 

The  SSM/I  also  makes  use  of  tbs  water  vapor  line  at  22.235  GHz. 

1.3  USE  OF  THE  183  GHz  ABSORPTION  LINE 

The  water  vapor  rotation  absorption  line  feature  at  183.31  GHz  is  about  20  times 
stronger  than  the  rotation  featiure  at  22.235  GHz  (Waters,  1976).  For  this  reason  some  * 

researchers  have  discussed,  in  both  abstract  and  definitive  terms,  its  potential  use  to  profile 
atmospheric  water  vapor  since  the  1970’s. 

Gaut  et  aL  (1975)  performed  a  study  of  a  variety  of  microwave  remote  sensing  sys-  ^ 

terns  to  satisfy  Air  Force  meteorolo^cal  data  requirements.  The  simulation  and  retrieval 
exercises  investigated  the  following  parameters:  (a)  total  integrated  water  vapor,  (b)  the 
vertical  distribution  of  water  vapor,  (c)  temperature  profile,  and  (d)  the  integrated  and  ^ 

vertical  distribution  of  doud  liquid  water.  While  emphasizing  channds  below  60  GHz, 
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higher  firequency  channels  (including  those  in  the  vicinity  of  183  GHz)  were  also  evalu¬ 
ated.  The  higher  frequency  channels  were  explored  because:  (a)  it  was  desired  to  remove 
the  effects  of  variable  siuface  emissivity  over  land  at  the  lower  frequencies  used  to  obtsun 
water  vapor  information,  (b)  it  was  hypothesized  that  the  higher  frequency  channels  could 
provide  information  on  cloud  vertical  structure,  and  (c)  higher  frequency  channels  were 
expected  to  be  more  sensitive  to  the  integrated  water  vapor  of  clouds  with  relatively  small 
liquid  water  content  (Isaacs,  1987). 

The  pioneering  work  in  the  theory  of  using  the  183  GHz  rotation  line  for  atmospheri  c 
remote  sensing  was  written  by  Schaerer  and  Wilheit  (1979).  Several  authors  (Kakar,  1983; 
Kakar  and  Lambrightsen,  1984;  and  Rosenkranz  et  al.,  1982)  have  since  discussed  various 
methods  for  the  retrieval  of  clear  sky  atmospheric  moisture  profiles  using  channels  around 
the  183  GHz  line. 

Subsequently,  the  theory  has  been  extended  to  include  clouds  in  the  instnunent  field  of 
view  (FOV).  Isaacs  and  Deblonde  (1987)  examined  some  of  the  effects  of  beam-filling  liquid 
water  clouds  on  millimeter  wave  moisture  retrievals.  They  compared  retrieved  soundings 
with  radiosondes  from  clear  sky  and  cloudy  conditions  in  simulation  studies.  WUheit 
(1990)  performed  simihir  simxilation  studies  expanding  on  the  work  he  did  previously  with 
Schaerer.  Both  of  these  simulation  studies  were  limited  by  the  constraint  of  only  one  cloud 
layer. 

Wang  et  oL  (1983)  showed  experimentally  that  the  183  GHz  line  could  be  used  for 
profiling  atmospheric  water  vapor.  They  used  an  instrument,  the  Advanced  Microwave 
Moisture  Soimder  (AMMS),  from  an  sdrbome  platform  and  this  work  was  for  clear  sky 
cases.  Lutz  et  oL  (1991)  used  the  algorithm  developed  by  Wilheit  for  profiling  atmospheric 
water  vapor  even  in  the  presence  of  clouds.  This  work  again  used  measurements  from  the 
AMMS.  They  showed  the  retrieved  profiles  were  in  general  agreement  with  those  from 
radiosonde  data.  Lutz  et  al.  also  noted  that  the  algorithm  did  not  perform  very  well  in 
the  vicinity  of  surface  fronts.  With  the  concept  proven  by  these  sdreraft  studies  the  next 
step  was  to  put  a  moisture  sounder  on  a  satellite. 


1.4  MOTIVATION,  OBJECTIVES,  AND  THESIS  OUTLINE 


1.4.1  Motivation 

The  attraction  for  working  with  SSM/T-2  data  is  the  opportunity  to  work  data  from 
the  newest  microwave  instrument,  also  the  first  with  frequencies  near  183  GHz,  on  an  Air 
Force  meteorological  satellite.  The  work  on  data  from  this  instrument  has  barely  begun 
and  the  ability  to  learn  about  some  of  the  c^abilities  of  the  SSM/T-2  is  a  special  oppor- 
trmity.  The  ability  to  retrieve  moisture  soundings  over  the  entire  globe  regularly  is  also 
appealing  to  the  needs  of  GEWBX  and  other  major  global  climate  research  programs. 
The  T-2  will  also  help  improve  the  output  of  numerical  weather  prediction  models.  Addi¬ 
tionally,  during  the  testing  of  the  capabilities  of  an  instrument  a  researcher  may  also  learn 
of  its  limitations,  whether  they  be  inherent  or  due  to  the  advancement  of  the  science. 

1.4.2  Objectives 

This  thesis  attempts  to  achieve  the  following  objectives: 

1.  Develop  an  imderstanding  of  many  of  the  sensitivity  factors  which  affect  the  mi¬ 
crowave  brightness  temperatures,  specifically  for  the  SSM/T-2  frequencies. 

2.  Relate  T-2  brightness  temperature  features  to  characteristics  of  the  general  circula¬ 
tion  of  the  atmosphere  and  characteristics  of  the  earth’s  surface. 

3.  Develop  a  simple  physical  retrieval  of  colunmar  water  vapor  which  may  be  used  to 
assist  in  the  interpretation  of  the  T-2  data  and  which  may  be  used  in  initializa¬ 
tions  of  more  complex  schemes.  This  retrieval  scheme  assumes  the  thermodynamic 
temperature  and  brightness  temperature  are  equal  at  a  level  of  approximately  equal 
water  vapor  overburden  (column  water  vapor  above  this  level)  irrespective  of  the 
temperature  and  moistiire  profiles. 

1.4.3  Thesis  Outline 

Chapter  2  overviews  the  characteristics  of  the  satellites  which  carry  the  SSM/T-2  and 
the  characteristics  of  the  instruments  flown  on  these  satellites.  Chapter  2  also  presents 
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brightness  temperature  maps  for  two  T-2  data  periods  and  an  analysis  of  these  data  is  pre¬ 
sented.  Chapter  2  concludes  with  scatter  plots  of  brightness  temperatures  from  two  regions 
which  demonstrate  distinctly  different  temperature  and  moisriire  properties.  Chapter  3 
develops  the  radiative  transfer  equation  as  it  applies  to  microwave  remote  sensing.  It  con¬ 
tinues  with  a  discussion  of  retrievals,  first  temperature  (briefly)  and  then  water  vapor  and 
introduces  the  concepts  of  weighting  functions.  Chapter  3  concludes  with  a  description  of 
various  models  (atmospheric  propagation,  atmospheric  temperature/pressure,  and  cloud) 
used  in  the  simulations  discussed  in  this  work.  Chapter  4  analyzes  the  response  of  the 
DMSP  SSM/T-2  channels  to  a  variety  of  factors  (what  is  called  the  forward  problem) 
including:  the  profile  of  atmospheric  moisture,  the  presence  of  clouds,  the  characteristics 
of  clouds,  and  the  bcickground  against  which  the  atmosphere  is  viewed.  The  analysis  of 
cloud  effects  here  is  limited  only  to  low  level  water  clouds.  The  possible  effects  of  ice 
scattering  at  183  GHz  is  not  considered  in  this  study  although  it  may  be  important  under 
some  circmnstances.  Chapter  5  introduces  the  concept  of  water  vapor  bmrden  and  through 
the  use  of  figures  and  tables,  it  is  shown  that  the  vapor  burden  is  (nearly)  constant  a  a 
given  frequency  near  the  183  GHz  absorption  line  and  is  only  slightly  dependent  on  the 
moisture  and  temperature  profile.  The  water  vapor  burden  is  mapped  onto  a  constant 
pressme  sruface  thus  providing  the  v/ater  vapor  burden  above  this  surface.  Chapter  6 
presents  conclusions  based  on  the  work  contained  herein  and  a  discussion  about  possible 
appUcations  of  the  SSM/T-2  instrument  and  the  possibilities  of  tmified  retrievals  using 
other  instruments  from  the  DMSP  microwave  suite. 


Chapter  2 


THE  DMSP  SSM/T-2  AND  T-1  INST.^IUMENTS  AND  DATA  ANALYSIS 
2.1  DMSP  Block  5D  Satellite  Characterisaf.s 

DMSP  Block  5D  Satellites  operate  in  a  sim  jyncbronous,  near-polar  orbit  with  a  I 

nominal  altitude  of  800  km.  The  inclination  angle  of  the  orbit  ia  98  .P^  and  an  orbit  period 
of  102.0  min.  The  orbit  produces  14.1  orbit  revolutions  of  the  satellite  per  day.  The 
combination  of  these  parameters  produces  data  void  areas,  diamond-shaped  ixi  appe?  lance  ^ 

near  the  equator,  in  a  24  hour  period.  This  missing  coverage  has  obvious  drawbacks  when 
the  data  are  to  be  applied  to  study  the  tropics  (such  as  in  the  use  of  SSM/a  to  monitor 
tropical  cyclone  intensity  and  position).  These  data  void  regions  shift  with  each  orbit  aurct 
complete  coverage  of  the  tropical  regions  is  achieved  after  72  hours.  The  orbit  inclination 
also  results  in  dradar  sectors  of  230  km  at  both  poles  that  are  never  sampled.  The 
swath  width  for  the  SSM/I  and  SSM/T  channels  is  1400  km  perpendicular  to  the  satellite 
subtrack.  The  SSM/I  is  a  conically  scanning  radiometer  (with  a  constant  scan  angle  of  45 
degrees  aft  of  the  satellite)  as  compared  to  the  SSM/T  radiometers  which  are  cross- track 
scanning  soimders  (which  scan  in  a  fashion  similar  to  the  OLS  described  below). 

The  primary  instnunent  on  the  DMSP  satellites  is  the  Operational  Linescan  System  ^ 

(OLS)  which  is  a  fom  channel  imager.  The  OLS  visible  channel  for  very  high  resolution  has 
a  resolution  of  600  m.  Since  the  bandwidth  of  this  channel  is  very  broad,  the  instrument 
receives  more  radiation  from  a  given  scene  (e.g.  as  compared  to  the  imagers  on  NOAA  I 

polar  orbiter  satellites)  and  hence  has  a  greater  sensitivity.  The  maximum  swath  width 
for  the  OLS  is  2500  km.  The  OLS  produces  the  visible  and  IR  im^ery  used  for  many 
operational  requirements  by  forecasters  in  the  Department  of  Defense  (DOD). 
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The  F8  DMSP  spacecraft,  launched  in  1987,  was  the  jfirst  satellite  with  the  Special 
Sensor  Microwave/Imager  (SSM/I).  The  SSM/1  is  a  seven-channel  passive  microwave  ra¬ 
diometer  which  operates  at  four  frequencies.  It  receives  both  vertically  and  horizontally 
linearized  radiation  at  19.3,  37.0,  and  85.5  GHz  and  vertical  only  at  22.2  GHz  (Hollinger 
et  cd.,  1990).  Data  fi^sm  this  instrument  are  used  for  the  determination  of  the  following 
environmental  parameters: 

1.  Ocean  surface  wind  speeds. 

2.  Ice  coverage,  age  and  extent. 

3.  Cloud  water  content. 

4.  Integrated  water  vapor  (predpitable  water). 

5.  Precipitation  over  water. 

6.  Soil  moistiure. 

7.  Land  surface  temperature. 

8.  Snow  water  content. 

9.  Cloud  amoimt. 

The  SSM/I  is  also  flown  on  Fll  together  with  the  SSM/T-2  and  SSM/T-1.  These 
instruments  being  on  the  same  satellite  offer  the  possibility  of  a  synergistic  (also  called  a 
unified)  approach  to  atmospheric  retrieval  where  deficiencies  of  one  instrument  perhaps 
can  be  diminished  or  removed  through  the  use  of  another  instrument. 

2.1.1  SSM/T-2  Channel  Characteristics 

The  Defease  Meteorological  Satellite  Program  (DMSP)  Special  Sensor  Microwave 
(also  called  Millimeter  Wave)  Water  Vapor  Soimder  (SSM/T-2)  currently  flies  on  the  F- 
11  satellite  which  wais  laimched  on  November  28, 1991.  The  SSM/T-2  possesses  5  channels: 
three  dual-pass  bands  located  on  the  183  GHz  (1.64  mm)  water  vapor  absorption  line,  one 
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on  the  line’s  ning  (150  GHz),  and  a  window  channel  (91.655  GHz).  The  local  descending 
nodal  crossing  time  of  F-11  is  approximately  0514.  The  T-2  is  a  total-power  radiometer 
which  measures  the  total  power  seen  by  the  antenna  over  the  signal  processing  integration 
time.  The  objective  of  the  SSM/T-2  is  to  provide  global  measiurements  of  the  vertical 
profile  of  water  vapor  to  support  Air  Force  applications  including  the  input  of  these  profiles 
to  numerical  weather  forecast  models.  The  characteristics  of  the  SSM/T-2  chaimels  are 
shown  in  Table  2.1  (including  the  noise  equivalent  temperature  uncertainty.  NEDT). 


Table  2.1:  SSM/T-2  Channel  Characteristics  (after  Griffin  et  cd,,  1993). 


Chan. 

No. 

Center 

Frequency 

(GHz) 

Nadir 

FOV 

(km) 

Beam 

Acceptance 

(degrees) 

Peak 

Altitude 

(hPa) 

Pol. 

NEDT 

(OK) 

Re.sponse 

HQH 

183.31  ±  3 

48 

3.3 

650 

■1 

water  vapor 

mM 

183.31  ±  1 

48 

3.3 

500 

D 

water  vapor 

183.31  ±  7 

48 

3.3 

800 

Bi 

0.6 

water  vapor 

ImM 

91.655 

88 

6.0 

Surface 

Bi 

0.6 

surface 

150.0 

54 

3.7 

1000 

0.6 

surface 

The  T-2  Weis  originally  designed  to  be  a  tour  channel  instrument  (minus  the  channel 
near  92  GHz).  The  work  of  Isaacs  and  Deblojide  (1985)  showed  that  retrievals  were 
not  very  accurate  under  some  conditions  (.specifically,  tropical  atmospheres)  without  the 
91.655  GHz  channel. 

The  T-2  views  both  an  internal  hot-ioad  target  {•'■'  300°K)  and  cosmic  background 
radiation  for  calibration  measurements.  The  T-2  also  underwent  a  calibration  study  led  by 
the  Geophysics  Directorate  of  the  Air  Force’s  Phillips  Laboratory.  This  calibration  study 
included  independent  measmements  €md  model  calculation  studies.  These  independent 
measurements  were  performed  by  the  NASA  Millimeter-wave  Imaging  Radiometer  (MIR) 
carried  on  the  NaSA  ER-2  aircraft  (Racette  et  al.,  1992).  The  MIR  and  T-2  contain 
essentially  the  same  chaimels  with  the  exception  of  the  window  channel  which  for  the  MIR 
is  located  at  89  GHz.  The  NASA  ER-2  equipped  vrith  the  MIR  undeflew  selected  DSMP 
satellite  passes  including  ocean,  coastal,  and  land  background  cases.  RMS  differences  of 
0.7-1.4°K  were  measured  between  the  T-2  and  the  MIR. 
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Radiative  transfer  model  calculations  >!srere  also  conducted  (at  coincident  times  and 
locations  with  the  model  of  Eyre  and  Woolf  (1988))  using  radiosonde  data  and  estimates 
of  the  STuface  emittances.  Although  the  comparisons  of  the  model  calculations  with  the 
T-2  data  showed  collocated  RMS  differences  of  6.5°K,  the  calibration  team  concluded 
based  on  the  data  they  had  obtained  that  the  SSM/T-2  suffers  no  significant  bias  in  its 
calibration  (Griffin  et  al.,  1993). 

2.1.2  T-2  Output  Parameters 

The  following  parameters  are  derived  firom  T-2  data  using  the  operational  retrieval 
retrieval  scheme  of  A.  Stogryn  (Boucher  et  aL,  1993): 

1.  Relative  hiunidity  at  1000,  850,  700,  500,  400,  and  300  mb  levels. 

2.  Specific  humidity  at  1000,  850,  700,  500,  400  and  300  mb  levels. 

3.  Water  vapor  mass  (kg/m^)  between  levels  surface-1000, 1000-850, 850-700, 700-500, 
500-300,  and  above  300  mb. 

Although  these  parameters  are  produced  operationally,  the  actual  retrieval  of  data 
on  water  vapor  from  T-2  data  is  complex  for  reasons  described  in  the  next  chapter  and 
there  renudns  considerable  scope  for  research  in  understanding  these  complexities. 

2.1.3  SSM/T-1  Chaimel  Characteristics 

The  Special  Sensor  Microwave/Temperatm-e  Soimder  (SSM/T-1)  was  developed  and 
built  by  Aerojet  Electrosystems  Company  for  the  U.S.  Air  Force  Space  Systems  Division 
Defense  Meteorological  Satellite  Program  (DMSP).  The  design  methodology  was  a  top- 
down  systems  development  approach.  The  SSM/T-1,  first  laimched  in  1978,  was  designed 
to  support  Air  Force  and  Department  of  Defense  operational  weather  forecasting,  specif¬ 
ically  the  data  required  by  the  Air  Force  Global  Weather  Central  (AFGWC)  numerical 
weather  prediction  (NWP)  models.  To  incorporate  the  T-1  data,  the  T-1  is  designed  to 
emulate  radiosonde  temperature  and  pressine  measurements,  i.e.  retrieve  temperatures  at 
the  standard  pressure  levels,  thicknesses  between  these  levels,  and  the  temperature  and 


pressure  of  the  tropopause.  Channel  selection  was  such  that  the  weighting  functions  are 
distributed  as  evenly  as  possible  throughout  the  atmosphere  (Falcone  and  Isaacs,  1987). 
Wdghting  functions  for  channels  1-4  of  the  T-1  were  shown  in  Figure  3.1.  The  channel 
characteristics  for  the  SSM/T-1  are  shown  in  Table  2.2. 


Table  2.2:  SSM/T-1  Channel  Characteristics. 


Center 

Frequency 

(GHz) 

Nadir 

FOV 

(km) 

Beam 

Acceptance 

(degrees) 

Peak 

Altitude 

(km) 

NEDT 

(»K) 

Response 

50.5 

204 

surface 

53.2 

204 

T  at  2  km 

54.35 

204 

T  at  6  km 

54.9 

204 

^Bf^B 

10 

T  at  10  km 

58.4 

204 

30 

T  at  30  km 

58.825 

204 

14 

16 

T  at  16  km 

59.4 

204 

14 

22 

T  at  22  km 

The  SSM/T-1  is  a  passive  ‘Dicke’  radiometer  which  scans  in  a  cross-track  fashion.  The 
radiometer  has  a  scan  time  of  32  seconds  and  a  dwell  time  for  each  scene  of  2.7  seconds;  the 
instrument  takes  0.3  seconds  between  each  scene  and  the  warm  and  cold  calibrations  occur 
during  the  remaining  time.  The  SSM/T-1  operates  in  an  oxygen  absorption  band  using 
seven  frequencies  between  50  and  60  GHz  and  is  similar  to  the  MSU  which  was  described 
above.  The  variation  of  the  dry  ^  absorption  at  each  of  these  frequencies  provides  the 
basis  for  the  retrieval  of  temperature  profiles  (Curtis  and  Shipley,  1987). 

2.1.4  SSM/T-2  and  T-1  Scan  Patterns 

The  scan  pattern  for  the  SSM/T-2  is  similar  to  the  SSM/T-1.  Since  the  SSM/T-2 
operates  at  a  longer  wavelength  than  the  SSM/T-1  and  hence  has  a  smaller  footprint, 
or  field  of  view,  the  SSM/T-2  samples  fom  times  as  many  scene  stations  (28  versus  7) 
per  scan  revolution,  and  miist  scan  four  times  as  fast  to  provide  coincident  samples.  The 
result  is  that  the  SSM/T-2  produces  16  times  as  many  total  stations  scenes.  The  footprint 
pattern  is  shown  in  Figure  2.1.  The  reason  for  coincident  samples  is  that  to  produce  water 


15 


( 


vapor  retrievals  the  T-2  operational  algorithm  also  requires  four  of  the  seven  measured  T- 
1  frequency  channels  (Channels  1-4)  to  provide  the  requisite  temperature  sounding.  The 
pattern  in  Figure  2.1  is  changed  by  the  presence  of  a  glare  obstructor  or  GLOB,  on  the 
satellite.  The  GLOB  is  designed  to  protect  the  OLS  from  solar  insolation  which  could 
degrade  the  instrument  but  this  also  means  that  1  scene  of  T-1  data  and  4  scenes  of  T-2 
data  on  the  sun  side  of  the  satellite  are  degraded  and  therefore  not  used.  * 

2.2  SSM/T-2  DATA 

The  SSM/T-2  data  used  for  this  thesis  was  kindly  provided  by  Mr.  Bruce  Thomas  ^ 

of  the  DMSP  Omaha  Field  Office  of  Aerospace  Corporation.  The  data  were  obttuned  for 
approximately  78  revolutions  of  satellite  F-11  in  March  1993.  The  data  package  received 
from  Aerospace  Corporation  consisted  of  T-2,  T-1,  and  SSM/I  data  in  a  VMS  back-up  save 

► 

set  format.  All  data  received  from  Aerospace  were  in  sensor  data  record  (SDR)  format 
(also  known  as  antenna  pattern  corrected  temperatmes). 

The  T-2  data  set  consists  of  13  words  in  each  record.  Word  1  contains  information  on 
the  scan  scene  (which  provides  the  scan  angle  for  the  sensor)  and  the  background  against  * 

which  the  sensor  is  viewing  (land,  ocean,  sea  ice,  or  coast).  The  reason  for  the  differen¬ 
tiation  between  backgroimds  is  to  have  a  mechanism  to  help  facilitate  the  determination 
of  the  surface  emission  for  each  viewing  scene.  A  discussion  of  emissions  for  various  back-  ) 

groimds  is  included  in  Chapter  4.  Word  2  cont^dns  the  Julian  horn  (referenced  to  00  UTC 
31  December  1967),  and  words  3  and  4  are  the  Julian  minute  and  second,  respectively. 

Word  5  is  the  longitude  (from  0.-360°  East)  and  word  6  is  the  latitude  (±  90°)  which  ^ 

together  determine  the  earth  location  of  the  scan  scene.  Word  7  provides  the  terrain 
height  (m)  and  word  8  is  the  lOOOhPa  height  (m)  from  the  operational  forecasting  model 
at  AFGWC.  Words  9-13  are  the  sensor  data  records  (SDRs)  for  the  five  T-2  channels 

» 

(Kelvin)  (B.  Thomas,  personal  communication). 

2.3  T-2  BRIGHTNESS  TEMPERATURE  MAPS 

The  T-2  brightness  temperatures  fields  for  the  available  data  from  March  1993  were  » 

mapped  and  contoured.  Since  the  4  scan  scenes  on  one  side  of  the  swath  are  tmusable 
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Figisr*  2.1:  Footprint  pattom  for  th*  SSH/T-2  and  T-1,  (D.  Rhudy,  personnel  conBmmication)  . 
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because  of  the  shielding  by  the  GLOB  as  mentioned  above  and  because  of  the  brightness 
temperature  changes  as  a  function  of  scan  scene  angle  (discussed  in  further  detail  in 

» 

Cluster  4),  the  outer  4  scenes  on  the  other  side  of  the  scan  were  also  removed.  This 

residts  in  a  scan  line  that  is  considerably  less  than  the  1400  km  swath  width  mentioned 

above  but  it  does  me^ul  a  data  set  which  is  more  representative  requiring  no  corrections 

to  the  T*2  brightness  temperatures  for  angle  effects.  The  result  of  this  is  that  there  are  ^ 

areas  during  the  two  periods  which  are  not  sampled  and  these  areas  are  represented  by 

gray  in  the  colored  diagrams. 

Data  were  binned  into  1  degree  latitude  x  longitude  bins  and  then  averaged.  Since  I 

resolution  at  nadir  is  about  50  km,  bins  twice  the  size  of  nadir  resolution  should  be 

representative  for  the  data. 

2.3.1  Channel  1  Brightness  Temperature  Maps  ^ 

Channel  1,  at  183.31  ±  3  GHz,  is  used  for  the  detemaination  of  mid-upper  tropospheric 
water  vapor  and  relative  humidity  information.  As  shown  below  (in  Figiu^s  4.1-4.3), 
the  peak  of  the  weighting  function  for  this  channel  lies  between  3-0  km  for  the  three  ^ 

McClatchey  atmospheres.  Figures  2.2  and  2.3  are  the  Channel  1  brightness  temperatmes 
for  the  periods  10-14  March  and  16-18  March  1993,  respectively. 

The  areas  in  red  which  are  located  both  north  and  south  of  the  equator  (e.g.  roughly  I 

20-25°  from  the  equator)  delineate  the  location  of  subsidence  associated  with  the  sub¬ 
tropical  highs.  This  prominent  feature  is  also  prevalent  in  the  distributions  for  other 
frequencies.  The  subtropical  high  in  the  North  Atlantic  extends  from  western  Africa  ^ 

westward  all  the  way  to  the  Caribbean  Sea. 

Figure  2.3  suggests  a  more  coherent  and  stronger  ITCZ  in  the  Pacific  Ocean  (the 
white  and  green  color  slmding  with  several  purple  pixels).  The  South  Pacific  Convergence 

I 

Zone  (SPCZ)  (from  the  central  South  Pacific  extending  northwest  to  east  of  New  Guinea) 
is  also  more  distinct  in  this  figxire.  Both  features  appear  cold  (and  thus  moist)  in  contrast 
to  the  the  warm  (and  thxis  dry)  subsiding  tur  poleward  of  these  features. 
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2.3.2  Channel  2  Brightness  Temperature  Maps 

The  weighting  functions  of  Channel  2  peak  highest  in  the  atmosphere  (Figures  4. 1-4.3  i 

with  peaks  between  5-8  km)  and  measurements  provided  by  Channel  2  thus  applies  to 
upper  tropospheric  moisture.  Figure  2.4  and  2.5  are  the  Channel  2  brightness  temperattire 
maps  for  the  two  time  periods  described  above.  Maps  for  this  channel  have  less  spread 
of  brightness  temperatures  than  for  Channel  1.  Chaimel  2  shows  several  areas  with  high 
brightness  temperatures  over  the  oceans  and  these  are:  (1)  in  the  central  North  Pacific 
near  10-15°  latitude  and  near  the  date  line,  (2)  in  the  eastern  North  Pacific  south  of  Baja, 

California,  (3)  in  the  eastern  North  Atlantic  off  the  west  coast  of  Afiica,  (4)  in  the  South 

Atlantic  off  the  coast  of  Brazil  and  (5)  near  Madagascar  in  the  Indian  Ocean.  These 

brightness  temperatures  indicate  that  the  atmosphere  is  very  dry  in  the  upper  levels  and 

these  regions  are  also  over  the  subsidence  re^ons  indicated  in  Figures  2.2  and  2.3.  The  ^ 

coldest  brightness  temperatures  over  the  tropical  oceans  occm  near  Fiji  in  the  western 

Pacific  and  in  the  central  Indian  Ocean.  The  cold  brightness  temperatures  indicate  deep 

convection  in  these  areas  and  elevated  amounts  of  upper  tropospheric  moisture.  I 

2.3.3  Channel  3  Brightness  Temperature  Maps 

Channel  3,  located  at  183.31  ±  7  GHz,  senses  levels  of  water  vapor  that  occiu  deeper 
into  the  atmosphere  according  to  the  location  of  the  weighting  function  maxima  for  this  ^ 

channel  that  are  depicted  in  Chapter  4.  Figure  2.6  and  2.7  are  the  brightness  temperature 
maps  for  10-14  March  and  16-18  March,  respectively.  Once  again,  there  exists  coherence  in 
the  subsidence  regions  such  that  the  broad  areas  of  red  color  are  the  regions  of  subtropical  » 

high  pressures  imderlying  those  dry  regions  noted  in  reference  to  the  maps  of  Channels 
1  and  2  brightness  temperatme.  Similarly,  the  moist  tropical  regions  where  the  deepest 
convection  occmrs  appear  relatively  cold  (indicated  by  white-green  colors).  ^ 

2.3.4  Chaimel  4  Brightness  Temperature  Maps 

Channel  4,  the  window  channel  at  91.655  GHz,  primarily  senses  the  surface  of  the 
earth.  Figures  2.8  and  2.9  are  show  the  brightness  temperattires  from  Channel  4  for  the  * 

time  periods  described  above.  The  red  color  across  the  land  masses  (especially  Australia) 
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in  the  southern  hemisphere  point  to  clear  skies  (at  least  predominantly),  a  warm  smface 
and  a- high  emissivity  for  these  land  masses.  This  is  not  surprising  when  considering  that 
the  land  has  had  the  entire  summer  to  become  warm  and  dry.  The  lower  emissivity  of 
the  adjacent  ocean  surfaces  results  in  much  lower  brightness  temperatures  and  it  is  easy 
to  differentiate  between  land  masses  and  ocean  surfaces. 

East  of  Japan  there  are  four  curved  features  with  brightness  temperatures  warmer 
than  the  siurounding  region.  Daily  maps  of  brightness  temperatures  (not  shown)  indicate 
this  feattire  is  a  mesoscale  storm  system  which  has  been  viewed  by  the  T-2  on  successive 

I 

satellite  passes  (over  several  days).  The  doud  features  associated  with  this  low-presstue 
system  show  up  as  w<uin  brightness  temperatures  over  the  cold  ocean  background,  resem¬ 
bling  a  comma  cloud. 

Regions  of  the  tropical  oceans  show  little  structure.  The  yellow  and  brown  colors  * 

extend  farther  south  in  the  Southern  Hemisphere  as  compared  to  the  northern  extent  in 
the  Northern  Hemisphere  due  to  warmer  sea  su'iace  temperatures  and  low  level  moisture 
over  these  oceans,  since  the  weighting  functions  corresponding  to  this  channel  peak  at  /near  ^ 

the  surface.  Since  SSTs  typically  lag  by  about  two  months,  SSTs  in  the  southern  oceans 
should  be  just  beginning  their  fall  from  the  yearly  temperature  maximum. 

2.3.5  Channel  5  Brightness  Temperature  Maps  I 

Channel  5,  at  150  GHz,  is  located  on  the  line  wing  of  the  absorption  line  at  183.31 
GHz.  As  shown  in  Chapter  4,  the  weighting  function  for  this  channel  peaks  in  the  lower 
troposphere  or  at  the  earth’s  smface  (depending  on  the  moisture  profile  of  the  atmosphere).  I 

Figures  2.10  and  2.11  are  the  brightness  temperature  maps  for  the  periods  10-14  and  16-18 
March  1993.  The  maps  hint  at  the  ITCZ  and  associated  convection  (with  the  light  yellow 
and  white  pixels)  in  the  tropical  oceans.  Otherwise,  there  is  not  a  great  deal  of  structure  ^ 

to  the  temperature  differences  that  hsisn’t  already  been  smd  in  the  previous  subsections. 

Figure  2.11  does  indicate  a  stronger  SPCZ  than  does  Figure  2.10.  Australia  has  high 
brightness  temperatures  over  the  entire  coimtry  likely  indicating  a  high  emissivity  for  the 
land  and  a  dry  atmosphere  above  it.  Agmu,  the  land  surfaces  are  easily  found  in  the 
brightness  temperature  maps  due  to  the  emissivity  contrasts  with  ocean  surfaces. 
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The  brightness  tempa*ature  maps  in  Figure  2.2-2.11  provide  an  opportunity  to  ex¬ 
amine  some  of  the  general  circulation  features  of  the  atmosphere  and  to  see  the  smface 

► 

characteristics  of  ocean  and  land  masses  at  the  frequencies  of  the  SSM/T-2. 

2.4  CHANNEL  COHERENCE 

A  helpful  way  of  analyzing  and  understanding  the  properties  of  the  T-2  channel  bright-  * 

ness  temperatures  is  to  consider  the  relationships  between  channels.  Here  we  consider  two 
re^pons  which  are  characterized  by  distinctly  different  large  scale  circulation  features  and 
thus  different  temperature  and  moisture  properties.  » 

The  two  areas  selected  are  located  over  the  Pacific  Ocean.  The  first  area  corresponds 
approximately  to  the  TOGA  COARE  (Coupled  Ocean- Atmosphere  Response  Ebcperiment) 
dommn  and  is  defined  as  the  area  160-170°  East  longitude  and  0-10°  South  latitude.  This  ^ 

area  is  largely  imder  the  control  of  deep  convection  and  ascent  as  part  of  the  Hadley  and 
Walker  circulations  and  represents  a  region  of  relatively  high  moisture  content  at  all  levels. 

The  brightness  temperatures  from  this  area  are  compared  and  contrasted  with  an  area  in 

> 

the  Eastern  Pacific  defined  by  the  coordinates  110-120°  West  longitude  and  10-20°  North 

latitude.  This  E^t  Pacific  area  is  located  in  a  region  dominated  by  a  subtropical  high 

with  its  (assximed)  associated  strong  subsidence  and  lower  column  water  vapor.  Time 

constraints  did  not  permit  a  detidled  analysis  of  the  cloudiness  over  these  regions  during  ^ 

this  study. 

The  brightness  temperatiu-es  for  these  regions  and  for  the  5  channels  are  plotted 
against  each  other  in  the  form  of  scatter  diagrams.  Figure  2.12  shows  Channel  5  versus  * 

Channel  4  for  the  two  Pacific  regions.  The  East  Pacific  (EP)  region  shows  a  very  small 
scatter  for  the  brightness  temperatme  data.  The  Channel  4  temperatmres  are  relatively 
low  (as  compared  to  the  Channel  5)  since  this  channel  is  designed  to  see  the  surface.  ^ 

As  discussed  in  Chapter  4,  the  emissivity  of  the  ocean  surface  is  relatively  low  (~  0.7) 
and  it  provides  a  cold  background.  Channel  5,  on  the  line  wing,  provides  both  surface 
effects  and  the  effects  of  low  level  water  vapor.  The  brightness  temperature  indicates  an  ^ 

emission  temperature  near  280  K.  The  West  Pacific  (WP)  area  shows  a  window  brightness 
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Figure  2.2;  Channel  1  brightness  temperature  map  for  10-14  March  1993. 


226-0  234-0  242-0  250-0  258-0  266-0  274-0  282-0 


» 


n 

<f) 

o> 


CO 

0 


CO  CO 

0 

CD 

1— H 

L 

CD 

0 

0 

c 

> 

c 

CE 

m 

0 

CO 

X 

C 

-r-J 

o 

•H 

1 

QQ 

CD 

Csl 

1 

0 

jaMW* 

0 

X 

Ns 

L 

0 

z: 

0) 

L 

CO 

0 

0 

CO 

O 

ST 

-r-H 


I 


I 


I 


» 


» 


» 


> 


» 


» 


March  10-14 »  1993 

Channel  2  brightness  temperature  map  for  10-14  March  1993. 
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Figure  2.S:  Channel  3  brightness  temperature  map  for  10-14  March  1993. 
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Figure  2.7:  Channel  3  brightness  temperature  map  for  16-18  March  1993. 
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Figure  2.9:  Channel  4  brightness  temperature  map  for  16-18  March  1993. 
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Figure  2.10:  Channel  5  brightness  temperature  map  for  30-14  March  1993. 
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Figure  2.11:  Channel  5  brightne.ss  temperature  map  for  16-18  March  1993. 
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teaaperature  of  10-15  degrees  higher  than  the  East  Pacific  area.  This  likely  occurs  due 
to  more  emission  by  the  water  vapor  (near  the  surface,  emitting  at  a  high  temperature) 

► 

in  the  wetter  atmosphere  in  the  West  Pacific.  The  scatter  of  data  fi-om  the  West  Pacific 
indicates  that  the  area  contains  many  inhomogeneities.  Possible  reasons  for  this  scatter 
will  be  discussed  later. 

Figure  2.13  depicts  the  brightness  temperatures  of  the  150  GHz  channel  and  the  183  ^ 

db  7  GHz  channel.  Since  these  two  channels  are  adjacent  to  each  other,  it  is  expected  that 
the  two  channels  would  be  well  correlated.  For  the  Bast  Pacific  area,  the  Channel  3  and 
5  temperatures  are  nearly  equal,  with  most  of  the  Chmmel  3  temperatures  within  10°  of  ► 

280  K.  For  the  WP  region,  Chai  nel  3  becomes  opaque  more  quickly  and  therefore,  most 
of  the  temperatures  are  below  275  K.,  The  vertical  correlation  of  brightness  temperature 
(and  therefore  water)  shows  up  very  well  for  the  WP  region.  An  almost  linear  relationship  I 

between  these  two  channels  is  inferred  by  the  data. 

Figure  2.14  is  the  scatter  plot  of  the  brightness  temperatures  of  Channels  1  and  5. 

The  EP  area  again  exhibits  higher  brightness  temperatures  for  channels  which  peak  above  ^ 

the  surface,  impi3dng  a  dry  upper  atmosphere.  Additionally,  there  is  not  much  scatter  of 
the  data  for  this  region,  implying  a  profil'*  which  is  nearly  invariant.  Channel  1  brightness 
teniperatores  for  tlvs  WP  rsg'cn  are  about  20°  colder  than  the  EP  region.  This  indicates 

\ 

that  this  channel  is  sensing  water  vapor  at  a  colder  temperature  which  means  a  higher 
altitude. 

Perhaps  the  best  depiction  of  the  differences;  in  the  two  areas  is  shown  in  Figure 
2.15.  This  figure  plots  the  window  channel  (Channel  4)  against  the  channel  closest  to  the  ^ 

absorption  line  (Channel  2).  Figure  2.15  shows  that  the  EP  area  has  a  high  Channel  2 
brightness  temperature,  implying  the  channel  peaks  well  down  in  the  atmosphere.  The 
WP  area  has  noticeably  lower  brightness  temperatures  for  Channel  2,  indicating  a  peak  of  > 

the  channel  higher  in  the  vertical  (and  consequently,  at  a  lower  thermometric  temperature 
than  for  the  EP).  Channel  4  brightness  temperatures  are  governed  by  the  low  emissivity 
of  the  cold  ocean  surface  in  the  EP  and  the  higher  thermal  temperature  of  the  surface  and  ^ 

lower  atmosphere  combined  for  the  WP  region. 
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Figure  2.14:  Same  as  Figure  2.12  except  Channel  5  vs.  Channel  1. 
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Figure  2.16  also  shows  a  distinct  difference  in  the  brightness  temperatures  for  the  two 
Pacific  areas.  This  figure  depicts  the  brightness  temperatures  for  Channel  4  versus  Channel 
1.  The  explanation  of  these  features  is  basically  the  same  as  the  previous  paragraph 
as  the  Channel  1,  mid-upper  tropospheric  peak  for  the  weighting  function,  brightness 
temperature  is  high  for  the  re^on  with  subsidence  and  is  low  for  the  region  with  an 
abundance  of  water  vapor.  Conversely,  the  Channel  4  brightness  temperatmes  are  high 
for  the  WP  and  and  relatively  low  for  the  £P.  SST  differences  certainly  play  a  part  in 
these  brightness  temperatme  differences,  with  higher  SSTs  (which  mefuis  higher  brightness 
temperatures)  for  the  WP.  Additionally,  water  vapor  near  the  surface  and  cloud  effects 
also  contribute  to  the  higher  brightness  temperatures  in  the  WP  region. 

The  effect  of  ice  and  mixed  phase  and  partial  cloud  cover  clouds  cannot  be  quanti¬ 
tatively  nor  qualitatively  examined  from  these  figures  or  from  the  brightness  temperature 
maps.  This  is  an  area  where  future  research  is  needed. 

In  this  section  it  has  been  shown  that: 

1.  Areas  with  distinctly  different  moisture  profiles  (such  as  regions  of  large  scale  subsi¬ 
dence  versus  moist  ascent)  can  be  uniquely  identified  through  the  use  of  brightness 
temperatme  plots. 

2.  Surface  characteristics  can  be  inferred  from  the  brightness  temperature  plots. 


► 


» 
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Figure  2.16:  Same  as  Figure  2.12  except  Channel  4  vs.  Channel  1. 
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chapter  3 


MICROWAVE  RADIATIVE  TRANSFER  AND  WEIGHTING 

FUNCTIONS 

Emission  from  the  atmosphere  and  the  Earth’s  surface  provide  the  source  of  radia¬ 
tion  received  at  satellite  altitude  by  microwave  instruments.  Radiation  detected  by  these 
instruments  is  changed  by  the  processes  of  absorption,  emission  and  scattering  within 
the  atmosphere.  These  processes  depend  on  the  properties  of  the  Earth’s  surface,  atmo¬ 
spheric  constituents  especially  water  vapor  and  the  concentration  of  hydrometeors  and 
information  about  these  properties  can  be  retrieved  from  measurements  of  this  radiation 
(Stephens,  1993).  As  stated,  detection  of  atmospheric  water  vapor,  and  specifically  its 
vertical  profile  is  the  objective  of  the  SSM/T-2  Millimeter  Wave  Moisture  Sounder. 

For  the  frequencies  of  interest,  clear  sky  absorption  arises  primatrily  from  water  vapor 
and  to  a  lesser  extent  oxygen  absorption  (Waters,  1976.)  Ozone  also  absorbs  in  this 
region,  but  to  a  much  lesser  extent  than  water  vapor  and  oxygen,  and  ozone  absorption 
has  a  negligible  effect  on  brightness  temperatures.  Cloud  droplets  also  absorb  in  the 
microwave/millimeter  wave  region. 

Figure  3.1  provides  a  physical  perspective  on  microwave  emisb..>n  and  shows  some  of 
the  effects  of  water  vapor,  surface  reflectance,  and  cloud  liquid  water  on  this  emission. 
Case  (2)  shows  the  emission  spectra  found  over  a  cold  surface  (e.g.  ocean)  versus  case  (4) 
which  presents  the  spectra  foimd  over  a  warm  surface  (e.g.  land).  The  difference  in  spectra 
between  case  1  and  case  2  shows  how  water  vapor  changes  the  brightness  temperatures. 
These  effects  lead  to  temperatme  decreases  over  land  (warm  background)and  increases 
over  ocean  (cold  background).  Similarly,  the  change  in  spectra  between  cases  2  and  3  show 
the  effects  of  cloud  droplet  (Rayleigh)  absorption.  Finally,  the  spectra  change  from  case 
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2  to  case  4  is  an  example  of  how  surface  characteristics  may  change  intensities  measured 
by  a  satellite  radiometer. 

3.1  THE  RADIATIVE  TRANSFER  EQUATION  FOR  MICROWAVE  RE¬ 
MOTE  SENSING 


The  radiative  transfer  equation  describes  how,  in  a  mathematical  sense,  radiation 
is  transferred  from  one  layer  to  another.  The  generalized  form  of  this  equation  which 
includes  the  effects  of  both  absorption  and  emission  of  energy  is  written  as: 

dl  =  -k^[I  -  B]ds,  (3.1) 

where  I  is  the  intensity,  ki,  is  the  volume  absorption  coefficient,  and  B  is  the  emission 
represented  by  Planck’s  function.  Equation  (3.1)  neglects  scattering,  which  is  common  for 
many  problems  in  microwave  radiative  transfer  and  it  also  assumes  the  condition  of  local 
thermodynamic  equilibrium.  The  principles  of  emission-based  sensing  are  used  for  many 
atmospheric  applications  including  the  remote  sensing  of  sea  surface  temperature  (SST), 
column  water  vapor,  cloud  liquid  water  content,  precipitation  rate,  temperature  profiles, 
and  clouds. 


3.1.1  Brightness  Temperatures 


At  the  wavelength  of  microwave/millimeter  wave  radiometers,  the  Rayleigh- Jeans 
distribution  applies  —  that  is,  radiant  energy  and  temperatme  are  synonymous.  This 
temperatme  is  referred  to  as  the  brightness  temperature,  and,  for  satellite  radiometers  is 
a  measme  of  the  upwelling  thermal  microwave  radiation  from  the  medium  (in  this  case  the 
atmosphere)  below.  Using  the  Rayleigh-Jeans  approximation,  assinning  specular  reflection 
and  for  a  nadir  viewing  radiometer,  this  brightness  temperature,  Tbiv),  can  be  expressed 


£is: 


r°° 

T  {v)  =  /  exp(-T(/i,  oo))7(i/,/i)r(/i)d/i  +  exp(-T(0,  oo))  X 
Jo 


(1  -  R)T,fc  +  Rr  exp(-T(0,  h))^{u,  h)T{h)dh 
Jo 


-t-  Rexp(-2r(0,  oo))Tcb,  (3.2) 


» 


► 


> 
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► 
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DMSP  INSTRUMENTS 


FREQUENCY  (GHz) 


Figure  3.1:  Brightness  temperature  spectra  from  1-200  GHz  for  cases  (1)  tropical  atmo¬ 
sphere  with  no  water  vapor  and  reflectivity  (R)  =  0.3;  (2)  tropical  atmosphere  (CWV  = 
41.2g/m^)  and  R  =  0.3;  (3)  tropical  atmosphere  with  a  cloud  (LWC  =  0.5g/m*)  from  1-2 
km  and  R  =  0.3  and  case  (4)  same  as  case  (2)  except  R  =  0.0.  The  locations  of  the  SSM/I 
and  S'SM/T-2  channels  are  indicated  by  “F  and  “2”,  respectively. 
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where  T{z,y)  =  “f{v,x)dx  represents  the  optical  depth  between  heights  z  and  y,  R 
is  the- reflectivity  of  the  surface  ((1-R)  is  the  emissi'wty  (e)  of  the  surface),  7(1/,  h)  the 
voltune  absorption  coefficient  of  the  atmosphere  at  a  height  h  for  a  frequency  v,  T(h)  the 
temperature  profile  of  the  atmosphere,  T,fc  the  thermometric  temperature  of  the  surface, 
and  Tc6  the  cosmic  backgroimd  (2.7°K)  temperature.  Each  term  on  the  right  hand  side  of 
(3.2)  has  a  related  physical  description.  The  first  term  is  the  contribution  to  the  brightness 
temperature  due  to  upward  emission  by  the  atmosphere.  The  second  term  is  the  ‘stirface 
term’  which  is  the  radiation  emitted  by  the  Earth’s  siirface  which  reaches  the  satellite. 
The  third  term  is  the  downwelling  radiation  of  the  atmosphere  which  is  reflected  by  the 
Earth’s  srirface  to  the  satellite.  The  last  term  is  the  cosmic  backgroimd  temperature  which 
is  reflected  by  the  earth’s  surface  to  the  satellite.  The  angular  dependence  of  quantities 
has  been  omitted  here  and  is  taken  to  be  understood. 

There  is  some  debate  as  to  the  necessity  of  the  cosmic  background  term.  Some 
researchers  state  that  this  term  can  be  ignored  whereas  Kakar  (1983)  argues  that  for  a 
channel  on  the  wing  of  the  absorption  line  (e.g.  such  as  the  150  GHz  channel  of  the  SSM/T- 
2)  that  this  term  would  add  approximate!}'  0.5°K  to  the  measured  brightness  temperature. 
The  term  is  retjuned  retained  in  the  radiative  transfer  model  introduced  below  and  used 
to  produce  the  results  of  subsequent  chapters. 

In  calculating  brightness  temperatures  from  (3.2),  the  state  of  the  polarization  of  ra¬ 
diation  has  not  been  considered.  The  m^  sources  of  polarization  arise  from  reflection  by 
the  surfaceand  by  scattering  from  irregular  particles  in  the  atmosphere.  The  latter  effects 
are  ignored  and  only  the  the  effects  of  polarization  due  to  surface  reflection  eire  consid¬ 
ered.  By  ignoring  any  cross  polarization  effects  within  the  atmosphere,  (3.2)  becomes  the 
radiative  transfer  equation  for  the  vertical  state  of  polarization  assuming  R  represents  the 
surface  reflectivity  for  vertical  polarization.  In  any  event,  the  effects  of  surface  reflectiv¬ 
ity  are  minimal  for  Channels  1  and  2,  and  thus  the  effects  of  polarization  are  of  little 
consequence  for  the  study  described  in  this  thesis. 
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3.2  WEIGHTING  FUNCTIONS 

For  atmospheric  temperattirc  profiling,  (3.1)  may  be  expressed  in  the  form: 

rf,(oo)=  r  K{u,h)T{h)dh+  A,  (3.3) 

Jq 

where  A  represents  terms  independent  of  the  temperature  profile  and  K(j/jh),  the  tempera¬ 
ture  weighting  function,  is  weakly  dependent  on  the  temperature  profile.  These  weighting 
functions  are  a  basis  for  selection  of  channels  for  temperature  profiling  instruments  and, 
in  some  cases,  a  basis  for  the  temperatme  profile  retrieval  itself. 

Smith  et  aL  (1979)  show  weighting  functions  for  the  TOYS  instrument  suite.  The 
MSU  operates  in  the  same  oxygen  absorption  band  as  the  SSM/T-1  instrument.  The 
depicted  weighting  functions  (shown  in  Figure  3.2)  give  a  demonstration  of  the  vertical 
resolution  of  each  of  the  channels.  The  weighting  functions  appear  in  the  integral  equations 
relating  temperature  to  the  measured  radiances  (as  in  (3.3)). 

3.2.1  Relative  Humidity-based  Weighting  Functions 

This  discussion  of  constituent  retrievals  follows  from  the  work  of  Schaerer  and 
Wilheit  (1979),  Constituent  retrievals  are  unfortunately  a  nonlinear  problem  because 
the  absorption  coefficient  appears  both  in  the  transmission  factor  (exp(~  /  ydx))  and  in 
the  radiation  factor  'f{h)T{H)  of  (3.2).  Therefore,  the  effect  of  'y{v,h)  at  one  height  is 
dependent  on  7(1/,  h)  at  aU  other  heights  and  so  the  problem  of  constituent  retrieval  is 
also  a  non-local  one. 

Similarly,  (3.1)  can  be  expressed  as  (Schaerer  and  Wilheit  (1979)): 

nioo)  =  r  G{h)j{h)dh  -f  B,  (3.4) 

Jo 

where  B  represents  terms  independent  of  7(h);  however,  no  insight  comparable  to  the  tem- 
peratme  profiling  example  evolves  since  G(h)  is  strongly  dependent  on  7(h).  Nonetheless, 
(3.4)  can  still  be  linearized.  A  small  change  in  the  brightness  temperatme  A  Tb{oo)  due 
to  a  small  chcinge  in  the  absorption  as  a  function  of  height,  i.e.  A  7(h),  can  be  expressed 
as: 

ATb(oo)=  /~G[h,7(h),T(/»)]A7(h)d/i. 

Jo 


(3.5) 
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After  some  manipulation  we  find  an  expression  for  G(h,  7(fi).  T(h)): 


G{hMh),T{h))  =  exp(-r(fi,c»))[r(/j)  - 


i?exp(-r(0,oo))exp(-r(0,/i))|T(fe)  -  Tbi{h)], 


(3.6) 


where  and  are  the  upwelling  and  downwelling  brightness  temperatures  at  height 
b.,  respectively.  The  weighting  function  with  respect  to  absolute  water  vapor,  p(h),  is: 

Mh) 


G,{h)  =  Gih) 

and  with  respect  to  relative  humidity,  BB(h): 


dp{h) 


(3.7) 


Gnnih)  - 


(3.8) 


where  ptat{h)  is  the  saturation  vapor  pressure  of  water  at  a  height  h.  In  this  way,  weighting 
functions  based  on  the  relative  humidity  profile  RH(h)  are  produced;  these  being  related 
to  brightness  temperature  changes  by 


A'rf,t=  r  GRH{h)ARH{h)dh. 
Jo 

3.2.2  Transmission-based  Weighting  Functions 


(3.9) 


Transmission-based  weighting  functions  assume  a  different  yet  comparable  form  and 
are  similar  to  those  used  for  temperature  retrievals.  Isaacs  and  Deblonde  (1987)  used  these 
weighting  functions  for  their  study  of  cloud  effects  on  183  GHz  retrievals.  These  weighting 
functions  are  also  the  same  as  those  used  for  infrared  measurement  techniques  where  both 
vertical  profiles  of  water  vapor  and  total  precipitable  water  are  derived  (Prabhakara  and 
Dalu,  1980). 

Define  T(r,p)  =  exp(-r/p)  as  the  transmittance  of  the  atmosphere.  Neglecting 
the  dependence  on  the  zenith  cmgle  (as  above),  the  derivative  of  the  transmittance  be¬ 
comes:  dT  =  exp(— T)dr  or  dT  =  exp(— r)  (—'fdk).  Using  the  Rayleigh-Jeans  limit,  (3.2) 
becomes: 

Tj,t(oo)  =  r(T„)r./c(l  -R)+  y'°T(r)dT(r)  -fT(r)  rJ^  T{T)dT{T)  +T^{To)RTd,iZ.10) 


» 


44 


where  7,,  is  the  optical  depth  fitom  the  satellite  to  the  stirface.  By  grouping  common  factors 

together,  Equation  3.10  can  be  simplified  to:  j 

n^ioo)  =  [r,ycl j  -R)  +  i2r(r„)T^]  T(r„)  +  [1  r(T)/2]  T(r)dT(r)  (3.11) 

where  there  is  now  only  one  integral  term. 

Defining  a  weighting  Aoction  W(2)  for  a  vertical  coordinate  2  as  W(2)  =  = 

T{2)‘^{2)  then  the  prenous  equation  can  be  written: 

2bT(^»)  =  ir('2^o)(l  -  R)  +Rr{2,)T^]Ti2o)  -^[1+T{2)R]  T{Z)W{2)d2.{Z.12)  ^ 

Weighting  functions  derived  in  this  manner  are  convenient  since  they  closely  resemble 
the  weighting  functions  often  used  for  temperature  retrievals.  Transmission  is  a  quantity 
easily  calculated  in  a  propagation  model,  and  these  models  permit  the  calculation  of  the  > 

derirative  of  transmission  with  height  and  hence,  the  weighting  functions  themselves. 

As  an  example  of  transmission-based  temperature  weighting  functions,  Figure  3,2 
shows  weighting  fimctions  plotted  for  channels  1-4  of  the  SSM/T-1  temperature  sounder  for  ^ 

a  tropical  atmosphere  with  a  thermcmetric  temperature  of  300  K.  Figure  3.2  is  comparable 
to  Figure  2,  Liou  et  al,  (1979)  and  the  MSU  channels  of  Figure  i,  Smith  et  al.  (1979). 

Weighting  functions  from  the  two  different  methods  appear  very  similar  when  plotted. 

Plots  of  the  transmission-based  weighting  functions  are  for  the  SSM/T-2  are  discussed  later 
in  Chapter  4. 

3.3  ATMOSPHERIC  PROFILES,  PROPAGATION  MODEL,  AND  CLOUD 

MODELS  * 

To  simiilate  microwave  brightness  temperatme  requires  the  use  of  models  to  represent 
various  aspects  of  the  atmosphere  md  the  propagation  of  microwave  energy  through  the 
atmosphere.  We  use  simulations  of  the  brightness  temperatures  to  understand  how  they 
vary  with  changing  atmospheric  properties  and  to  study  these  effects  in  isolation.  The 
selection  of  the  models  described  below  was  made  for  a  variety  of  reasons  including  the 
use  of  a  model(s)  in  similar  studies  done  previously,  the  most  current  model  avsdlable,  the  ^ 

computational  efficiency  of  the  model,  and  its  simplicity. 


SSM/T-1  Weighting  Functions 


Weight  (dTr/dz) 


Figure  3.2:  Clear  sky  weighting  functions  for  the  SSM/T-1  channels  1-4  for  a  tropical 
atmosphere,  reflectivity  =  0.03,  and  a  colunm  water  vapor  (CWV)  =  41.2  kg/m^. 
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3.3.1  Atmospheric  Profiles 

McClatchey  et  al.  (1972)  standard  atmosphere  profiles  of  mid-latitude  winter,  * 

mid-latitude  summer,  and  tropical  atmospheres  were  med  as  the  reference  for  height  (km), 
pressure  (mb  or  hPa),  temperature  (“K  or  K),  density  (gm~®)  and  water  vapor  density 
(gm"^).  The  radiative  transfer  code  also  included  a  feature  which  would  allow  the  water  > 

vapor  density  to  be  increased  Ox*  decreased  by  a  common  factor  throughout  the  atmosphere. 

This  is  used  to  increase  the  munber  of  different  atmospheres  for  the  simulations  discussed 

in  Chapter  5.  ^ 

3.3.2  The  Propagation  Model 

The  Liebe  millimeter  wave  propagation  model  (MPM92)  (which  is  an  updated  version 
of  MPM89  (Liebe,  1989)  and  Liebe,  personal  communication)  allows  the  specification  of  an  > 

input  atmospheric  profile.  The  atmospheric  profile  provides  the  necessary  meteorological 
parameters  to  allow  the  computation  of  the  volume  absorption  coefficient,  7.  Knowledge 
of  the  volume  absorption  coefficient  then  allows  the  computation  of  the  optical  depth,  r.  I 

From  the  optical  depth  the  calculation  of  the  transmittance  is  strrdghtforward  and  the 
transmission-based  weighting  function  follows  as  described  above. 

The  MPM92  model  accoimts  for  dry  ^  attenuation,  oxygen  and  water  vapor  line  ^ 

absorption,  water  vapor  continuum,  and  Rayleigh  cloud  droplet  absorption.  The  Zeeman 
effect  for  oxygen  for  pressimes  less  than  0.7  kPa  and  Doppler  broadening  for  water  vapor 
are  also  included.  One  distinct  advantage  of  MPM92  is  the  small  number  of  lines  used 

I 

in  the  model.  The  model  uses  44  O2  Juid  30  H2O  local  line  contributions  versus  the  over 

500  contributing  spectrcd  lines  for  frequencies  below  1  THz.  The  small  number  of  lines 

in  MPM92  require  considerably  less  computation  time  than  do  other  codes  (i.e.  RAD- 

TRAN).  Versions  of  MPM  have  been  used  by  various  researchers  (e.g.  Wang  et  al.  (1992)  * 

and  Fleming  et  al.  (1991))  to  study  atmospheric  propagation  at  millimeter/microwave 

wavelengths. 


3.3.3  Cloud  Models 


To  simulate  cloud  effects  within  the  field  of  view  (FOV)  of  the  radiometer,  a  set  of 
cloud  models  was  chosen  to  be  incorporated  within  the  atmospheric  attenuation  profiles. 
The  cloud  models  chosen  (after  Isaacs  and  Deblonde,  1987)  are  from  the  Air  Force  Geo¬ 
physics  Laboratory  FASCODE  model  (Falcone  et  al.,  1979).  Cloud  parameters  for  the  five 
cloud  models  are  given  in  Table  3.1. 


Table  3.1:  Cloud  Type  Characteristics  (after  Isaacs  and  Deblonde,  1987) 


Cloud 

Model 

Type 

Liquid 

Water 

Content, 

gm-3 

Vertical 

Extent, 

km 

1 

Stratus 

2.7 

0.15 

0.5-2.0 

2 

Ciunulus 

6.0 

1.00 

1.0-3.5 

3 

Altostratus 

4.5 

0.40 

2.5-3.0 

4 

Stratocumulus 

6.25 

0.55 

0.5-1.0 

5 

Nimbostratus 

3.0 

0.61 

0.5-2.5 

Cloud  attenuation  was  calculated  using  the  Rayleigh  absorption  approximation  of 
Mie’s  forward  scattering  function  (Liebe,  1989).  Cloud  liquid  water  content  {w,  gm“®)  and 
frequency  are  the  required  variables.  Cloud  attenuation  (calculated  by  Liebe’s  MPM92 
model)  was  zidded  to  the  clear  sky  absorption  for  each  layer  of  the  atmosphere  for  the 
entire  cloud  layer  as  specified  by  the  cloud  model.  The  MPM  model  assumes  that  the 
cloud  fills  the  entire  FOV  and  that  there  is  only  one  doud  layer.  These  assiunptions  will 
tend  to  overestimate  and  underestimate  the  effects  of  douds,  respectivdy. 


Chapter  4 


SENSITIVITY  EFFECTS  FOR  THE  DMSP  SSM/T-2 

This  chapter  is  devoted  primarily  to  the  discussion  of  forward  radiative  transfer  sim¬ 
ulations  of  the  SSM/T-2  brightness  temperatures  and  how  these  temperatures  respond  to  I 

a  variety  of  different  atmospheric  parameters.  The  approach  is  to  build  an  imderstanding 
of  this  radiative  transfer  problem  before  tackling  the  inverse  problem  —  that  of  retrieving 
moisture  information  from  brightness  temperature  measurements.  There  are  a  nmnber  ^ 

of  factors  which  influence  brightness  temperatures  and  hence,  the  retrieval  of  moisture 
profiles  using  3  chaimels  near  183  GHz  and  a  window  channel  near  92  and  a  line  wing 
channel  at  150  GHz.  These  factors  will  be  referred  to  as  sensitivity  factors  and  include: 

► 

the  thermometric  temperature  and  its  vertical  profile,  the  smface  reflectivity  (emissivity), 
the  hiunidity  (moisture)  profile,  the  charactoistics  of  clouds  (both  water  clouds  and  per¬ 
haps  ice  clouds),  the  surface  elevation,  the  channel  noise,  and  scan  angle.  The  influence 
of  some  of  these  factors  on  T-2  simulated  brightness  temperatures  is  now  examined.  ^ 

4.1  DEPENDENCE  ON  MOISTURE  PROFILES 

For  the  simulation  studies  described  below,  three  basic  atmospheric  profiles  are  used:  I 

a  tropical  (T)  profile,  a  midlatitude  (MLS)  summer  profile,  and  a  midlatitude  winter 
(MLW)  profile  (McClatchey  et  aL,  1972).  The  water  vapor  density  data  for  each  profile 
will  also  be  multiplied  by  a  constant  factor  at  each  level  to  provide  a  method  of  varying  ^ 

the  coltmm  water  vapor  (CWV)  amount.  Since  each  moisture  profile  is  characterized 
by  a  different  colimm  water  vapor  amoimt  and  since  each  atmosphere  has  a  different 
temperature  profile,  the  simulated  channels  of  the  T-2  will  respond  differently  to  each 
profile. 
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The  clear  sky  weighting  hinction  profiles  for  an  ocean  background  with  an  emissivity 
of  0.7  are  shown  in  Figures  4. 1-4.3  with  altitude  as  the  ordinate  and  the  weighting  function  ^ 

^  (or  ^  as  designated  in  the  figures)  as  the  abscissa.  These  profiles  are  plotted  for  the 
nadir  angle  and  for  values  the  column  water  vapor  CWV  (kg  m~^)  as  given. 

Comparing  these  figures  it  can  easily  be  seen  how  the  amoimt  of  moisture  in  the 
atmosphere  affects  the  location  of  the  peaks  of  the  weighting  functions.  For  example,  for 
the  tropical  and  midlatitude  summer  atmospheres  the  182.31  GHz  channel  peaks  between 
7  and  8  km  whereas  for  the  midlatitude  winter  atmosphere  the  channel  peaks  between  5 
and  6  km.  Similar  effects  are  seen  for  the  other  channels  where  the  more  moisture  in  the  ^ 

atmosphere  the  higher  in  the  atmosphere  the  channel  peaks.  Even  the  ‘window’  channel 
(91.655  GHz)  sees  the  effects  of  moisture  in  the  atmosphere  and  this  serves  as  a  reminder 
that  the  ‘window’  is  not  entirely  transparent.  > 

It  can  also  be  seen  from  Figures  4.1-4.3  the  broadness  of  the  weighting  functions  as 
depicted  for  these  cases.  This  broadness  is  an  obvious  limitation  of  the  system  in  that  the 
resolution  for  retrievals  is  limited  by  the  lack  of  sharpness  of  these  weighting  functions.  ^ 

The  resolution  of  the  T-2  was  expected  to  be  on  the  order  of  2-3  km  based  on  the  weighting 
functions  of  this  type  (e.g.  Schaerer  and  Wilheit,  1979). 

Additionally,  from  Figures  4. 1-4.3  it  is  relatively  easy  to  develop  some  intuition  as 

► 

to  the  pmpose  of  each  channel,  at  least  for  dear-sky  conditions.  The  two  channels  dos- 

est  to  the  183  GHz  line  rdate  to  the  measmement  upper  tropospheric  moistme;  the  line 

wing  channel  (150  GHz)  to  low  levd  moisture;  and  the  window  channd  senses  low  levd 

moisture  and  siurface  characteiistics,  induding  the  emissivity  and  perhaps  the  thermomet-  ^ 

ric  (‘skin’)  temperature.  The  remaining  channd  is  rdevant  to  midtropospheric  moisture 

determinations . 

As  a  final  ex£imination  of  how  moisture  changes  the  wdghting  functions  of  the  SSM/T-  * 

2,  the  midlatitude  summer  profile  will  agdn  be  used.  In  this  case,  the  specific  humidity 
is  increased  by  a  factor  of  two  throughout  the  entire  atmosphere  (20  km).  Figure  4.4, 
compared  to  Figure  4.2,  shows  how  the  increased  moisture  at  all  levds  moisture  changes  I 

the  weighting  functions  significantly. 
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Weight  (dTr/dz) 


Figure  4.1:  Clear  sky  weighting  functions  for  the  SSM/T-2  channels  for  a  tropical  at¬ 
mosphere  over  an  ocean  backgroimd  with  a  reflectivity  factor  of  0.3  and  a  thermometric 
temperature  of  300  K. 
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Figure  4.2:  The  same  as  Figure  4.1  except  for  a  midlatitude  summer  atmosphere. 
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Figure  4.4:  The  same  as  Figure  4.2  except  with  2  x  water  vapor  profile. 
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These  four  figures  graphically  show  how  the  moisture  profiles  affect  the  weighting 
functions  when  it  is  plotted  against  height  (km).  Certainly  there  is  a  comprehensive 
change  in  the  weighting  functions  when  the  atmospheric  water  vapor  is  doubled.  This 
goes  back  to  the  discussion  earlier  about  the  nonlinearities  in  the  process  caused  by  the 
absorption  factor. 

Another  way  of  gr^hicaUy  showing  where  the  atmospheric  constituents  diminish  the 
emission  of  the  atmosphere  and  siuface  is  to  plot  transmission  versus  altitude.  Clear  sky 
atmosphere  plots  of  the  transmission  corresponding  to  Figures  4.1^.3  are  shown  in  Figures 
4.5-4.7.  Noteworthy  is  the  generally  similar  appearance  of  each  of  the  transmission  curves 
and  how  the  levels  of  the  most  rapid  decrease  in  transmission  with  decreasing  altitude 
correspond  to  the  peaks  of  the  weighting  functions. 

4.2  DEPENDENCE  ON  BACKGROUND 

The  surface  reflectivity  R  (or  emissivity  c  =  1-R)  depends  on  the  characteristics  of 
the  surface  of  the  earth  as  viewed  by  the  satellite  radiometer.  Land  and  coastal  regions 
have  large  variations  in  surface  emissivity  and  hence,  the  retrieval  of  moisture  profiles  can 
be  expected  to  sxiffer  degradation  for  these  regions  as  compared  with  ocean  surfaces. 

4.2.1  Ocean  Background 

The  emissivity  of  a  smooth  ocean  surface  can  be  calculated  from  the  dielectric  con¬ 
stant  of  water.  Calculations  using  the  dielectric  data  of  Ray  (1972)  result  in  a  value  of 
~0.7  for  the  emissivity  of  calm  ocean  water  at  183  GHz.  Lutz  et  al.  (1990)  have  devel¬ 
oped  an  algorithm  (from  Wilheit,  1990)  with  an  iterative  procedure  for  determining  the 
emissivity  for  wind-roughened  ocean  surfaces  or  land  surfaces.  The  effect  of  wind  rough¬ 
ening  procedure  degrades  the  ability  to  retrieve  moisture  profiles  (versus  retrievals  over  a 
constant  background). 

The  low  emissivity  of  ocean  surfaces  is  one  important  aspect  which  helps  with  the 
retrieval  of  low  level  moistme.  This  is  because  the  radiometer  senses  emission  from  the 
warm  atmosphere  in  contrast  to  a  cold  background  which  results  in  a  large  brightness 
temperature  contrast.  Additionally,  the  emissivity  of  a  water  surface  is  thought  to  be 
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Figure  4.5:  Clear  sky  transnussion  for  the  tropical  atmosphere. 
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Figure  4.6;  Clear  sky  transmission  for  the  midlatitude  atmosphere. 
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Figure  4.7:  Clear  sky  transmission  for  the  midlatitude  winter  atmosphere. 
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rather  stable,  especially  when  compared  to  land  smfaces.  One  problem  with  the  emissivity 
of  the  ocean  is  that  it  is  a  weak  function  of  the  surface  wind  speed  above  it.  Varioiis 
models  (e.g.  Wilheit,  1979  and  Pandey  and  Kakar,  1982)  were  developed  for  determining 
the  emissivity  of  a  wind-roughened  ocean  surface  for  earlier  instnunents  which  operated 
at  lower  fi%quencies  (e.g.  SMMR  and  SSM/I).  Cturently,  there  is  a  dearth  of  models  which 
are  applicable  for  the  region  near  183  GHz. 

To  exanxine  the  effect  of  the  emissivity  of  the  ocean’s  surface,  the  emissivity  of  an 
ocean  surface  \mderlying  a  clear  sky  nxidlatitude  winter  atmosphere  was  varied  between 
0.82  and  0.64.  These  results  are  shown  in  Table  4.1. 


Table  4.1:  Brightness  temperatures  (K)  for  emissivities  from  0.82  to  0.66  (reflectivities 
from  0.34  to  0.18)  for  a  midlatitude  winter  atmosphere. 


Midlatitude  Winter  Brightness  Temperattires 

FREQ. 

(GHz) 

Reflectivity 

0.34 

0.32 

0.30 

0.28 

0.26 

0.24 

0.22 

0.20 

0.18 

180.31 

254.99 

255.06 

255.13 

255.20 

255.27 

255.34 

255.41 

255.48 

255.55 

182.31 

245,77 

245.77 

245.77 

245.77 

245.77 

245.77 

245.77 

245.77 

245.77 

176.31 

243.01 

244.40 

245.80 

247.19 

248.59 

249.98 

251.37 

252.77 

254.16 

150.00 

203.81 

208.08 

212.35 

216.62 

220.90 

225.17 

229.44 

233.72 

237.99 

91.655 

193.83 

198.78 

203.74 

208.69 

213.64 

218.60 

223.55 

228.50 

233.46 

The  channel  closest  to  the  absorption  line  (182.31  GHz)  has  the  same  brightness 
temperature  for  all  emissivities  since  the  atmosphere  is  so  opaque  that  microwave  thermal 
energy  emitted  and  reflected  by  the  siuface  does  not  reach  the  satellite  (i.e.  the  smdace 
term  has  no  contribution).  The  180.31  GHz  channel  is  also  basically  imaffected  by  the 
surface  reflectivity.  The  line  wing  channel  (150  GHz)  and  the  window  channel  (91.655 
GHz)  brightness  temperatmes  change  appreciably  as  the  emissivity  changes  according  to 
the  values  depicted  in  the  table.  This  result  will  be  discussed  later  with  respect  to  the 
operational  algorithm  in  use. 

The  root  mean  square  (RMS)  error  caicidated  pre-launch  performance  of  the  T-2  with 
the  a  priori  regression  coefficients  iising  the  multiple  linear  regression  technique  employed 
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operationally  showed  the  best  retrieval  statistics  occurrs  over  ocean  backgrounds  with 
their  relatively  uniform  sinface  emissivity  (Boucher  et  aL,  1993). 

4.2.2  Land  Background 

The  emissivity  of  dry  land  backgrounds  is  thought  to  vary  from  0.8  to  1.0  over  the 
range  of  frequencies  of  140  to  200  GHz  (Schaerer  and  Wilheit,  1979).  These  early  values 
should  be  compared  with  those  of  Wang  et  al.  (1992)  below.  As  mentioned  above,  the 
Wilheit  (1990)  algorithm,  as  amended  by  Lutz  et  al.,  retrieves  the  land  emissivity  by  an 
iterative  process  as  part  of  its  computations. 

The  presence  of  vegetation  cover,  however,  can  change  the  reflectivity  of  land  surfaces 
as  vegetation  generally  is  an  absorptive  medium  and  therefore  decreases  reflectivity  values. 
Frozen  soils  also  generally  have  low  reflectivity  values  (~  0.1).  Soils  with  moisture  contents 
of  less  than  50%  have  dielectric  constants  where  the  real  part  is  less  than  6  which  yields 
a  value  of  R  of  less  than  0.15  (Wang  and  Schmugge,  1990). 

Snow  cover  presents  a  different  challenge.  Wang  et  al.  (1992),  using  the  AMMS 
described  in  Chapter  1,  (.served  emissivities  as  low  as  0.58  when  making  flights  over 
Alaska.  They  sxumised  that  the  reason  for  such  low  emissivities  is  due  to  scattering  by 
snow  particles  on  the  grotmd.  Wang  et  aL  found  the  emissivities  on  these  flights  generally 
r2inged  from  0.6-O.9  which  differ  from  the  estimates  of  Schaerer  and  Wilheit  of  more  than 
a  decade  ago. 

Conceptually,  low  level  moisture  retrievals  would  be  expected  to  be  more  difficult  over 
land  beccuise  of  the  high  and  variable  emissivity  of  this  surface.  This  presents  the  possible 
situation  where  emission  from  the  lower  level(s)  of  the  atmosphere  and  the  emission  from 
earth’s  smface  may  be  nearly  at 

the  same  brightness  tempei-aliu-e.  Isaacs  and  Deblonde  (1985)  noted  that  surface 
emissivity  affects  the  retrieval  accuracy  at  all.  levels  and  not  just  at  the  layers  near  the 
surface  based  on  the  statistical  retrieval  method  they  used. 

Brightness  temperatures  for  an  ocean  surface  with  an  emissivity  of  0.7  and  for  a  land 
surface  with  an  emissivity  of  1,0  are  compared  in  Table  4.2. 
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The  three  channels  nearest  the  absorption  line  (182.31  GHz  channel  not  shown)  are 
virtually  unaffected  by  the  surface  emissivity.  However,  this  is  certainly  not  the  case 
for  the  150  and  91.655  GHz  channels.  Radically  different  brightness  temperatures  are 
computed  for  cases  with  emissivities  of  1.0  and  0.7  for  these  channels.  This  is  simply  the 
factor  of  a  warm  background  and  a  cold  background,  respectively. 

As  mentioned  above,  for  the  ocean  background  case  in  Table  4.1  brightness  tempera- 
tTires  for  the  150  and  91.655  GHz  changed  appreciably  for  small  emissivity  changes.  The 
operational  algorithm  retrieves  the  s\uface  emissivity  using  the  91.655  GHz  channel  and 
a  physical  retrieval  model  (Griihn  et  al.,  1993). 

4.3  DEPENDENCE  ON  CLOUDS 

Water  clouds  affect  microwave  radiation  primarily  through  the  processes  of  absorption 
and  emission.  As  discussed  in  the  previous  chapter,  Rayleigh  absorption  is  assumed  by 
the  cloud  model  (MPM92)  being  used.  The  process  of  scattering  is  usually  considered 
negligible  since  the  single  scatter  albedo,  w,  at  the  microwave/millimeter  wavelengths  is 
very  small. 

To  depict  the  effects  of  clouds  on  weighting  functions.  Figures  4.8-4.12  show  the  five 
cloud  models  for  a  tropical  atmosphere  over  an  ocean  (R  =  0.3  or  e  =  0.7)  backgroimd.  For 
the  cloud  models  used,  and  from  the  study  of  the  figures  shown,  the  following  conclusions 
are  deduced: 

1.  The  channel  closest  to  the  absorption  line  does  not  ‘see’  any  of  the  low  level  water 
clouds. 

2.  The  180.31  GHz  channel  sees  only  those  douds  with  cloud  tops  above  2  km. 

3.  The  weighting  functions  for  the  other  annels  are  radically  affected  by  clouds 

in  the  FOV. 

4.  The  weighting  fimctions  peak  sharply  at  the  doud  top. 


61 


I 


Weight  (clTr/dz) 


Figure  4.8:  Weighting  functions  as  in  Figure  4.1  except  for  stratus  cloud  (w  =  O.lSgm  ®) 
from  0.5-2.0  km  (doud  model  1). 
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Figure  4.9:  The  same  as  Figure  4.8  except  for  oimulus  cloud  (u;  =  1.0  gm  from  1.0-3.5 
km  (cloud  model  2). 
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Figure  4.11:  The  same  as  Figure  4.8  except  for  stratocumulus  cloud  {w  =  0.55  gm~®)  from 
0.5- 1.0  km  (cloud  model  4). 
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4.3.1  Dependence  on  Layer  Thickness 

The  microwave  properties  of  clouds  in  the  Rayleigh  absorption  re^me  are  defined  by 
several  factors.  One  factor  is  the  thickness  of  the  cloud  layer.  This  .s  because  for  a  very 
thin  cloud,  much  of  the  microwave  energy  from  below  is  transmitted  through  the  cloud 
layer.  To  determine  the  effect  of  varying  cloud  thicknesses  on  the  brightness  temperatxires 
at  satellite  altitude,  clouds  were  inserted  into  the  tropical  atmosphere  with  a  constant 
cloud  base  of  1  km.  Cloud  thicknesses  were  varied  from  1  to  5  km  by  1  km  increments. 
Table  4.3  shows  simulated  brightness  temperatures  for  clouds  of  variable  thickness  and 
liquid  water  content.  According  to  the  simidations  (siunmarized  in  this  table),  the  150 
GHz  channel  and  to  a  lesser  extent  the  92  GHz  are  the  most  responsive  to  the  presence  of 
clouds.  The  responsiveness  in  the  150  GHz  channel  is  expected  since  it  is  designed  to  peak 
near  the  surface  as  depicted  in  Figure  4.1.  Note  that  brightness  temperatures  decrease  for 
increasing  LWC  and  increasing  cloud  thickness  (A  z)  with  a  land  background.  The  reason 
for  this  is  that  the  clouds  diminish  the  effects  of  the  smface  term. 

4.3.2  Dependence  on  Layer  Location 

It  has  been  seen  that  the  weighting  frmctions  for  atmospheres  containing  clouds  tend 
to  peak  near  the  cloud  top.  This  gives  a  good  indication  of  the  effects  of  where  a  cloud  is 
located  in  the  atmosphere  and  its  influence  on  brightness  temperatures. 

To  determine  the  influence  of  the  layer  location  on  brightness  temperature.  Table  4.4 
shows  the  simulated  brightness  temperatures  for  a  tropical  atmosphere  with  a  constant 
cloud  thickness,  reflectivity,  and  a  varying  cloud  base  and  liquid  water  content.  The  data 
in  this  table  indicates  how  the  brightness  temperatures  decrease  as  a  cloud  is  moved  higher 
into  the  atmosphere.  This  is  expected  since  the  temperatiu-e  profile  decreases  with  height 
and  the  cloud  limits  the  nucrowave  emissions  from  below  it. 

4.3.3  Dependence  on  Cloud  Liquid  Water  Content 

Figures  4.13-4.16  depict  brightness  temperatmes  for  the  five  T-2  channels  for  several 
cloud  models  with  the  LWC  varied  from  0.0  to  0.50  gm~®.  The  brightness  temperatures 
for  the  3  channels  near  183  GHz  are  basically  constant  whereas  for  the  150  and  91.655 
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Table  4.2:  Comparison  of  brightness  temperatures  for  land  (e  =  1.0)  and  ocean(e  =  0.7) 
stufaces. 


Brightness  Temperature 

Atmospheric 

case 

180.31  GHz 

176.3] 

GHz 

150.00  GHz 

91.655  GHz 

Land 

Ocean 

Land 

Ocean 

Land 

Ocean 

Land 

Ocean 

Clear  tropical 

264.71 

264.71 

277.58 

277.36 

291.19 

265.12 

295.95 

237.72 

263.54 

263.54 

276.32 

275.08 

290.91 

252.62 

295.84 

229.72 

Clear  ML  winter 

256.44 

255.31 

271.77 

249.34 

291.92 

223.19 

295.93 

216.26 

Cloud  model  1 

256.89 

256.89 

269.82 

269.81 

284.18 

272.76 

292.04 

248.02 

Cloud  model  2 

264.67 

264.67 

276.89 

276.87 

289.13 

277.67 

294.28 

253.08 

Cloud  model  3 

262.32 

262.32 

271.23 

271.23 

278.64 

278.51 

283.51 

279.32 

Cloud  model  4 

263.92 

263.92 

275.58 

275.52 

286.91 

273.42 

292.59 

250.31 

Cloud  model  5 

264.71 

264.71  i  277.22 

277.18 

289.55 

278.92 

294.77 

255.77 

Table  4.3:  Brightness  temperatures  for  variable  cloud  thickness  (A  z  (km))  and  LWC 
(gm"^)  for  a  fixed  cloud  base  (1  km),  emissivity  (0.97),  and  a  tropical  atmosphere. 


BR] 

[GHTNESS  TEMPERATURES  (K) 

Cloud 

A  z 

FREQUENCY  (GHz) 

LWC 

(km) 

180.31 

176.31 

150.00 

91.655 

Cloud 

1 

264.67 

276.92 

287.85 

289.91 

Model 

2 

264.08 

276.00 

286.20 

288.91 

1 

3 

263.11 

274.26 

283.54 

287.14 

LWC  = 

4 

261.99 

271.69 

280.13 

284.67 

0.15 

5 

260.36 

268.58 

276.24 

281.69 

Cloud 

1 

264.62 

275.24 

284.30 

288.59 

Model 

2 

263.42 

272.68 

280.90 

285.34 

2 

3 

261.25 

269.58 

275.91 

280.37 

LWC  = 

4 

258.78 

265.44 

269.86 

274.15 

1.00 

5 

255.71 

260.65 

263.89 

267.92 

Cloud 

1 

264.65 

276.50 

286.97 

289.59 

Model 

2 

263.90 

275.16 

284.43 

287.69 

3 

3 

262.63 

272.95 

280.49 

284.44 

LWC  = 

4 

261.17 

269.61 

275.41 

279.97 

0.40 

5 

259.13 

265.53 

269.92 

274.90 

» 
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Table  4.4:  Brightness  temperatures  for  variable  cloud  base  and  LWC  for  a  fixed  emissivity 
(c  =  1.0)  and  atmosphere  (tropical). 


BRIGH'] 

DNESS  TEMPERATURE  (K) 

Cloud 
Base  (km) 

LWC 

(gm~^) 

FREQUENCY  (GHz) 

182.31 

180.31 

176.31 

150.00 

91.655 

1 

0.15 

251.24 

264.67 

276.92 

289.34 

294.49 

1.5 

0.15 

251.24 

264.50 

276.57 

288.69 

294.02 

2 

0.15 

251.23 

264.08 

276.14 

287.99 

293.55 

WSB 

0.15 

251.11 

263.13 

274.65 

286.24 

294.43 

0.15 

249.51 

260.88 

271.16 

283.08 

290.25 

1 

1.00 

251.24 

264.62 

275.24 

284.54 

290.43 

1.5 

1.00 

251.24 

264.26 

273.99 

282.95 

288.81 

2 

1.00 

251.23 

263.42 

272.73 

281.44 

287.32 

3 

1.00 

251.04 

261.25 

269.71 

276.77 

283.11 

5 

1.00 

247.74 

255.82 

261.20 

265.86 

273.27 

GHz  channels  brightness  temperatures  vary  with  liquid  water  content,  cloud  location,  and 
cloud  thickness. 

Figures  4.15  and  4.16  are  identical  except  for  the  emissivity  (reflectivity).  The  latter  I 

figure  shows  the  effects  of  the  clouds  by  decreasing  the  91.655  and  150  GHz  brightness 
temperatures  over  a  300  K  surface  as  the  cloud  liquid  water  cont<ait  is  increased.  Contrast 
this  to  Figure  4.15  where  the  brightness  temperature  for  a  low  emissivity  surface  (cold  ^ 

background  at  210  K,  as  discussed  above).  In  this  case,  the  brightness  temperatures  of 
the  window  channels  systematically  increases  as  the  liquid  water  path  (LWP)  is  increased. 

This  yields  the  observation  that  clouds  generally  increase  brightness  temperatures 
over  oceans  and  decrease  brightness  temperatures  over  land.  These  effects,  however,  ap¬ 
pear  only  in  the  window  channels  and  to  a  lesser  extent  in  the  183  ±  7  GHZ  channel, 
depending  on  cloud  altitude.  This  response  to  clouds  can  also  be  seen  in  Table  4.2.  The 
near  linear  response  of  the  window  (especially)  and  line  wing  chaimel  (as  shown  in  Figure  ^ 

4.13)  to  changes  in  LWC  offers  a  simple  way  of  retrieving  this  cloud  property  by  using 
these  channels.  Several  authors  have  considered  this  (including  Isaacs  et  al.,  1985). 
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Figure  4.13:  Brightness  temperatures  for  Cloud  Model  1  and  a  tropical  atmosphere  with 
LWC  from  f*  ''t  to  0.50  gm~®. 
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Figure  4.14:  Brightness  temperatures  as  in  Figure  4.13  except  for  Cloud  Model  2. 
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Figure  4.15:  Brightness  temperatures  as  in  Figure  4.13  except  for  Cloud  Model  3. 
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Figure  4.16:  Brightness  temperatures  as  in  Figure  4.15  except  for  R  =  0.0. 


4.4  SCAN  ANGLE  EFFECTS 


The  effects  of  scan  angle  on  brightness  temperatures  can  be  considered  through  the  use 
of  MPM92.  Scan  angle  affects  brightness  temperatures  since  it  changes  the  optiad  path. 
As  shown  in  Figure  3.1,  the  scan  angle  (or  scene  station  view  angle  as  in  Figure  3.1)  on  the 
SSM/T-2  varies  between  1.5°  and  40.5°  as  the  sensor  scans  across  a  line.  Unfortunately, 
there  is  no  easy  way  to  correct  for  this  change  in  brightness  temperatme  for  the  change 
in  the  atmospheric  path  which  is  viewed.  Figures  4.17  and  4.18  show'  that  the  brightness 
temperature  can  both  increase  or  decrease  as  the  sensor  moves  away  from  mtdir.  There 
are  competing  effects  here;  limb  darkening,  occurs  when  the  sensor  moves  away  from  nadir 
and  the  path  length  increases  and  the  weighting  fimctions  move  upwards  vertically  in  the 
atmosphere.  Limb  brightening  occurs  for  the  wing  and  window  chaimels  which  sense  near 
the  surface.  The  temperature  of  a  cold  sruface  is  replaced  by  higher  temperatmes  as  the 
path  is  extended  by  viewing  away  from  nadir.  This  complex  problem  would  likely  be 
handled  operationally  through  the  use  of  regression  equations.  For  example,  regression 
equations  are  used  to  solve  the  scan  angle  differences  for  the  MSU  (Grody,  1983). 

4.5  INTERPRETATION  OF  RESULTS 

The  SSM/T-2  responds  to  a  variety  of  atmospheric  phenomenon  in  a  maimer  much 
different  than  those  instruments  which  preceded  it  in  the  remote  sensing  of  the  atmosphere. 
In  this  chapter  it  has  been  showm  that: 

1.  Weighting  functions  and  brightness  temperatures  are  dependent  on  the  moisture  and 
temperature  profile  of  the  atmosphere. 

2.  The  background  and  its  emissivity  have  a  significant  effect  on  the  brightness  tem¬ 
peratures  of  the  wrindow  channels. 

3.  Backgrounds  vnth  higher  emissivities  result  in  higher  brightness  temperatures. 

4.  Liquid  water  clouds  that  form  lew  in  the  at' nosphere  have  a  pnofoimd  influence  on 
brightness  temperatures  in  the  window  channels. 
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Figure  4.17:  Brigihtnesi  teiiq>«raturet  ta  a  function  of  obiervation  angle  for  a  tT(9)ical 
atmosphere. 
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Figure  4.18:  The  same  as  Figure  4.17  except  for  a  midlatitude  sixminer  atmosphere  and 
=  285  K. 
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(a)  Increasiiig  doud  UNC  decreases  brightness  temperatures. 

(b)  Higher  douds  correspond  to  lower  brightness  temperatures.  ^ 

(c)  Liquid  water  contcoit  affects  the  91  and  150  GHZ  channels  brightness  temper¬ 
atures  greatly. 

(d)  Clouds  gen*  rally  increase  brightness  temperatures  over  ocean  surfaces  and  de-  , 

crease  brightness  temperatures  over  land. 

Certainly  water  douds  affect  brightness  temperatures  in  a  systematic  manner.  Clouds 
become  more  difficult  to  quantify  when  one  considers  the  possibility  of  mixed  phase  douds,  i 

ice  douds,  and  partial  doud  cover. 

i 

t 

I 
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Ciiapter  5 

WATER  VAPOR  BURDEN 

In  this  chapter  a  simple  method  is  described  for  deriving  the  column  vapor  bur¬ 
den  from  brightness  temperatures.  It  will  be  demonstrated  how  the  measured  brightness  t 

temperatures,  especially  those  of  183  ±  i  GHz  (Channel  2)  and  183  ±  3  GHz  (Channel 
1),  correspond  to  a  level  of  constant  water  vapor  overburden.  This  result  together  with 
temperature-pressure  data  are  used  to  establish  the  burden  at  a  constant  pressure.  The  , 

method  is  demonstrated  using  only  Channel  2  brightness  temperatures. 

5.1  WATER  VAPOR  BURDEN  AND  THE  SSM/T-2  FREQUENCIES 

The  approach  followed  is  from  the  study  of  RosenJcranz  et  al.  (1982)  who  describe  the 
temperat\ue  weighting  functions  for  both  the  oxygen  band  (e.g.  SSM/T-1  frequencies)  and 
for  the  water  vapor  absorption  line  at  183.31  GHz.  They  define  the  integrated  colunm 
water  vapor  above  a  given  level  as  the  water  vapor  burden  (U)  which  is  proportional  * 

to  r,  the  optical  depth.  Rosenkranz  et  al.  adso  noted  that  when  weighting  functions 
are  expressed  as  functions  of  water  vapor  burden  rather  than  as  functions  of  pressure  or 
temperatme,  these  weighting  functions  are  nearly  symmetric  in  shape  and  their  peaJis  > 

are  invariant  with  respect  to  changes  in  the  water  vapor  profile.  These  functions  are  also 
displaced  downward  as  the  frequency  varies  away  from  the  absorption  line  et  183  (3Hz. 

Another  way  of  stating  this  principle  is  that  the  water  vapor  burden  may  be  expressed  ^ 

as: 

t/(h)=  rp{h')dh',  (5.1) 

Jh 

where  p  (h)  is  the  water  vapor  density  as  a  function  of  height,  h.  This  means  that  the  ^ 

satellite  brightness  temperature  is  a  weighted  average  of  the  temperature  profile.  The 
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weighting  is  defined  by  the  weighting  function,  which,  when  expressed  in  terms  of  U, 
is  approximately  invariant  with  respect  to  changes  in  p(h).  Therefore,  for  a  temperature 

> 

profile  with  (say)  a  linear  lapse  rate  near  the  peak  of  the  weighting  function,  the  brightness 
temperatiu^  and  thermometric  temperature  of  the  atmosphere  will  be  equal  at  the  peak 
of  the  weighting  function  (Wilheit  and  al  Khalaf,  1993).  This  leads  to  the  hypothesis  that 
the  water  vapor  burden  at  the  level  where  the  atmospheric  temperature  (Tt)  equals  the 
brightness  temperature  (Tj,)  will  be  constant  for  a  given  frequency  independent  of  the 
water  vapor  profile. 

To  test  this  hypothesis,  the  3  atmospheric  profiles  described  in  Chapter  3  were  used  for  ^ 

clear  sky  temperatme/pressure  profiles  and  various  water  vapor  profiles  (using  the  water 
vapor  multiplicative  factor)  were  employed  to  yield  47  cases.  The  water  vapor  burden  was 
calctdated  using  (5.1)  where  the  height  h  is  the  altitude  where  Tt  =  Tb.  Table  5.1  shows  I 

the  results  for  these  water  vapor  burdens  for  the  brightness  temperatures  at  183.31  ±  1 
GHz  for  these  47  different  atmospheric  cases. 

Nine  cases  resiilt  in  water  vapor  burdens  of  less  than  0.35  kg/m^  for  the  47  cases  in  ^ 

Table  5.1.  Certainly,  one  can  argue  whether  case  7  is  a  physically  realistic  atmosphere. 

The  other  cases  with  low  water  vapor  burdens  are  for  small  water  vapor  multiplicative 
factors.  Nonetheless,  the  hypothesis  appears  realistic  for  the  183.31  ±  1  GHz  channel.  For 
the  cases  considered  in  Table  5.1  a  mean  water  vapor  burden  of  0.3735  kg/m^  and  a  root 
mean  square  error  (RMSE)  value  of  0.0256  kg/m^. 

Water  vapor  burdens  were  also  calcidated  for  the  same  47  cases  for  183.31  ±  3  GHz. 

These  results  are  listed  in  Table  5.2.  Again  the  values  of  water  vapor  biirden  are  reasonably  ^ 

constant  with  a  mean  and  RMS  values  of  1.2739  kg/m*  and  0.3709  kg/m*,  respectively. 

The  cases  with  very  low  water  vapor  burdens  are  characterized  by  optical  depths  of  less 

than  4,  meaning  the  surface  emissivities  are  iofluencing  the  brightness  temperatures.  Com-  > 

puting  these  statistics  for  only  the  40  highest  cases  fields  a  mean  and  RMSE  of  1.4097 

and  0.1486  kg/m*,  respectively. 

In  Table  5.2,  case  8  has  a  s-urface  relative  humidity  of  about  14%  with  the  surface  I 

temperature  nem-  20°C.  Cases  16,  17,  and  18  have  relative  humidities  of  25-47%  with  a 
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Tkbk  6.1:  Water  V«^r  Overburdens  (U)  at  18S  ±  1  Ghs  for  Tb  =  Tf 


CASE 

HjO  MULT. 
FACTOR 

TEMPERATURE 

PROFILE 

WATER  VAPOR 
PROFILE 

U 

(kg/m“) 

cwv 

(kg/m’) 

1 

0.50 

Ttopical 

IVopical 

.4083 

20.61 

2 

0.76 

Trop>Al 

.3954 

30.91 

3 

1.0 

Tropical 

.3910 

41.21 

4 

1.25 

IVopicai 

.3806 

51.61 

5 

1.50 

Tropical 

.3639 

61.82 

6 

1.76 

liopical 

Tropical 

.3505 

72.12 

7 

2.00 

Tropical 

Topical 

.3389 

82.42 

3 

0.25 

ML  summer 

ML  summer 

.3962 

7.32 

9 

0.50 

ML  summer 

ML  summer 

.3844 

14.64 

0.76 

ML  summer 

ML  summer 

.3849 

21.95 

11 

1.0 

ML  summer 

ML  summer 

.3917 

29.27 

12 

1.25 

ML  summer 

ML  summer 

.3947 

36.59 

13 

1.50 

ML  summer 

ML  summer 

.3886 

43.91 

14 

1.75 

ML  summer 

ML  sununer 

.4057 

51.23 

IS 

2.00 

ML  summer 

ML  summer 

.3774 

58.54 

16 

0.40 

ML  winter 

ML  winter 

.3178 

3.42 

17 

0.50 

ML  winter 

ML  winter 

.3384 

4.27 

18 

0.75 
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stirface  temperature  of -1°C  and  cases  32-34  have  relative  humidities  of  5-11%.  Certainly, 
cases  32-34  would  be  extremely  rare,  a  surface  temperature  of  about  20°C  and  a  relative 
hiunidity  near  10%  (although  these  sturface  conditions  are  occasionally  observed  in  Col¬ 
orado,  but  the  surface  pressure  certainly  isn’t  1000  mb).  The  question  of  whether  these 
atmospheres  are  realistic  or  not  is  debatable  but  if  they  happen,  it  will  be  very  rarely. 

There  is  also  more  scatter  of  the  data  for  the  ±  3  GHz  channel  compared  to  the  ± 
1  GHz  channel.  The  results  in  Tables  5.1  and  5.2  correspond  well  with  results  of  similar 
studies  conducted  by  Wilheit  and  al  Khalaf  (1993). 

Figures  5.1  and  5.2  are  plots  of  the  water  vapor  burden  for  the  47  cases  for  183.31 
±  1  and  3  GHz,  respectively  (these  are  the  data  given  in  Tables  5.1  and  5.2).  Figure  5.1 
shows  how  the  vapor  burdens  corresponding  to  the  level  at  Tj,  =  Tt  only  slightly  varies 
from  atmosphere  to  atmosphere  whereas  Figure  5.2  shows  that  the  water  vapor  burden  is 
much  more  spread  although  the  outliers  ^ure  of  a  debatable,  physically,  as  stated  above. 

Figures  5.3-5.6  show  weighting  functions  plotted  in  terms  of  the  burden  U  for  the 
183.31  ±  1,  3,  7  GHz  and  150.00  GEbs  channels,  respectively.  These  figures  graphically 
portray  the  essence  of  the  h3rpothesis  stated  above.  While  the  magnitudes  of  the  weighting 
function  change  for  each  atmospheric  profile,  the  peaks  of  these  functions  and  their  shape 
are  remarkably  similar  for  these  different  atmospheric  profiles.  This  relative  independence 
of  the  shape  of  W(U)  on  both  the  temperature  and  moistme  profiles  provides  a  useful  way 
of  interpreting  the  SSM/T-2  brightness  temperature  distributions  such  as  those  shown  in 
Figures  2.2-2.11. 

5.2  CHANNEL  2  WATER  VAPOR  BURDEN  MAP 

Biised  on  the  discussion  and  results  of  the  previous  section,  this  section  presents  a 
retrieval  of  the  water  vapor  burden  for  Chaimel  2  (183.31  ±  1  GHz)  on  a  constant  pres¬ 
sure  stirface.  The  data  used  in  this  retrieval  consisted  of  the  T-2  brightness  temperatiues 
described  above  and  European  Centre  for  Medium  Range  Weather  Forecasts  (ECMWF) 
monthly  mean  temperature  and  pressure  analyses  for  March  1992  (used  to  simulate  the 
temperature  and  pressure  data  which  will  normally  be  available  from  the  SSM/T-1  Tem¬ 
perature  Sounder).  ECMWF  monthly  mean  data  were  used  because  temperature  data 
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Table  5.2:  Water  Wapor  Overburden  (U)  at  183.31  ±  3  GHz  for 
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Figure  5.1:  Water  Vapor  Overburdens  at  183.31  ±  1  GHz  for  47  Atmospheric  Cases. 
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Figure  5.2:  Water  Vapor  Overburden  as  in  Figure  5.1  except  at  180.31  GHz. 
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Figure  5.6:  The  same  as  Figure  5.3  except  for  150  GHz, 
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coincidenl  with  the  time  of  the  SSM/T-2  observatioiis  was  unavailable  at  the  time  of  this 
research.  The  results  presented  here,  especially  in  the  tropics,  where  synoptic  temperatxire 
variations  are  small,  will  not  greatly  alter  the  results  presented. 

The  retrieval  starts  from  the  assumption  that  the  brightness  temperature  and  the 
thermometric  temperature  of  the  atmosphere  are  equal  when  the  vapor  burden  is  0.3735 
kg/m^  (the  mean  value  from  Table  5.1).  Using  the  temperature  and  pressure  data  of 
ECMWF  fit>m  March  1992,  the  pressure  p*  corresponding  to  the  level  where  Tf  =  Tj,  is 
derived  by  interpolation.  Maps  of  p*  represent  the  contoxu^  of  the  constant  U  surface. 
High  values  of  p*  correspond  to  dry  conditicns  whereas  low  values  of  p*  correspond  to  a 
moist  atmosphere. 

A  plot  of  the  interpolated  p*  derived  in  this  way  is  presented  in  Figure  5.7.  The 
pressure  varies  from  high  values  in  the  mid-to-high  latitudes  (varying  from  600-500  mb) 
to  re^ons  assodated  with  deep  convection  in  the  equatorial  regions  that  have  p*  <  300 
mb.  These  areas  are  found  in  North-central  South  America,  the  West  Pacific  and  the 
South  Pacific  Convergence  Zone  (SPCZ),  Indonesia  extending  west  into  the  Indian  Ocean, 
and  portions  of  Central  Africa.  These  regions  may  be  considered  to  be  tbe  ascent  region 
of  the  Hadley  circulation  which  pumps  water  vapor  from  low  in  the  atmosphere  to  near 
the  tropopause.  The  subtropical  high  regions  are  not  eis  distinct  on  this  map  as  compared 
to  the  brightness  temperature  maps  although  some  regions  do  show  distinct  maxima  in  p* 
such  as  the  Bay  of  Bengal,  and  in  regions  across  the  northern  and  southern  Pacific.  This 
also  is  not  siirprising  since  the  upper  level  water  vapor  is  carried  from  the  equjitorial  region 
at  upper  levels  tow€irds  the  poles  so  this  water  vapor  at  upper  levels  would  be  expected 
to  only  drop  off  slowly.  Pressiures  increase  when  moving  away  from  the  equator.  This 
pressure  difference  between  the  tropics  and  the  arctic  regions  is  very  large  in  the  Southern 
Hemisphere  whereas  this  difference  is  not  nearly  as  great  in  the  Northern  Hemisphere. 

After  calculating  this  pressure  at  each  grid  point,  the  mean  pressure,  p*,  is  computed. 
A  simple  method  for  interpolating  the  water  vapor  information  contmned  in  Figure  5.7 
to  a  constant  pressure  surface  follows.  This  surface  is  selected  as  the  me€ui  pressure  p* 
between  40°  North  and  South  latitude  which  for  this  case  was  393  mb. 


Figure  5.7; 


87 


) 


The  water  vapor  burden  for  Channel  2,  Uj  has  a  value  of  0.37  kg  m  *  and  may  be 
expressed  as  the  integral; 

U2  =  -  !  (5-2) 

g  Jp* 

where  g  is  the  gravitational  constant  and  r  is  the  mi-iring  ratio.  Assuming  the  mixing  ratio 
can  be  expressed  by:  r  =  and  making  the  substitution  P  =  ^  where  po  is  the 

surface  pressure  and  r,  the  saturation  mixing  ratio,  (5.2)  can  be  rewritten  as: 

U2  =  ^  J^  p^dp.  (5.3) 

.  i 

Finally,  making  the  substitution  p*  =  |-  and  integrating  results  in  an  expression  for  U2: 

po 

The  difference  between  U  (p*)  and  U  (p*)  at  each  grid  point  is  now  determined  and  * 

used  to  adjust  the  column  overburden.  This  difference,  A  U  can  be  written  as  (from  (5.2)): 

Mf-^r  p^dp,  (5.5) 

9  Jp‘  > 

which,  on  integration,  becomes: 

(5-') 

On  rearrangement  using  (5.4)  for  U2,  (5.6)  can  be  written: 

(s-T) 

► 

Finally,  recalling  the  definitions  for  p  and  p  and  rearranging  results  in  an  equation  for  A  U 
which  is  added  to  the  value  of  U2  at  each  grid  point: 

(5-8)  , 

Since  p*,  the  pressiue  onto  which  U  is  to  be  mapped,  is  given  (393  mb),  p*  is  derived 
from  the  T-2  brightness  temperature,  and  U2  is  known  (0.37  kg/m*),  then  A  U  can  be 
readily  determined  and  the  burden  U(p*)  =  U  +  A  U  then  follows.  The  water  vapor  > 

overburden  U  at  p*  were  derived  in  this  way  and  mapped  (Figure  5.8).  The  water  vapor 
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burden  on  the  constant  piressure  surface  is  very  high  {~  0.7-0.8  kg/m^)  in  portions  of  the 
trc^ics  and  extremely  low  (below  0.2  kg/m^)  in  most  areas  north  of  Antarctica.  These 
results  correlate  well  with  the  previous  figure.  The  water  vapor  burden  at  ~400  mb  would 
be  expected  to  be  larger  at  the  tropics  than  towards  the  poles.  Stuprising  is  the  magnitude 
of  the  difierences  between  these  regions. 

This  process  is  instructive  when  one  considers  that  part  of  the  operational  output 
&om  the  T-2  is  layer  water  vapor,  as  stated  in  Chapter  2.  Performing  calculations  which 
yield  the  water  v^or  burden  on  several  levels  means  the  layer  water  vapor  mass  can  also 
be  carried  out  in  this  way.  Cert^nly  the  simplified  approach  shown  above  would  not  be  * 

used  operationally  but  it  does  provide  a  fi-amework  to  understand  how,  in  a  general  sense, 
the  layer  vary  In  the  atmosphere  and  may  provide  a  valuable  way  of  initializing  more 
sophisticated  retrieval  approaches.  ^ 
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Figure  5.8:  Water  vapor  burden  on  a  constant  pressure  (393  mb)  surface. 


Chapter  6 


CONCLUSIONS  AND  SUMMARY 


6.1  SUMMARY 

The  opportunity  to  be  among  the  first  researchers  to  work  with  data  firom  the  SSM/T- 
2  hficrowave  Water  Vapor  Sounder  presented  a  unique  opportunity  to  examine  some  of  the 
capabilities  of  this  instrument.  Since  the  SSM/T-2  is  a  new  instrument,  the  characteristics 
of  this  instrument  and  the  other  instruments  fiown  on  the  DMSP  spacecraft  F-11  were 
reviewed.  The  DMSP  microwave  suite,  (SSM/I,  SSM/T-1,  and  SSM/T-2  instruments), 
and  the  visible  and  infrared  scanners  result  in  a  multi-spectral  platform  from  which  obser¬ 
vations  of  the  earth  may  be  taken,  including  the  first  satellite  platform  with  frequencies 
near  the  183  GHz  absorption  line. 

Remote  sensing  in  the  microwave  regime  has  occurred  since  the  1960’s.  Techniques 
for  remote  soimding  in  the  microwave  region  have  been  proposed  since  the  1970’s  but 
only  since  the  laimch  of  the  SSM/T-2  has  it  become  possible  to  profile  the  distribution 
of  water  vapor.  Since  this  sounding  approach  is  solved  using  weighting  functions  the 
properties  of  these  functions  for  the  SSM/T-2  fi-equencies  are  examined  for  a  number  of 
different  circrimstances. 

Understanding  the  response  of  the  measured  SSM/T-2  brightness  temperatures  to  a 
variety  of  factors  (called  sensitivity  effects,  herein)  was  achieved  using  a  number  of  models, 
such  as  a  model  of  the  propagation  of  microwave  radiation  through  the  atmosphere,  models 
of  the  meteorolo^cal  state  variables  of  the  atmosphere,  and  simple  models  of  the  charac¬ 
teristics  of  clouds  in  the  atmosphere.  Examining  how  brightness  temperatures  respond  to 
different  atmospheric  parameters  provides  conceptual  framework  necessary  to  understand 
the  properties  of  the  SSM/T-2  brightness  temperatures  and  the  weighting  functions  that 
characterize  these  temperatures 


The  T-2  data  were  analyzed  by  averajpng  these  data  in  1°  latitude  y  ionptude  bins 
over  two  time  periods.  Because  of  the  small  swath  'jjidth,  the  result  of  discarding  data  &om 
the  4  scan  srenes  furthest  from  nadir  on  each  side,  some  regions  are  not  sampled  during 
the  two  time  periods  in  March  1993  (10-14  and  16-18).  The  global  maps  of  SSM/T-2 
brightness  temperabjres  derived  for  these  two  time  periods  were  examined  and  related  to 
the  gross  characteristics  of  the  general  circulation  of  the  atmosphere.  Additionally,  some 
information  ab'^ut  the  surface  characteristics  of  earth  and  the  moisture  profiles  above  the 
land  slufaces  are  discernible  from  these  brightness  temperatures. 

Comparisons  of  brightness  temperatures  from  two  regions  of  the  tropic^d  Pacific  were 
presented.  One  of  these  re^ons  was  located  within  the  TOGA  COARE  domain  and  is 
characterized  by  deep  convection  and  moist  conditions.  The  second  region  was  over  the 
Eastern  Pacific  imder  prevailing  subsidence  and  thus  represents  conditions  of  drying. 

The  concept  of  the  invariance  of  the  weighting  functions  invariant  with  respect  to 
water  vapor  overburden  is  tested.  It  is  demonstrated  that  the  level  where  Channel  2  and 
Channel  X  brightness  temperatures  match  the  thermodynamic  temperattne  are  levels  of 
^proximately  equal  overburden  irrespective  of  the  atmospheric  profiles  of  temperature 
and  moisture. 

The  pressure  p*,  which  coirespcnds  to  the  level  where  the  brightness  temperature 
for  Channel  2  equals  the  thermodynamic  temperature,  was  derived  by  interpolation  from 
temperatme-pressure  data.  The  water  vapor  burden,  for  Channel  2,  U2,  was  then  mapped 
on  a  constant  pressrrre  smface  p*  =  393  mb. 

6.2  CONCLUSIONS 

Prom  the  analyses  of  both  model  simulations  and  observations  performed  in  this  study, 
che  following  conclusions  were  made: 

1.  Weighting  functions  at  T-2  frequencies  (and  therefore  brightness  temperatures)  are 
strongly  dependent  on  the  profiles  of  water  vapor.  The  dependence  on  temperature 
is  weaker. 
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2.  Background  (strrface)  emissivities  strongly  affect  brightness  temperatures  at  91  GHz, 

150  GHz,  and  at  183  +  7  GHs.  Scenes  characterized  by  backgrounds  with  high 
emissivities  are  scenes  of  high  brightuesr  temperatures- 

3.  Liquid  water  clouds  profoundly  alfect  brightness  tempexatures: 

(a)  Increased  cloud  liquid  water  content  (LWC)  decreases  brightness  temperatures 
over  high  eniisd^dty  backgrounds. 

(b)  Cloudi!  generally  increase  brightness  temperatures  over  low  emissivity  (i.e. 

ocean)  surfaces  and  decrease  bri^tness  temperatures  over  high  (i.e.  land)  etuis-  * 

sivity  surfaces. 

4.  An  observational  study  of  the  T-2  brightness  temperature  maps  results  in  the  abil- 

) 

ity  to  see  some  characteristics  of  the  known  general  circulation  of  the  atmosphere. 

The  subtropical  re^ons  with  thdr  associated  subsidence  and  hence,  high  bright¬ 
ness  tempa'atures,  ai-e  striking  features  on  the  T-2  brightness  temperature  maps 
for  all  channels  on  the  183  GHz  absorption  line.  The  .South  Padfic  Convergence  ^ 

Zone  (SPCZ)  and  the  Intertropical  Convergence  Zone  (ITCZ)  appear  on  the  maps 
of  Channels  1  and  2  which  are  indicative  of  high  mid-to-upper  tropospheric  water 
vapor.  > 

5.  Surface  characteristics  can  also  be  observed  in  the  T-2  brightness  temperature  maps 
for  the  window  channels  (92  and  150  GHz).  Arid  land  areas  possess  large  emissivities 

and  high  brightness  temperatures.  All  the  land  masses  in  the  Southern  Hemisphere  ^ 

appear  warn!,  in  contrast  to  the  relatively  cold  adjacent  ocean  tenons. 

6.  Water  vapor  burdens  are  relatively  constant  for  the  two  SSM/T-2  frequencies  nearest 
the  absorption  line.  The  water  vapor  burden  is  0.37  ±  0.02  kg/m^  for  Channel  2 
and  1.4  ir  015  kg/m*  for  Channel  1.  The  pressure  of  the  sxirface  of  constant  vapor 
overburden,  introduced  here  as  p*,  varies  from  aroxmd  600  mb  in  high  latitudes  to 

pressures  less  than  300  mb  in  regions  associated  with  moist  convection.  * 


» 
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7.  Brightness  temperatisres  and  temperature-pressure  profiles  obtained  from  ECMWF 

analyses,  combined  with  the  above  result,  were  then  mapped  onto  a  constant  pressure  ^ 

surface.  This  provides  a  simple  method  for  retrieving  vapor  overbtirden.  The  method 
was  illustrated  for  a  constant  393  mb  pressure  surface. 

6.3  DISCUSSION  OF  FUTURE  WORK  , 

Obvioijsly,  since  the  SSM/T-2  represents  a  new  microwave  radiometer  on  a  satellite 
there  are  a  number  of  sensitivity  issues  which  still  need  to  be  examined.  The  first  issue 
that  comes  to  mind  is  that  both  mixed  phase  and  ice  clouds  are  not  discussed  V'^ry  much  * 

in  the  literature,  except  that  the  influence  of  ice  clouds  is  negligible  at  the  frequencies  of 
the  SSM/T-2.  j^Rxed-phase  clouds  also  represent  a  significant  hurdle  in  the  interpretation 
of  brightness  tooiperatures.  Certainly,  the  sdiility  to  model  this  problem  and  quantify  the  ^ 

effects  of  ice  and  mixed  phase  clouds  would  increase  the  ability  to  retrieve  water  vapor 
information. 

A  validation  test  of  the  retrieval  scheme  introduced  here  should  also  be  imdertaken. 

This  cotdd  be  accomplished  using  data  from  special  intensive  observation  periods  (e.g. 

TOGA  COARE)  where  water  vapor  information  derived  from  radiosonde  observations 
could  be  compared  with  coincident  T-2  data. 

» 

Improved  models  of  the  surface  emissivity  at  the  183  GHz  region  are  also  needed. 

These  models  should  address  the  enussivities  ocean  surfaces  and  sea  ice  and  snow  surfaces. 

The  possibility  of  using  the  T-2  for  snow  cover  mapping  should  also  be  investigated. 

The  183  GHz  re^on  is  under  study  by  NASA  in  a  high  priority  program  to  place  a  Mi-  * 

crowave  Imager  and  Soimder  on  a  future  geostationary  satellite  -  the  GEO-PLATFORM. 

A  microwave  sounder  of  this  t3rpe  would  certadnly  improved  the  sampling  capabilities  for 
mesoscale  forecasting.  I 

I  feel  the  immediate  future  (less  than  five  years)  will  see  an  enormous  amoimt  of 
research  along  the  lines  of  a  synergistic  approach/unified  retrieval.  The  scanners  and 
radiometers  on  DMSP  spacecraft  F-11  provide  a  data  set  which  is  as  comprehensive  as  ^ 

any.  Certainly  the  SSM/I  could  be  used  for  the  ocean  surface  wind  speeds  and  this 
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information  could  be  used  to  help  in  the  determination  of  the  background  emissivity  over 
ocean  surfaces  for  the  183  GHz  channels.  Additionally,  *'he  SSM/I  and  SSM/T-2  can  be 
used  together  for  column  water  vapor  and  cloud  liquid  water  content.  Here  it  may  be 
possible  to  explore  column  water  vapor  over  both  land  and  ocean  studaces.  Present  SSM/I 
analysis  of  this  column  quantity  is  restricted  to  ocean  areas  only  The  Adsible  and  infrared 
data  from  the  OLS  can  be  used  for  cloud  type  determination  (and  cloud  cover)  and  the 
responses  of  the  SSM/T-2  channels  can  be  evaluated  to  improve  the  understanding  of  the 
effects  of  clouds  in  the  183  GHz  region. 

The  SSM/T-2  will  serve  to  only  expand  the  research  opportunities  to  those  with  a  * 

keen  interest  in  water  vapor  and  it  will  undoubtedly  improve  the  ability  of  NWP  models  to 
forecast  the  state  of  the  atmosphere.  While  a  great  deal  of  work  remains  to  be  done  on  the 
use  of  the  183  GHz  absorption  line  for  water  vapor  profiling,  it  has  been  an  enlightening  I 

opportunity  to  use  the  T-2  data. 

> 

> 

> 
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ABSTRACT  OF  THESIS 

AN  ANALYSIS  OF  SIMULATED  AND  ACTUAL  DMSP  SSM/T-2  BRIGHTNESS 

TEMPERATURES 

Defense  Meteorolo^cal  Satellite  Program  (DMSP)  Special  Sensor  Microwave/T-2 
{SSM/T-2  or  T-2)  global  brightness  temperature  maps  are  related  to  some  of  the  known  ob¬ 
served  general  circulation  featmes  of  the  atmosphere  and  surface  characteristics  of  Earth. 
The  brightness  temperatures  from  two  re^ons  with  different  large  scale  circulations  are 
compared  and  contrasted. 

The  concept  of  the  invariance  of  the  weighting  functions  with  respect  to  water  va¬ 
por  overburden  is  tested.  It  is  demonstrated  that  the  levels  where  the  183  ±  1  and  3 
GHz  brightness  temperatures  match  the  thermodynamic  temperature  are  levels  of  ap¬ 
proximately  equal  overbrnden.  Water  vapor  overburden  (integrated  water  vapor  above  a 
pven  level)  was  retrieved  and  mapped  on  a  constant  pressure  surface  from  T-2  brightness 
temperatures  and  atmospheric  tempaature-pressure  data.  A  simple  physical  retrieval 
scheme  based  on  the  approximately  equal  overburden,  irrespective  of  the  moisture  profile, 
for  the  183  ±  1  GHz  frequency  was  developed  and  demonstrated.  This  simple  retrieval 
could  be  used  to  initialize  more  sophisticated  retrievals. 

Simtilation  studies  of  some  factors  which  affect  brightness  temperatures  at  the 
SSM/T-2  frequencies  were  conducted  and  the  results  discussed.  These  factors  included 
moisture  profiles,  backgrounds  with  different  emissivities,  and  low  level  water  clouds. 
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Chapter  1 


INTRODUCTION 


1.1  BACKGROUND 

Water  vapor  plays  a  vital  role  in  establishing  our  present-day  climate.  It  is  the 
most  significant  contributor  to  the  observed  greenhouse  effect  of  the  planet  (Stephens  and 
Tjemkes,  1992).  Water  vapor  links  the  processes  of  evaporation,  cloud  formation,  and 
precipitation.  Latent  heat  release  plays  a  large  role  in  the  dynamical  flow  fields  of  the 
tropical  atmosphere  and  knowledge  of  the  distribution,  transport,  and  divergence  of  water 
in  all  phases  is  critical  to  work  pertaining  to  the  difficult  problem  of  cloud  feedback  and 
its  role  in  climate  change  (Stephens,  personal  conununication). 

The  source  of  water  vapor  is  primarily  at  the  earth’s  surface  whereas  its  sink  is  in 
tbci  process  of  precipitation.  Combining  these  somces  and  sinks  with  the  atmospheric 
circulation  means  the  distribution  of  water  vapor  in  the  atmosphere  is  highly  variable 
both  in  spatial  and  temporal  scales  (Prabhakara  and  Dalu,  1980).  Thus,  it  is  necessary 
to  establish  a  dense  network  for  observing  water  vapor  and  its  distribution.  Present  day 
observations  are  largely  based  in  the  use  of  radiosonde  data  and  these  data  lack  in  the 
necessary  global  coverage  needed  for  tmderstanding  the  role  of  water  vapor  in  climate. 

Water  vapor  also  affects  the  exchange  of  radiation  between  the  atmosphere  and  space 
directly  by  its  influence  in  atmospheric  emission  and  indirectly  by  its  role  in  cloud  for¬ 
mation  processes  (Bates  and  Stephens,  1991).  GEWEX  (the  Global  Energy  and  Water 
Cycle  Experiment,  part  of  the  World  Climate  Research  Program)  has  the  basic  theme  of 
developing  a  quantitative  understanding  of  how  water  shapes  the  energy  budget  of  the  cli¬ 
mate  system  of  earth.  The  connections  and  relationships  between  the  global  hydrological 
cycle  and  the  energy  budget  are  critical  to  the  question  (and  potential  problem)  of  climate 


change  (Stephens  and  Tjemkes,  1992).  The  typical  climate  change  scenario  (2  x  CO2) 
predicts  a  global  warming  of  approximately  twice  (or  treble,  depending  on  the  model)  the 
C02-induced  warming  with  water  vapor  feedback  than  without  this  feedback  (Houghton 
et  aL,  1990).  The  first  part  of  the  water  vapor  feedback  is  due  to  the  direct  greenhouse 
effects  of  doubled  CO2  which  directly  leads  to  increased  sea  surface  temperature  (SST). 
The  increased  SST  increases  evaporation  in  the  boimdary  laj^er  which  increases  humidity 
in  the  boundary  layer.  The  increase  of  water  vapor  in  the  boundary  layer  through  its 
greenhouse  effect  increases  the  SST.  The  second  part  of  the  feedback  process  was  hy¬ 
pothesized  by  Lindzen  (1990).  Water  vapor  feedback  is  also  thought  to  be  modulated  by 
the  changes  in  the  vertical  distribution  of  water  vapor.  Lindzen  argues  that  the  increased 
warming  near  the  ground  results  in  increased  and  deeper  cumulus  convection.  This  in¬ 
creased  cumulus  convection,  he  claims,  leads  to  drying  of  the  upper  troposphere  (above 
5  km)  and  elevation  of  the  altitude  at  which  convected  heat  is  deposited.  Since  green¬ 
house  absorption  is  most  important  above  5  km  these  cumulus  effects  are  negative  and 
hence,  should  diminish  the  effect  of  CO2  warming.  In  contrast  to  this  view,  present  day 
climate  models  actually  produce  a  moistening  of  the  upper  troposphere  associated  with 
this  convection.  Both  hypotheses  are  depicted  in  Figure  1.1  (Stephens,  personal  commu¬ 
nication).  Attempts  to  study  quantitatively  the  distribution  and  transport  of  water  vapor 
in  the  atmosphere  have  been  made  using  the  global  radiosonde  network  (e.g.  Peixoto  and 
Oort,  1983).  Unfortunately,  these  studies  have  met  limited  success  since  these  radiosonde 
observations  are  confined  nudnly  to  land  areas  since  this  is  where  the  radiosonde  stations 
are  primarily  located.  Important  atmospheric  phenomenon  restricted  to  oceanic  areas  are 
thus  largely  ignored  (e.g.  many  aspects  of  El  Nino,  including  westerly  wind  bursts  and  the 
location  of  moisture  convergence/divergence).  Early  work  on  transports  based  on  in&ared 
(IR)  water  vapor  data  have  been  reported  by  Wittmeyer  (1990).  However,  Wittmeyer 
and  Vonder  Haar  (1993)  and  others  report  serious  adverse  cloud  impacts  in  baroclinic 
zones.  Schwartz  and  Doswell  (1991)  contend  that  more  observations  are  essential  in  both 
time  and  space  to  improve  mesoscale  weather  forecasting  and  further,  that  remote  sensing 
systems  have  yet  to  approach  the  radiosonde’s  capability  to  resolve  the  vertical  thermody¬ 
namic  structure  of  the  atmosphere.  While  the  later  may  or  may  not  remain  true  in  the  age 


Figure  1.1:  Hypotheses  on  water  vapor  feedback  mechanisms  (Stephens,  persona]  com¬ 
munication).  The  left  portion  of  the  dit^ram  is  the  feedback  as  presently  perceived,  with 
increasing  temperature  leading  to  increased  water  vapor  largely  in  the  boundary  layer.  The 
right  portion  of  the  diagram  is  the  additional  modificat  ion  to  the  feedback  by  increasing 
the  upper  tropospheric  moisture  through  convection. 


of  new  microwave  (millimeter  wave)  instruments,  the  former  may  be  improved  drastically 
by  these  new  remote  sensing  systems.  Certainly,  the  future  lies,  in  many  respects,  in  the 
advancement  of  the  capabilities  of  remote  sensing  and  the  application  of  these  capabilities 
to  the  measurement  of  atmospheric  phenomena. 

Satellites  are  used  for  many  applications  of  remote  sensing.  The  theory  and  pr€ictice 
of  remote  sensing  spans  many  disciplines  and  is  intimately  related  to  inversion  theory  (see 
Twomey,  1977).  The  deployment  of  satellites  expanded  the  opportunities  for  deriving  the 
distribution  of  water  vapor  over  the  enth-e  globe  (as  opposed  to  the  limitations  with  land- 
and  ship-based  rsidiosondes).  Moreover,  many  of  the  satellites  launched  the  1970s  and 
1980s  were  restricted  to  measurements  of  the  total  precipitable  water.  The  latmch  of  the 
Special  Sensor  Microwave/Water  Vapor  So'onder  (SSM/T-2  or  T-2  which  is  described  in 
more  detail  in  the  next  chapter)  heralds  a  new  era  in  water  vapor  studies  with  its  apparent 
capability  of  profiling  atmospheric  water  vapor  although  not  with  the  vertical  resolution 
of  current  radiosondes. 

Observing  system  simidation  experiments  (OSSEs)  wer^  carried  out  by  Hoffinan  et  ol. 
(1990)  to  attempt  to  determine  the  impact  of  the  Special  Sensor  Microwave/Temperatitre 
(SSM/T)  sensors,  in  addition  to  other  sensors.  In  these  experiments  the  addition  of  the 
SSM/T-2  and  SSM/T-l  data  improves  moissmre  analysis,  especially  in  the  tropics  and  the 
exlratiopics  of  the  Southern  Hemisphere.  Root  mean  square  (RMS)  errors  for  relative 
humidity  forecasts  were  decreased  from  29%  to  22%  in  the  OSSEs  with  the  SSM/T-2 
data  included.  In  the  Southern  Hemisphae,  the  SSM  data  Improved  the  500  hPa  height 
forecasts  by  12  hours. 

1.2  MICROWAVE/MILLIMETER  SATELLITE  RETRIEVALS 

The  potential  benefits  of  using  passive  microwave  sensors  to  measure  atmospheric, 
param.etm  was  first  recognized  in  the  1940’s.  Nonetheless,  it  was  only  during  the  early 
1960’s  wheii  the  first  microwave  radiometer  was  used  on  a  spacecraft  for  the  purpose 
of  atmospheric  remote  sensing.  The  U.S.  Mariner  2  Venus  Probe  m  1962  took  the  first 
satellite  microwave  emission  observations  of  a  planetary  atmosphere.  The  two-channel 
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radiometer  determined  an  upper  limit  on  water  vapor  in  Venw’  atm.osphere  and  retrieved 
a  planetary  surface  temperature  near  600°K  (Barath  et  al.,  1964).  The  first  U.S.  satellite 
to  make  microwave  observations  of  the  earth's  atmosphere  was  Nimbus  5  in  1972. 

Microwave  temperature  retrieval  on  satellites  operated  by  National  Oceanic  and  At¬ 
mospheric  Administration  (NOAA)  began  with  the  Nimbus  5  (Nimbus  E)  Satellite  Mi¬ 
crowave  Spectrometer  and  continued  with  the  Nimbus  6  Scanning  Microwave  Spectrometer 
(SCAMS).  The  SCAMS  used  three  channels  (52.85  GHz,  53.85  GHz,  and  55.45  GHz)  in 
the  oxygen  absorption  band  between  50  and  60  GHz.  The  follow-on  to  the  SCAMS  is  the 
Microwave  Sounding  Unit  (MSU).  It  serves  as  a  complement  to  the  infrared  temperature 
sounding  unit  called  the  High  Resolution  Infrared  Sounder  2  (HiRS/2).  The  MSU  has  four 
channels  (50.30  GHz,  53.74  GHz,  54.96  GHz,  and  57.95  GHz)  in  the  oxygen  absorption 
band.  The  MSU  is  a  ‘Dicke’  radiometer  where  the  output  of  the  receiver  is  proportional 
to  the  difference  between  the  brightness  temperature  of  the  scene  being  viewed  and  the 
temperature  of  the  internal  radiation  soiurce.  The  primary  purpose  of  the  MSU  is  to  make 
temperatme  sovmdings  in  the  presence  of  clouds  although  MSU  brightness  temperature  , 

data  are  also  currently  being  analyzed  to  understand  certain  aspects  of  the  Earth’s  climate 
(Spencer  et  oL,  1990). 

Water  vapor  retrieval  from  satellites  occurs  primarily  in  two  regions  of  the  electro- 

I 

magnetic  spectrum,  the  infrared  and  microwave/millimeter  wave  regions.  The  HlRS/2 
component  of  the  polar  orbiting  TIROS-N  operational  vertical  sounder  (TOYS)  (Smith  et 
al.,  1979;  Hillger  tuid  Yonder  Haar,  1981)  and  the  YISSR  atmospheric  sounder  (YAS)  on 

I 

Geostationary  Operational  Environmental  Satellites  (GOES)  (Smith,  1983)  are  examples 
of  infrared  instruments.  These  instruments  use  the  absorption  in  the  6.7  /mi  infrared  vi¬ 
bration  band.  They  have  had  limited  success  in  obtmning  operationally  useful  water  vapor 
profiles,  in  part  because  of  the  limitation  of  all  infrared  sensors  that  moisture  soimding  is  * 

impossible  at  levels  below  cloud  top. 

The  appe^  of  microwave/millimeter  wave  sounding  is  that  temperature  retrievals 
staffer  limited,  if'  aaiy,  degradation  due  to  the  presence  of  clouds.  For  years  the  concept  > 

of  tising  microwave  remote  sensing  for  v/ater  vjtpor  retrieval  was  discussed.  The  first 
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efforts  foc'ised  on  the  use  of  the  weak  HaO  rotational  line  at  22.235  GHz  for  sensing  total 
precipitable  water.  The  Nimbus  5  (Nimbus  E)  Satellite  Microwave  Spectrometer  (NEMS) 
iised  channels  at  22.235  and  31.4  GHz  to  measure  atmospheric  water  vapor  and  liquid 
water  over  the  oceans,  even  in  the  presence  of  many  types  of  clouds.  NEMS  also  measured 
the  atmospheric  temperature  profile  (as  mentioned  above)  from  0-20  km  (Staelin  et  aL, 

» 

1976). 

Millimeter  wavelengths  also  possess  the  advantage  of  low  emissivity  values  for  the 
ocean  smface  versus  infrared  wavelengths.  The  low  emissivity  values  provide  a  relatively 
cold  background  that  provides  enough  contrast  to  allow  the  retrieval  of  low  level  moisture  * 

fields  (Hoffinan  et  oL,  1990).  The  accuracy  of  these  low  level  retrieved  values  decreases 
with  high  surface  wind  speeds  due  to  the  formation  of  foam  on  the  ocean  siuface.  This 
foam  causes  the  emissivity  of  the  ocean  to  increase,  thus  effectively  reducing  the  contrast  I 

of  this  background  surface.  Additionally,  millimeter  wavelengths  are  able  to  sense  through 
clouds  so  we  can  measure  water  vapor  in  the  atmosphere  below  the  clouds. 

The  Special  Sensor  Microwave/Imager  (SSM/I)  which  has  been  carried  on  DMSP  ^ 

satellites  since  1987  is  also  used  to  retrieve  precipitable  water.  Jackson  (1992)  presents  a 
review  of  several  of  the  methods  used  for  retrieval  of  precipitable  water  with  the  SSM/I. 

The  SSM/I  also  ii.akes  use  of  the  water  vapor  line  at  22.235  GHz. 

► 

1.3  USE  OF  THE  183  GHz  ABSORPTION  LINE 

The  water  vapor  rotation  absorption  line  feature  at  183.31  GHz  is  about  20  times 
stronger  than  the  rotation  feature  at  22.235  GHz  (Waters,  1976).  For  this  reason  some  * 

researchers  have  (fiscussed,  in  both  abstract  and  definitive  terms,  its  potential  use  to  profile 
atmospheric  water  vapor  since  the  1970’s. 

Gaut  et  oL  (1975)  performed  a  study  of  a  variety  of  microwave  remote  sensing  sys-  I 

terns  to  satisfy  Air  Force  meteorological  data  requirements.  The  simulation  and  retrieval 
exercises  investigated  the  following  parameters:  (a)  total  integrated  water  vapor,  (b)  the 
vertical  distribution  of  water  vapor,  (c)  temperature  profile,  and  (d)  the  integrated  and  ^ 

vertical  distribution  of  cloud  liquid  water.  While  emphasizing  channels  below  60  GHz, 
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higher  frequency  channels  (including  those  in  the  vicinity  of  183  GHz)  were  also  evalu¬ 
ated.  The  higher  frequency  channels  were  explored  because:  (a)  it  was  desired  to  remove 
the  effects  of  variable  surface  emissivity  over  land  at  the  lower  frequencies  used  to  obtain 
water  V£q>or  information,  (b)  it  was  hypothesized  that  the  higher  frequency  channels  could 
provide  information  on  cloud  vertical  structtire,  and  (c)  higher  frequency  channels  were 

» 

expected  to  be  more  sensitive  to  the  integrated  water  vapor  of  clouds  with  relatively  small 
liquid  water  content  (Isaacs,  1987). 

The  pioneering  work  in  the  theory  of  using  the  183  GHz  rotation  line  for  atmospheric 
remote  sensing  was  written  by  Schaerer  and  Wilheit  (1979).  Several  authors  (Kakar,  1983;  * 

Kakar  and  Lambrightsen,  1984;  and  Rosenkranz  et  aL,  1982)  have  since  discussed  various 
methods  for  the  retrieval  of  clear  sky  atmospheric  moisture  profiles  using  channels  around 
the  183  GHz  line.  > 

Subsequently,  the  theory  has  been  extended  to  include  clouds  in  the  instrument  field  of 
view  (FOV).  Isaacs  and  Deblonde  (1987)  examined  some  of  the  effects  of  beam-filling  liquid 
water  clouds  on  millimeter  wave  moisture  retrievals.  They  compared  retrieved  soundings  ^ 

with  radiosondes  from  clear  sky  and  cloudy  conditions  in  simulation  studies.  WUheit 
(1990)  performed  similar  simulation  studies  expanding  on  the  work  he  did  previously  with 
Schaerer.  Both  of  these  simulation  studies  were  limited  by  the  constrmnt  of  only  one  cloud 
layer. 

Wang  et  aL  (1983)  showed  experimentally  that  the  183  GHz  line  could  be  used  for 
profiling  atmospheric  water  vapor.  They  used  an  instnnnent,  the  Advanced  Microwave 
Moisture  Sounder  (AMMS),  from  an  airborne  platform  and  this  work  was  for  clear  sky 
cases.  Lutz  et  al.  (1991)  used  the  algorithm  developed  by  Wilheit  for  profiling  atmospheric 
water  vapor  even  in  the  presence  of  clouds.  This  work  again  used  measurements  from  the 
AMMS.  They  showed  the  retrieved  profiles  were  in  general  agreement  with  those  from  ^ 

radiosonde  data.  Lutz  et  al.  also  noted  that  the  algorithm  did  not  perform  very  well  in 
the  vidnity  of  surface  fronts.  With  the  concept  proven  by  these  aircraft  studies  the  next 
step  was  to  put  a  moisture  soimder  on  a  satellite.  I 


» 


1.4  MOTIVATION,  OBJECTIVES,  AND  THESIS  OUTLINE 


1.4.1  Motivation 

The  attraction  for  working  with  SSM/T'2  data  is  the  opportunity  to  work  data  from 
the  newest  microwave  instrument,  also  the  first  with  frequencies  near  183  GHz,  on  an  Air 
Force  meteorological  satellite.  The  work  on  data  from  this  instrument  has  barely  begun 
and  the  ability  to  learn  about  some  of  the  capabilities  of  the  SSM/T-2  is  a  special  oppor- 
timity.  The  ability  to  retrieve  moisture  soundings  over  the  entire  globe  regul^u'ly  is  also 
appealing  to  the  needs  of  GEWEX  and  other  major  global  climate  research  programs. 
The  T-2  will  also  help  improve  the  output  of  numerical  weather  prediction  models.  Addi- 
tionall}',  during  the  testing  of  the  capabilities  of  an  instrument  a  researcher  may  also  learn 
of  its  limitations,  whether  they  be  inherent  or  due  to  the  advancement  of  the  science. 

1.4.2  Objectives 

This  thesis  attempts  to  achieve  the  following  objectives: 

1.  Develop  an  imderstan<hng  of  nsany  of  the  sensitivity  factors  which  affect  the  mi¬ 
crowave  brightness  temperatures,  specifically  for  the  SSM/T-2  frequencies. 

2.  Relate  T-2  brightness  temperature  features  to  characteristics  of  the  general  circula¬ 
tion  of  the  atmosphere  and  characteristics  of  the  earth’s  surface. 

3.  Develop  a  simple  physical  retrieval  of  columnar  water  vapor  which  may  be  used  to 
assist  in  the  interpretation  of  the  T-2  data  and  which  may  be  used  in  initializa¬ 
tions  of  more  complex  schemes.  This  retrieval  scheme  assumes  the  thermodynamic 
temperature  and  brightness  temperature  are  equal  at  a  level  of  approximately  equal 
water  vapor  overburden  (column  water  vapor  above  this  level)  irrespective  of  the 
temperature  and  moisture  profiles. 

1.4.3  Thesis  Outline 

Chapter  2  overviews  the  characteristics  of  the  satellites  which  carry  the  SSM/T-2  and 
the  characteristics  of  the  instruments  flown  on  these  satellites.  Chapter  2  also  presents 
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brightness  temperature  maps  for  two  T-2  data  periods  and  an  analysis  of  these  data  is  pre- 
sentedv  Chapter  2  concludes  with  scatter  plots  of  brightness  temperatures  from  two  regions 
which  demonstrate  distinctly  different  temperature  and  moisture  properties.  Chapter  3 
develops  the  radiative  transfer  equation  as  it  applies  to  microwave  remote  sensing.  It  con¬ 
tinues  with  a  discussion  of  retrievals,  first  temperatme  (briefly)  and  then  water  vapor  and 
introduces  the  concepts  of  weighting  functions.  Chapter  3  concludes  with  a  description  of 
various  models  (atmospheric  propagation,  atmospheric  temperature/pressure,  and  cloud) 
used  in  the  simiilations  discussed  in  this  work.  Chapter  4  analyzes  the  response  of  the 
DMSP  SSM/T-2  channels  to  a  variety  of  factors  (what  is  called  the  forward  problem) 
including:  the  profile  of  atmospheric  moisture,  the  presence  of  clouds,  the  characteristics 
of  clouds,  and  the  backgroimd  against  which  the  atmosphere  is  viewed.  The  analysis  of 
cloud  effects  here  is  limited  only  to  low  level  water  clouds.  The  possible  effects  of  ice 
scattering  at  183  GElz  is  not  considered  in  this  study  althou^  it  may  be  important  imder 
some  circumstances.  Chapter  5  introduces  the  concept  of  water  vapor  burden  and  through 
the  use  of  figures  and  tables,  it  is  shown  that  the  vapor  burden  is  (nearly)  constant  a  a 
^ven  frequency  near  the  183  GHz  absorption  line  and  is  only  slightly  dependent  on  the 
moisture  and  temperatme  profile.  The  water  vapor  burden  is  mapped  onto  a  constant 
pressure  surface  thus  providing  the  water  vapor  burden  above  this  surface.  Chapter  6 
presents  conclusions  based  on  the  work  contained  herein  and  a  discussion  about  possible 
applications  of  the  SSM/T-2  instrmnent  and  the  possibilities  of  unified  retrievals  tising 
other  instruments  from  the  DMSP  microwave  suite. 
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Chapter  2 

THE  DMSP  SSM/T-2  AND  T-1  INSTRUMENTS  AND  DATA  ANALYSIS  ' 

2.1  DMSP  Block  5D  Satellite  Characteristics 

DMSP  Block  5D  Satellites  operate  in  a  sun-synchronous,  near-polar  orbit  with  a  ) 

nominal  altitude  of  800  km.  The  inclination  an^e  of  the  orbit  is  98.8°  and  an  orbit  period 
of  102.0  min.  The  orbit  produces  14.1  orbit  revolutions  of  the  satellite  per  day.  The 
combination  of  these  parameters  produces  data  void  areas,  diamond-shaped  in  appearance  ^ 

near  the  equator,  in  a  24  hour  period.  This  missing  coverage  has  obvious  drawbacks  when 
the  data  are  to  be  applied  to  study  the  trcq>ics  (such  as  in  the  use  of  SSM/I  to  monitor 
tropical  cyclone  intensity  and  position).  These  data  void  regions  shift  with  each  orbit  and 
complete  coverage  of  the  tropical  regions  is  achieved  after  72  hours.  The  orbit  inclination 
also  results  in  circular  sectors  of  280  km  at  both  poles  that  are  never  sampled.  The 
swath  vddth  for  the  SSM/I  and  SSM/T  channels  is  1400  km  perpendicular  to  the  satellite 
subtrack.  The  SSM/I  is  a  conically  scanning  radiometer  (with  a  constant  scan  angle  of  45 
degrees  aft  of  the  satellite)  as  compared  to  the  SSM/T  radiometers  which  are  cross-track 
scanning  sounders  (which  scan  in  a  fashion  similar  to  the  OLS  described  below). 

The  primary  instrument  on  the  DMSP  satellites  is  the  Operational  Linescan  System  ^ 

(OLS)  which  is  a  fom:  channel  imager.  The  OLS  visible  channel  for  very  high  resolution  has 
a  resolution  of  600  m.  Since  the  bandwidth  of  this  channel  is  very  broad,  the  instrument 
receives  more  radiation  from  a  given  scene  (e.g.  as  compared  to  the  imagers  on  NOAA  I 

polar  orbiter  satellites)  and  hence  has  a  greater  sensitivity.  The  maximtnn  swath  width 
for  the  OLS  is  2500  km.  The  OLS  produces  the  visible  and  IR  imagery  used  for  many 
operational  requirements  by  forecasters  in  the  Department  of  Defense  (DOD). 
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The  F8  DMSP  spscecrail,  launched  in  1987,  was  the  firs;  satellite  with  the  Special 
Sensor  Microwave/Imagar  (SSM/I).  The  SSM/I  is  a  seven-channd  passive  microwave  ra¬ 
diometer  which  operates  at  four  frequencies.  It  receives  both  vertically  and  horizontally 
linearized  radiation  at  19.3,  37.0,  and  85.5  GHz  and  vertical  only  at  22.2  GHz  (HoUinger 
et  cd.,  1990).  Data  from  tlus  instrument  are  used  for  the  determination  of  the  following 
environmental  parameters: 

1.  Ocean  surface  wind  speeds. 

2.  Ice  coverage,  age  and  extent. 

3.  Cloud  water  content. 

4.  Integrated  water  vapor  (precipitable  water). 

5.  Precipitation  over  water. 

6.  Soil  moistiupe. 

7.  Land  surface  temperature. 

8.  Snow  water  content. 

9.  Cloud  amoimt. 

The  SSM/I  is  also  flown  on  Fll  together  with  the  SSM/T-2  and  SSM/T-1.  These 
instruments  being  on  the  same  satellite  offer  the  possibility  of  a  synergistic  (also  called  a 
unified)  approach  to  atmospheric  retrieval  where  deficiencies  of  one  instrument  perhaps 
can  be  diminished  or  removed  through  the  use  of  another  instrument. 

2.1.1  SSM/T-2  Channel  Characteristics 

The  Defense  Meteorological  Satellite  Program  (DMSP)  Special  Sensor  Microwave 
(also  called  Millimeter  Wave)  Watt?  Vapor  Soimder  (SSM/T-2)  currently  flies  on  the  F- 
11  satellite  which  was  launched  on  November  28, 1991.  The  SSM/T-2  possesses  5  chaimels: 
three  dual-pass  bands  located  on  the  183  GHz  (1.64  mm)  water  vapor  absorption  line,  one 
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on  the  line’s  wing  (150  GHz),  and  a  window  channel  (91.655  GHz).  The  local  descending 
nodal  crossing  time  of  F-11  is  approximately  0514.  The  T-2  is  a  tottd-power  radiometer 
which  measures  the  total  power  seen  by  the  antenna  over  the  signal  processing  integration 
time.  The  objective  of  the  SSM/T-2  is  to  provide  global  measurements  of  the  vertical 
profile  of  water  vapor  to  support  Air  Force  applications  including  the  input  of  these  profiles 
to  munerical  weather  forecast  models.  The  characteristics  of  the  SSM/T-2  channels  are 
shown  in  Table  2.1  (including  the  noise  equivalent  temperature  imcertainty,  NEDT). 


Table  2.1:  SSM/T-2  Channel  Characteristics  (after  Griffin  et  aL,  1993). 


Chan. 

No. 

Center 

Frequency 

(GHz) 

Nadir 

FOV 

(km) 

Beam 

Acceptance 

(degrees) 

Peak 

Altitude 

(hPa) 

Pol. 

NEDT 

(°K) 

Response 

183.31  ±  3 

48 

3.3 

650 

■9 

water  vapor 

183.31  ±  1 

48 

3.3 

500 

mm 

water  vapor 

183.31  ±  7 

48 

3.3 

800 

■i 

water  vapor 

91.655 

88 

6.0 

Smface 

0.6 

surface 

5 

150.0 

54 

3.7 

_ 1 

1000 

■I 

surface 

The  T-2  wa-s  originally  designed  to  be  a  four  channel  instrument  (minus  the  channel 
near  92  GHz).  The  work  of  Isaacs  and  Deblonde  (1985)  showed  that  retrievals  were 
not  very  accurate  imder  some  conditions  (specifically,  tropical  atmospheres)  without  the  I 

91.655  GHz  channel. 

The  T-2  views  both  an  internal  hot-load  target  (~  SO'^'^K)  and  cosmic  background 
radiation  for  calibration  measurements.  The  T-2  also  underwent  a  calibration  study  led  by 

► 

the  Geophysics  Directorate  of  the  Air  Force’s  Phillips  Laboratory.  This  calibration  study 

included  independent  measurements  and  model  calculation  studies.  These  independent 

measurements  were  performed  by  the  NASA  Millimeter-wave  Imaging  Radiometer  (MIR) 

carried  on  the  NASA  ER-2  sdrcraft  (Racette  et  al.,  1992).  The  MIR  and  T-2  contsun 

essentially  the  same  chaimels  with  the  exception  of  the  window  channel  which  for  the  MIR 

is  located  at  89  GHz.  The  NASA  ER-2  equipped  with  the  MIR  imdeflew  selected  DSMP 

satellite  passes  including  ocean,  coastal,  and  land  background  cases.  RMS  differences  of  ^ 

0.7- 1.4 °K  were  measured  between  the  T-2  and  the  MIR. 
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Radiative  tra!7.sfer  model  calculations  were  also  conducted  (at  coincident  times  and 
locations  with  the  model  of  Eyre  and  Woolf  (1988))  using  radiosonde  data  and  estimates 

i 

of  the  surface  emlttances.  Although  the  comparisons  of  the  model  calculations  with  the 
T-2  data  showed  collocated  RMS  difr(n'ences  of  6.5°K,  the  calibration  team  concluded 
based  on  the  data  they  luid  obtained  that  the  SSM/T-2  suffers  no  significant  bias  in  its 
calibration  (Griffin  et  al.,  1993). 

2.1.2  T-2  Output  Parameters 

The  following  parameters  are  derived  from  T-2  data  using  the  operational  retrieval  ^ 

retrieval  scheme  of  A.  Stogryn  (Boucher  ef  aL,  1993): 

1.  Relative  humidity  at  1000.  ??50,  700,  500,  400,  and  300  mb  levels. 

2.  Specific  humidity  at  1000,  830,  700,  500,  400  and  300  mb  levels. 

3.  Water  vapor  mass  (kg/m^)  between  levels  surface-1000, 1000-850,  850-700,  700-500, 

500-300,  and  above  300  mb.  ^ 

Although  these  parameters  ai-e  produced  operationally,  the  actual  retrieval  of  data 
on  water  vapor  from  T-2  data  is  complex  for  reasons  described  in  the  next  chapter  £md 
there  remains  considerable  scope  for  research  in  understanding  these  complexities.  ( 

2.1.3  SSM/T-1  Chaimel  Characteriaiics 

The  Special  Sensor  Microwave/Temperature  Soimder  (SSM/T-1)  was  developed  and 
built  by  Aerojet  Electrosystems  Company  for  the  T.S-  Air  Force  Space  Systems  Division 
Defense  Meteorolo^cal  Satellite  Program  (DMSP).  The  design  methodology  was  a  top- 
down  systems  development  approach.  The  SSM/T-1,  first  launched  in  1978,  was  designed 
to  support  Air  Force  and  Department  of  Defense  operational  weather  forecasting,  specif-  ^ 

ically  the  data  required  by  the  Air  Force  Global  Weather  Central  (AFGWC)  numerical 
weather  prediction  (NWP)  models.  To  incorporate  the  T-1  data,  the  T-1  is  designed  to 
emulate  radiosonde  temperature  and  pressure  measurements,  i.e.  retrieve  temperatu’*es  at  I 

the  standard  pressiu^  levels,  thicknesses  between  these  levels,  and  the  temperature  and 


pressure  of  the  tropopause.  Channel  selection  was  such  that  the  weighting  functions  are 
distributed  as  evenly  as  possible  throughout  the  atmosphere  (Falcone  and  Isaacs,  1987). 
Weighting  functions  for  channels  1-4  of  the  T-1  were  shown  in  Figure  3.1.  The  channel 
characteristics  for  the  SSM/T-1  are  shown  iu  Table  2.2. 


Table  2.2:  SSM/T-1  Channel  Characteristics. 


Center 

R-equency 

(GHz) 

Nadir 

FOV 

(km) 

Beam 

Acceptance 

(degrees) 

Peak 

Altitude 

(km) 

NEDT 

(°K) 

Response 

50.5 

204 

14 

surface 

53.2 

204 

14 

2 

T  at  2  km 

54.35 

204 

14 

6 

0.4 

T  at  6  km 

54.9 

204 

14 

10 

0.4 

T  at  10  km 

58.4 

204 

14 

30 

0.5 

T  at  30  km 

58.825 

204 

14 

16 

0.4 

T  at  16  km 

59.4 

204 

14 

22 

0.4 

T  at  22  km 

The  SSM/T-1  is  a  passive  ‘Dicke’  radiometer  which  scans  in  a  cross-track  fashion.  The 
radiometer  has  a  scan  time  of  32  seconds  and  a  dwell  time  for  each  scene  of  2.7  seconds;  the 
instrument  takes  0.3  seconds  between  each  scene  and  the  warm  ^md  cold  calibrations  occm* 
during  the  remaining  time.  The  SSM/T-1  operates  in  an  oxygen  absorption  band  using 
seven  frequencies  between  50  and  60  GHz  and  is  similar  to  the  MSU  which  was  described 
above.  The  variation  of  the  dry  sdr  absorption  at  each  of  these  frequencies  provides  the 
basis  for  the  retrieval  of  temperature  profiles  (Curtis  and  Shipley,  1987). 

2.1.4  SSM/T-2  and  T-1  Scan  Patterns 

The  scan  pattern  for  the  SSM/T-2  is  similar  to  the  SSM/T-1.  Since  the  SSM/T-2 
operates  at  a  longer  wavelength  than  the  SSM/T-1  and  hence  has  a  smaller  footprint, 
or  field  of  view,  the  SSM/T-2  samples  four  times  as  many  scene  stations  (28  versus  7) 
per  scan  revolution,  and  must  scan  four  times  as  fast  to  provide  coincident  samples.  The 
resiilt  is  that  the  SSM/T-2  produces  16  times  as  many  total  stations  scenes.  The  footprint 
pattern  is  shown  in  Figure  2.1.  The  reason  for  coincident  samples  is  that  to  produce  water 
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vapor  retrievals  the  T-2  operational  algorithm  also  requires  four  of  the  seven  measured  T- 
1  frequency  channels  (Channels  1-4)  to  provide  the  requisite  temperature  soxmding.  The 

I 

pattern  in  Figure  2.1  is  changed  by  the  presence  of  a  glare  obstructor  or  GLOB,  on  the 
satellite.  The  GLOB  is  designed  to  protect  the  OLS  from  solar  insolation  which  coiild 
degrade  the  instrument  but  this  also  means  that  1  scene  of  T-1  data  and  4  scenes  of  T-2 
data  on  the  sun  side  of  the  satellite  are  degraded  and  therefore  not  used.  ' 

2.2  SSM/T-2  DATA 

The  SSM/T-2  data  used  for  this  thesis  was  kindly  provided  by  Mr.  Bruce  Thomas  ^ 

of  the  DMSP  Omaha  Field  Office  of  Aerospace  Corporation.  The  data  were  obt^dned  for 
approximately  78  revolutions  of  satellite  F-11  in  March  1993.  The  data  package  received 
from  Aerospace  Corporation  consisted  of  T-2,  T-1,  and  SSM/I  data  in  a  VMS  back-up  save 

» 

set  format.  All  data  received  from  Aerospace  were  in  sensor  data  record  (SDR)  format 
(also  known  as  antenna  pattern  corrected  temperatures). 

The  T-2  data  set  consists  of  13  words  in  each  record.  Word  1  contzuns  information  on 
the  scan  scene  (which  provides  the  scan  angle  for  the  sensor)  and  the  background  against  ^ 

which  the  sensor  is  viewing  (land,  ocean,  sea  ice,  or  coast).  The  reason  for  the  differen¬ 
tiation  between  backgrounds  is  to  have  a  mechanism  to  help  facilitate  the  determination 
of  the  surface  emission  for  each  viewing  scene.  A  discussion  of  emissions  for  varioiis  back-  ^ 

grounds  is  included  in  Chapter  4.  Word  2  cont^s  the  Julian  hour  (referenced  to  00  UTC 
31  December  1967),  and  words  3  and  4  are  the  Julian  minute  €uid  second,  respectively. 

Word  5  is  the  longitude  (from  0.-360°  East)  and  word  6  is  the  latitude  (±  90°)  which  ^ 

together  determine  the  earth  location  of  the  scan  scene.  Word  7  provides  the  terrain 
height  (m)  and  word  8  is  the  lOOOhPa  height  (m)  from  the  operational  forecasting  model 
at  AFGWC.  Words  9-13  are  the  sensor  data  records  (SDRs)  for  the  five  T-2  channels 
(Kelvin)  (B.  Thomas,  personal  commxmication). 

2.3  T-2  BRIGHTNESS  TEMPERATURE  MAPS 

The  T-2  brightness  temperatures  fields  for  the  available  data  from  March  1993  were  ► 

mapped  and  contoured.  Since  the  4  scan  scenes  on  one  side  of  the  swath  are  unusable 
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because  of  the  shieiding  by  the  GLOB  as  meatioued  above  aad  because  of  the  brightness 
temperature  changes  as  a  fiinction  of  scan  scene  angle  (discussed  in  fsjrther  detail  in 
Chapter  4),  the  outer  4  scenes  on  the  othe?  side  of  the  scan  were  also  removed.  This 
results  in  a  scan  line  that  is  considerably  less  than  the  1400  km  swath  width  mejitloned 
above  but  it  does  mean  a  data  set  which  is  more  representative  requiring  no  corrections 
to  the  T-2  brightness  temperatures  for  angle  effects.  The  result  of  this  is  that  there  are 
areas  during  the  two  periods  which  are  not  sampled  and  these  areas  aie  represented  by 
gray  in  the  colored  diagrams. 

Data  were  binned  into  1  degree  latitude  x  longitude  bins  and  then  averaged.  Since 
resolution  at  nadir  is  about  50  km,  bins  twice  the  size  of  nadir  resolution  should  be 
representative  for  the  data. 

2.3,1  Channel  1  Brightness  Temperatui-e  Maps 

Channel  1,  at  183.31  ±  3  GHz,  is  used  for  the  determination  of  mid-upper  tropospheric 
water  vapor  and  relative  humidity  information.  As  shovm  below  (bi  Figures  4.i-4.3), 
the  peak  of  the  weighting  function  for  this  channel  lies  between  3-6  km  for  the  three 
McCiatchey  atmospheres.  Figures  2.2  and  2.3  are  the  Channel  1  brightness  temperatures 
for  the  periods  10-14  March  and  16-18  March  1993,  respectively. 

The  areas  in  red  which  are  located  both  north  and  south  of  the  equator  (e.g.  roughly 
20-25°  from  the  equator)  delineate  the  location  of  subsidence  associated  with  the  sub¬ 
tropical  highs.  This  prominent  feature  is  also  prevalent  in  the  distributions  for  other 
frequencies.  The  subtropical  high  in  the  North  Atlantic  extends  from  western  Africa 
westward  all  the  way  to  the  Caribbean  Sea. 

Figure  2.3  suggests  a  more  coherent  and  stronger  ITCZ  in  the  Pacific  Ocean  (the 
white  and  green  color  shading  with  several  purple  pixels).  The  South  Pacific  Convergence 
Zone  (SPCZ)  (from  the  central  South  Pacific  extending  northwest  to  east  of  New  Guinea) 
is  also  more  distinct  in  this  figure.  Both  features  appear  cold  (and  thus  moist)  in  contrast 
to  the  the  warm  (and  thus  dry)  subsiding  air  poleward  of  these  features. 
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2.3.2  Channel  2  Brightness  Temperature  Maps 

The  weighting  functions  of  Channel  2  peak  highest  in  the  atmosphere  (Figux*es  4. 1-4.3  > 

with  peaks  between  5-8  km)  and  measurements  provided  by  Chaimel  2  thtis  applies  to 
upper  tropospheric  moisture.  Figure  2.4  aad  2.5  are  the  Channel  2  brightness  temperature 
maps  for  the  two  time  periods  described  above.  Maps  for  this  channel  have  less  spread  ^ 

of  brightness  temperatures  than  for  Channel  1.  Channel  2  shows  several  areas  with  high 
brightness  temperatures  over  the  oceans  and  these  are:  (1)  in  the  central  North  Pacific 
near  10-15’  latitude  and  near  the  date  line,  (2)  in  the  eastern  North  Pacific  south  of  Baja, 

California,  (3)  in  the  eastern  North  Atlantic  off  the  west  coast  of  Afidca,  (4)  in  the  South 

Atlantic  off  the  coast  of  Brazil  and  (5)  near  Madagascar  in  the  Indian  Ocean.  These 

brightness  temperatures  indicate  that  the  atmosphere  is  very  dry  in  the  upper  levels  and 

these  regions  are  also  over  the  subsidence  re^ons  indicated  in  Figures  2.2  and  2.3.  The  ^ 

coldest  brightness  temperatures  over  the  tropical  oceans  occur  near  Fiji  in  the  western 

Pacific  and  in  the  central  Indian  Ocean.  The  cold  brightness  temperatures  indicate  deep 

convection  in  these  areas  and  elevated  amounts  of  upper  tropospheric  moisture.  I 

2.3.3  Channol  3  Brightness  Temperature  Maps 

CK'Ujnei  3,  lo-^ated  at  183.31  d:  7  GRz,  senses  levels  of  water  vapor  that  occur  deeper 
into  the  atmosphpre  accovrimg  to  tho  location  of  the  weighting  function  maxima  for  this  * 

channel  that  are  depicted  in  Chapter  4.  Figun'  2.6  and  2.7  are  the  brightness  temperatme 
maps  for  lt}-14  Mard’  and  16-18  March,  respectively.  Once  again,  there  exists  coherence  in 
the  subsidoice  regions  such  that  ths  broad  areas  of  red  color  tU-e  the  regions  of  subtropiaU  I 

high  pres.suies  imderlying  those  dry  regions  noted  in  reference  to  the  maps  of  Channels 
1  and  2  brightness  temperature.  Similarly,  the  moist  tropical  regions  where  the  deepest 
convection  occurs  appear  relatively  cold  (indicated  by  white-green  colors).  ^ 

2.3.4  Channel  4  Brightness  Temperature  Maps 

Channel  4,  the  window  charinel  at  91.655  GHz,  primarily  senses  the  surface  of  the 
earth.  Figures  2.8  and  2.9  are  show  the  brightness  temperatures  from  Channel  4  for  the  > 

time  periods  described  above.  The  red  color  across  the  land  masses  (especially  Australia) 


» 


in  the  southern  hemisphere  poict  to  clear  skies  (at  least  predominantly),  a  warm  siudace 
and  a  high  emissivity  for  these  land  masses.  This  is  not  surprising  when  considering  that 
the  land  has  had  the  entire  siiramer  to  become  warm  and  dry.  The  lower  emissivity  of 
the  adjacent  ocean  surfaces  resiilts  in  much  lower  brightness  temperatures  and  it  is  easy 
to  diiferentiate  between  land  masses  and  ocean  surfaces. 

East  of  Japan  there  are  four  curved  features  with  brightness  temperatures  warmer 
than  the  surroimding  region.  D^y  maps  of  brightness  temperattires  (not  shown)  indicate 
this  feature  is  a  mesoscale  storm  system  which  has  been  viewed  by  the  T-2  on  successive 
satellite  passes  (over  several  days).  The  cloud  features  associated  writh  this  low-pressure 
system  show  up  as  warm  brightness  temperatures  over  the  cold  ocean  backgroimd,  resem¬ 
bling  a  comma  cloud. 

Regions  of  the  tropical  oceans  show  little  structure.  The  yellow  and  brown  colors 
extend  farther  south  in  the  Southern  Hemisphere  as  compared  to  the  northern  extent  in 
the  Northern  Hemisphere  due  to  warmer  sea  smface  temperatmes  and  low  level  moisture 
over  these  oceans,  since  the  weighting  functions  corresponding  to  this  channel  peak  at  /near 
the  surface.  Since  SSTs  typically  lag  by  about  two  months,  SSTs  in  the  southern  oceans 
should  be  just  beginning  their  fall  from  the  yearly  temperature  maYiTnuTn. 

2.3.5  Channel  5  Brightness  Temperature  Maps 

Channel  5,  at  150  GHz,  is  located  on  the  line  wing  of  the  absorption  line  at.  183.31 
GHz.  As  shown  in  Chapter  4,  the  weighting  function  for  this  channel  pe£iks  in  the  lower 
troposphere  or  at  the  earth’s  surface  (depending  on  the  moistme  profile  of  the  atmosphere). 
Figures  2.10  and  2.11  are  the  brightness  temperature  maps  for  the  periods  10-14  and  16-18 
March  1993.  The  maps  hint  at  the  ITCZ  and  associated  convection  (with  the  light  yellow 
and  white  pixels)  in  the  tropical  oceans.  Otherwise,  there  is  not  a  great  deal  of  structure 
to  the  temperature  differences  that  hasn’t  already  been  said  in  the  previous  subsections. 

Figure  2.11  does  indicate  a  stronger  SPCZ  than  does  Figure  2.10.  Australia  has  high 
brightness  temperatures  over  the  entire  coimtry  likely  indicating  a  high  emissivity  for  the 
land  and  a  dry  atmosphere  above  it.  Agmn,  the  land  surfaces  are  easily  found  in  the 
brightness  temperature  maps  due  to  the  emissivity  contrasts  with  ocean  stufaces. 


The  brightness  temperature  maps  in  Figure  2.2-2.11  provide  {m  opportunity  to  ex¬ 
amine  some  of  the  general  circulation  features  of  the  atmosphere  and  to  see  the  surface 
characteristics  of  ocean  and  land  masses  at  the  frequencies  of  the  SSM/T-2. 

2.4  CHANNEL  COHERENCE 

A  helpful  way  of  analyzing  and  understanding  the  properties  of  the  T-2  channel  bright¬ 
ness  temperatures  is  to  consider  the  relationships  between  channels.  Here  we  consider  two 
regions  which  are  characterized  by  distinctly  different  large  scale  circulation  features  and 
thus  different  temperattire  and  moisture  properties. 

The  two  areas  selected  are  located  over  the  Pacific  Ocean.  The  first  area  corresponds 
approximately  to  the  TOGA  CO  ARE  (Coupled  Ocean- Atmosphere  Response  Experiment) 
donudn  and  is  defined  as  the  area  160-170°  East  longitude  and  0-10°  South  latitude.  This 
area  is  largely  under  the  control  of  deep  convection  and  ascent  as  part  of  the  Hadley  and 
Walker  circulations  and  represents  a  re^on  of  relatively  high  moistme  content  at  all  levels. 
The  brightness  temperatures  from  this  area  are  compared  and  contrasted  with  an  area  in 
the  Eastern  Pacific  defined  by  the  coordinates  110-120°  West  longitude  and  10-20°  North 
latitude.  This  East  Pacific  area  is  located  in  a  region  dominated  by  a  subtropical  high 
with  its  (assumed)  assodated  strong  subsidence  and  lower  colmnn  water  vapor.  Time 
constrmnts  did  not  permit  a  detailed  analysis  of  the  doudiness  over  these  regions  during 
this  study. 

The  brightness  temperatures  for  these  regions  and  for  the  5  channels  are  plotted 
against  each  other  in  the  form  of  scatter  dif^ams.  Figure  2.12  shows  Channel  5  versus 
Channel  4  for  the  two  Pacific  regions.  The  East  Padfic  (EP)  region  shows  a  very  small 
scatter  for  the  brightness  temperature  data.  The  Channel  4  temperatures  are  relatively 
low  (as  compared  to  the  Channel  5)  since  this  channel  is  designed  to  see  the  surface. 
As  discussed  in  Chapter  4,  the  emissivity  of  the  ocean  surface  is  relatively  low  (~  0.7) 
and  it  provides  a  cold  backgroimd.  Channel  5,  on  the  line  wing,  provides  both  surface 
effects  and  the  effects  of  low  level  water  vapor.  The  brightness  temperature  indicates  an 
emission  temperature  near  280  K.  The  West  Padfic  (WP)  area  shows  a  window  brightness 
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Figure  2.2:  Channel  1  brightness  temperatiu-e  map  for  10-14  March  1993. 
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Figure  2.3:  Channel  1  brightness  temperature  map  for  16-18  March  1993. 
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Figure  2.4:  Channel  2  brightness  temperature  map  for  10-14  March  1993. 
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Figure  2.5:  Chauinel  2  brightness  temperature  map  for  16-18  Meuch  1993. 


25 


I 


1  uegnee  Din  Hveroges 
March  10-14 »  1993 

Figure  2.6:  Channel  3  brightness  temperature  map  for  10-14  March  1993. 
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temperature  of  10-15  degrees  higher  than  the  East  Pacific  area.  This  likely  occurs  due 
to  more  ^nission  by  the  water  vapor  (near  the  surface,  enutting  at  a  high  temperature) 
in  the  wetter  atmosphere  in  the  West  Pacific.  The  scatter  of  data  from  the  ’West  Pacific 
indicates  that  the  area  contains  many  inhomogeneities.  Possible  reasons  for  this  scatter 
will  be  discussed  later. 

Figure  2.13  depicts  the  brightness  temperatures  of  the  150  GHz  chaimel  and  the  183 
±  7  GEfr  channel.  Since  these  two  channels  are  ac^acent  to  each  other,  it  is  expected  that 
the  two  channels  would  be  well  correlated.  For  the  East  Pacific  area,  the  Channel  3  and 
5  temperatures  are  nearly  equal,  with  most  of  the  Channel  3  temperatures  within  10°  of 
280  K.  For  the  WP  region,  Channel  3  becomes  opaque  more  quickly  and  therefore,  most 
of  the  temperatures  are  below  275  K.  The  vertical  correlation  of  brightness  temperature 
(and  therefore  water)  shows  up  very  well  for  the  WP  repon.  An  almost  linear  relationship 
between  these  two  channels  is  inferred  by  the  data. 

Figure  2.14  is  the  scatter  plot  of  the  brightness  temperatures  of  Channels  1  and  5. 
The  EP  area  again  exhibits  higher  brightness  temperatures  for  channels  which  peak  above 
the  surface,  implying  a  dry  upper  atmosphere.  Additionally,  there  is  not  much  scatter  of 
the  data  for  tlus  region,  implying  a  profile  which  is  nearly  invariant.  Channel  1  brightness 
temperatures  for  the  WP  region  are  about  20°  colder  than  the  EP  region.  This  indicates 
that  this  chsmnel  is  sensing  water  vapor  at  a  colder  temperature  which  means  a  higher 
altitude. 

Perhaps  the  best  depiction  of  the  differences  in  the  two  areas  is  shown  in  Figure 
2.15.  This  figiire  plots  the  window  channel  (Channel  4)  against  the  channel  closest  to  the 
absorption  line  (Chaimel  2).  Figure  2.15  shows  that  the  EP  area  has  a  high  Channel  2 
brightness  temperature,  implying  the  channel  peaks  well  down  in  the  atmosphere.  The 
WP  area  has  noticeably  lower  brightness  temperatures  for  Channel  2,  indicating  a  peak  of 
the  channel  higher  in  the  vertical  (and  consequently,  at  a  lower  thermometric  temperature 
than  for  the  EP).  Channel  4  brightness  temperatures  are  governed  by  the  low  emissivity 
of  the  cold  ocean  surface  in  the  EP  and  the  higher  thermal  temperature  of  the  surface  and 
lower  atmosphere  combined  for  the  WP  region. 
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Figure  2.12:  Channel  5  versus  Cheinnel  4  brightness  temperatures  for  the  two  Pacific  areas. 
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Figure  2.13:  Channel  5  versus  Channel  3  brightness  temperatures  for  the  two  Pacific  areas 
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Figure  2.14:  Same  as  Figure  2.12  except  Channel  5  vs.  Channel  1. 
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Figure  2.15:  Channel  4  versus  Channel  2  brightness  temperatures  for  the  two  Pacific  areas 
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Figure  2.16  also  shows  a  distinct  difference  in  the  brightness  temperatures  for  the  two 
Pacific  areas.  This  figure  depicts  the  brightness  temperatures  for  Channel  4  versus  Channel  ^ 

1.  The  explanation  of  these  featmes  is  basically  the  same  as  the  previous  paragraph 
as  the  Channel  1,  mid-upper  tropospheric  peak  for  the  weighting  function,  brightness 
temperature  is  high  for  the  re^on  with  subsidence  and  is  low  for  the  region  with  an 

► 

abundance  of  water  vapor.  Conversely,  the  Channel  4  brightness  temperatures  are  high 
for  the  WP  and  and  relatively  low  for  the  EP.  SST  differences  certainly  play  a  part  in 
these  brightness  temperatme  differences,  with  higher  SSTs  (which  means  higher  brightness 
temperatures)  for  the  WP.  Additionally,  water  vapor  near  the  surface  and  cloud  effects  ^ 

also  contribute  to  the  higher  brightness  temperatures  in  the  WP  region. 

The  effect  of  ice  and  mixed  phase  and  partial  cloud  cover  clouds  cannot  be  quanti¬ 
tatively  nor  qualitatively  examined  from  these  figures  or  from  the  brightness  temperature  I 

maps.  This  is  an  area  where  future  research  is  needed. 

In  this  section  it  has  been  shown  that: 

1.  Areas  with  distinctly  Afferent  moisture  profiles  (such  as  regions  of  large  scale  subsi-  ^ 

dence  verstis  moist  ascent)  can  be  uniquely  identified  through  the  use  of  brightness 
temperature  plots. 

2.  Surface  charew:teristics  can  be  inferred  from  the  brightness  temperature  plots. 


Channel  4  Brightness  Temperature 
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Figure  2.16:  Same  as  Figure  2.12  except  Channel  4  vs.  Channel  1. 


> 


» 


Chapter  3 

MICROWAVE  RADIATIVE  TRANSFER  AND  WEIGHTING 

FUNCTIONS 

Emission  from  the  atmosphere  and  the  Hearth’s  surface  provide  the  somce  of  radia¬ 
tion  received  at  satellite  altitude  by  microwave  instruments.  Radiation  detected  by  these 
instruments  is  changed  by  the  processes  of  absorption,  emission  and  scattering  within 
the  atmosphere.  These  processes  depend  on  the  properties  of  the  Earth’s  surface,  atmo¬ 
spheric  constituents  especially  water  vapor  and  the  concentration  of  hydrometeors  and 
information  about  these  properties  can  be  retrieved  from  measurements  of  this  radiation 
(Stephens,  1993).  As  stated,  detection  of  atmospheric  water  vapor,  and  specifically  its 
vertical  profile  is  the  objective  of  the  SSM/T-2  Millimeter  Wave  Moisture  Sounder. 

For  the  frequencies  of  interest,  clear  sky  absorption  arises  primarily  from  water  vapor 
and  to  a  lesser  extent  oxygen  absorption  (Waters,  1976.)  Ozone  also  absorbs  in  this 
region,  but  to  a  much  lesser  extent  than  water  vapor  and  oxygen,  and  ozone  absorption 
has  a  negligible  effect  on  brightness  temperatures.  Cloud  droplets  also  absorb  in  the 
microwave/millimeter  wave  region. 

Figure  3.1  provides  a  physical  perspective  on  microwave  emission  and  shows  some  of 
the  effects  of  water  vapor,  s\uface  reflectance,  and  cloud  liquid  water  on  this  emission. 
Case  (2)  shows  the  emission  spectra  foimd  over  a  cold  surface  (e.g.  ocean)  versus  case  (4) 
which  presents  the  spectra  foimd  over  a  warm  surfrw:e  (e.g.  land).  The  difference  in  spectra 
between  case  1  and  case  2  shows  how  water  vapor  changes  the  brightness  temperatmes. 
These  effects  lead  to  temperature  decreases  over  land  (warm  background)and  increases 
over  ocean  (cold  background).  Similarly,  the  change  in  spectra  between  cases  2  and  3  show 
the  effects  of  cloud  droplet  (Rayleigh)  absorption.  Finally,  the  spectra  change  from  case 
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2  to  case  4  is  an  example  of  how  surface  char{K;teri8tics  may  change  intensities  measured 
by  a  satellite  radiometer. 

3.1  THE  RADIATIVE  TRANSFER  EQUATION  FOR  MICROWAVE  RE¬ 
MOTE  SENSING 

The  radiative  transfer  equation  describes  how,  in  a  mathematical  sense,  radiation 
is  transferred  from  one  layer  to  another.  The  generalized  form  of  this  equation  which 
includes  the  effects  of  both  absorption  and  emission  of  energy  is  written  as: 

dl  = -K[I  -  B]da,  (3.1) 

where  I  is  the  intensity,  k^,  is  the  voliune  absorption  coefHcient,  and  B  is  the  emission 
represented  by  Planck’s  function.  Equation  (3.1)  neglects  scattering,  which  is  common  for 
many  problems  in  microwave  radiative  transfer  and  it  also  ass\unes  the  condition  of  local 
thermodynamic  equilibriiun.  The  principles  of  emission-based  sensing  are  used  for  many 
atmospheric  applications  including  the  remote  sensing  of  sea  sxuface  temperature  (SST), 
column  water  vapor,  cloud  liquid  water  content,  precipitation  rate,  temperature  profiles, 
and  clouds. 

3.1.1  Brightness  Temperatures 

At  the  wavelength  of  microwave/millimeter  wave  radiometers,  the  Rayleigh-Jeans 
distribution  applies  —  that  is,  radiant  energy  and  temperature  are  synonymous.  This 
temperature  is  referred  to  as  the  brightness  temperature,  and,  for  satellite  radiometers  is 
a  measure  of  the  upwelling  thermal  microwave  radiation  from  the  mediiun  (in  this  case  the 
atmosphere)  below.  Using  the  Rayleigh-Jeans  approximation,  assuming  specular  reflection 
and  for  a  nadir  viewing  radiometer,  this  brightness  temperature,  Tb{v),  ran  be  expressed 


Tb{v)  =  f  exp(-r(h,  oo))y{v,  h)T{h)dh  -|-  exp(-r(0,  oo))  X 
Jo 


(1  -  R)Ttfc  +  R  f  exp(-T(0,  h))‘y{v,  h)T{h)dh  ilexp(-2r(0,  oo))Tc6,  (3.2) 
Jo 
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as: 
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DMSP  INSTRUMENTS 


FREQUENCY  (GHz) 


Figure  3.1:  Brightness  temperature  spectra  from  1-200  GHz  for  cases  (1)  tropical  atmo¬ 
sphere  with  no  water  vapor  and  reflectivity  (R)  =  0.3;  (2)  tropical  atmosphere  (CWV  = 
41.2g/m*)  and  R  =  0.3;  (3)  tropical  atmosphere  with  a  cloud  (LWC  =  0.5g/m^)  from  1-2 
km  and  R  =  0.3  and  case  (4)  same  as  case  (2)  except  R  =  0.0.  The  locations  of  the  SSM/I 
and  SSM/T-2  channels  are  indicated  by  “P  and  “2”,  respectively. 
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where  r(z,  y)  =  7(*'>  x)dx  represents  the  optical  depth  between  heights  z  and  y,  R 

is  the  reflectivity  of  the  surface  ((1-R)  is  the  emissivity  (e)  of  the  surface),  7(1/,  h)  the 
volume  absorption  coefflcient  of  the  atmosphere  at  a  height  h  for  a  frequency  v,  T(h)  the 
temperature  profile  of  the  atmosphere,  T./c  the  thermometric  temperature  of  the  surface, 
and  Td,  the  cosmic  background  (2.7°K)  temperature.  Each  term  on  the  right  hand  side  of 
(3.2)  has  a  related  physical  description.  The  first  term  is  ';he  contribution  to  the  brightness 
temperatme  due  to  upward  emission  by  the  atmosphere.  The  second  term  is  the  ‘surface 
term’  which  is  the  radiation  emitted  by  the  Earth’s  surface  which  reaches  the  satellite. 
The  third  term  is  the  downwelling  radiation  of  the  atmosphere  which  is  reflected  by  the 
Elarth’s  surface  to  the  satellite.  The  last  term  is  the  cosmic  background  temperatme  which 
is  reflected  by  the  earth’s  surface  to  the  satellite.  The  angular  dependence  of  quantities 
has  been  omitted  here  and  is  taken  to  be  understood. 

There  is  some  debate  as  to  the  necessity  of  the  cosmic  background  term.  Some 
researchers  state  that  this  term  can  be  ignored  whereas  Kakar  (1983)  argues  that  for  a 
channel  on  the  wing  of  the  absorption  line  (e.g.  such  as  the  150  GHz  chaonel  of  the  SSM/T- 
2)  that  this  term  would  add  approximately  O.S^K  to  the  measured  brightness  temperature. 
The  term  is  ret€uned  retained  in  the  radiative  transfer  model  introduced  below  and  used 
to  produce  the  results  of  subsequent  chapters. 

In  calailating  brightness  temperatures  from  (3.2),  the  state  of  the  polarization  of  ra¬ 
diation  has  not  been  considered.  The  mcdn  sources  of  polarization  arise  from  reflection  by 
the  surfaceand  by  scattering  from  irregular  particles  in  the  atmosphere.  The  latter  effects 
are  ignored  and  only  the  the  effects  of  polarization  due  to  surface  reflection  are  consid¬ 
ered.  By  ignoring  any  cross  polarization  effects  within  the  atmosphere,  (3.2)  becomes  the 
radiative  transfer  equation  for  the  vertical  state  of  polarization  assuming  R  represents  the 
siirface  reflectivity  for  vertical  polarization.  In  any  event,  the  effects  of  siuface  reflectiv¬ 
ity  are  minimal  for  Channels  1  and  2,  md  thus  the  effects  of  polarization  are  of  little 
consequence  for  the  study  described  in  this  thesis. 
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3.2  WEIGHTING  FUNCTIONS 

For  atmospheric  temperature  profiling,  (3.1)  may  be  expressed  in  the  form: 

11(00)  =  r  K{u,  h)T{h)dh  +  A,  (3.3) 

Jo 

where  A  represents  terms  independent  of  the  temperature  profile  and  K(i/4i)t  the  tempera¬ 
ture  weighting  function,  is  weakly  dependent  on  the  temperature  profile.  These  weighting 
functions  are  a  basis  for  selection  of  channels  for  temperature  profiling  instruments  and, 
in  some  cases,  a  basis  for  the  temperature  profile  rctAicVdl  itself. 

Smith  et  oL  (1979)  show  weighting  functions  for  the  TOYS  instrument  suite.  The 
MSU  operates  in  the  same  oxygen  absorption  band  as  the  SSM/T-1  instrument.  The 
depicted  weighting  functions  (shown  in  Figure  3.2)  give  a  demonstration  of  the  vertical 
resolution  of  each  of  the  channels.  The  weighting  functions  appear  in  the  integral  equations 
relating  temperature  to  the  measured  radiances  (as  in  (3.3)). 

3.2.1  Relative  Humidity-based  Weighting  Functions 

This  discussion  of  constituent  retrievals  follows  from  the  work  of  Schaerer  and 
Wilheit  (1979).  Constituent  retrievals  are  unfortunately  a  nonlinear  problem  because 
the  absorption  coefficient  appears  both  in  the  transmission  factor  (exp(—  / ~fdx))  and  in 
the  radiation  factor  'y{h)T{H)  of  (3.2).  Therefore,  the  effect  of  7(1/,  h)  at  one  height  is 
dependent  on  7(1/,  h)  at  all  other  heights  and  so  the  problem  of  constituent  retrieval  is 
also  a  non-local  one. 

Similarly,  (3.1)  can  be  expressed  as  (Schaerer  and  Wilheit  (1979)): 

n(oo)  =  r  G{hh{h)dh  -f  B,  (3.4) 

Jo 

where  B  represents  terms  independent  of  7(h);  however,  no  insight  comparable  to  the  tem¬ 
perature  profiling  example  evolves  since  G(h)  is  strongly  dependent  on  7(h).  Nonetheless, 
(3.4)  can  still  be  linearized.  A  small  change  in  the  brightness  temperature  A  Tt,(oo)  due 
to  a  small  chcinge  in  the  absorption  as  a  function  of  height,  i.e.  A  7(h),  can  be  expressed 
as: 

Ar6(oo)=  r  G[h,'r{h),T{h)]A'yih)dh. 

Jo 


(3.5) 
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After  some  manipulation  we  find  an  expression  for  G(h,  7(h),  T(h)): 
G{h,y{h),T{h))  =  exp(-r(fi,  oo))[T{h)  -  (/»)]+ 


i2exp(-T(0,oo))exp(-r(0,fi))[r(/>)  -  rfri(fi)],  (3.6) 

where  and  Tb^  are  the  upwelling  and  downwelling  brightness  temperatures  at  height 
h,  respectively.  The  weighting  fimction  with  respect  to  absolute  water  vapor,  p(h),  is: 


G,{h)  =  G{h) 

and  with  respect  to  relative  humidity,  RH(h): 


dp{h) 


where  pjat(b)  is  the  saturation  vapor  pressure  of  water  at  a  height  h.  In  this  way,  weighting 
functions  based  on  the  relative  hiunidity  profile  RH(h)  are  produced;  these  being  related 
to  brightness  temperature  changes  by 


ArbT=  r  GnHih)/iRHih)dh. 
Jo 


3.2.2  Transmission-based  Weighting  Functions 

Transmission-based  weighting  functions  assmne  a  different  yet  comparable  form  and 
are  similar  to  those  used  for  temperatme  retrievals.  Isaacs  and  Deblonde  (1987)  used  these 
weighting  functions  for  their  study  of  cloud  effects  on  183  GHz  retrievals.  These  weighting 
functions  are  also  the  same  as  those  used  for  infirared  measurement  techniques  where  both 
vertical  profiles  of  water  vapor  and  total  precipitable  water  are  derived  (Prabhalmra  and 
Dalu,  1980). 

Define  T(r,p)  =  exp(— r/p)  as  the  transmittance  of  the  atmosphere.  Neglecting 
the  dependence  on  the  zenith  angle  (as  above),  the  derivative  of  the  transmittance  be¬ 
comes:  dT  =  exp(-r)dr  or  dT  =  exp(-r)  {—'ydh).  Using  the  Rayleigh-Jeans  limit,  (3.2) 
becomes: 

rbT(oo)  =  T(r„)T./e(l  -  R)  +^°T(r)dr(r)  +r(r)i?^  T(r)dr(r)  +T^{to)RTa{3.10) 
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where  is  the  optical  depth  from  the  satellite  to  the  surface.  By  grouping  common  factors 
together,  liquation  3.10  can  be  simplified  to: 

TiT(oo)  =  [r./c(l  -R)  +  RTirM  r{To)  +  [1  +  r(r)iJ]  r(r)dr(r)  (3.11) 


where  there  is  now  only  one  integral  term. 

Defining  a  weighting  function  W(2)  for  a  vertical  coordinate  Z  as  W{Z)  =  = 

T{Z)‘y{2)  then  the  previous  equation  can  be  written: 


TfcTC^oo)  =  [T{Zo){l  -  R)  +Rr{Zo)T^]T{Zo)  +[l+T{Z)R]J^°^T{Z)W{Z)dZ.{3.12) 


Weighting  functions  derived  in  this  manner  are  convenient  since  they  closely  resemble 
the  weighting  functions  often  used  for  temperature  retrievals.  Transmission  is  a  quantity 
easily  calculated  in  a  propagation  model,  and  these  models  permit  the  calculation  of  the 
derivative  of  transmission  with  height  and  hence,  the  weighting  ftmctions  themselves. 

As  an  example  of  transmission-based  temperature  weighting  functions.  Figure  3.2 
shows  weighting  functions  plotted  for  channds  1-4  of  the  SSM/T-1  temperature  sounder  for 
a  tropical  atmosphere  with  a  thermometric  temperature  of  300  K.  Figure  3.2  is  comparable 
to  Figure  2,  Liou  et  oL  (1979)  and  the  MSU  channels  of  Figure  1,  Smith  et  al.  (1979). 

Weighting  functions  from  the  two  different  methods  appear  very  similar  when  plotted. 
Plots  of  the  transmission-based  weighting  functions  are  for  the  SSM/T-2  are  discussed  later 
in  Chapter  4. 


3.3  ATMOSPHERIC  PROFILES,  PROPAGATION  MODEL,  AND  CLOUD 
MODELS 


To  simulate  microwave  brightness  temperature  requires  the  use  of  models  to  represent 
various  aspects  of  the  atmosphere  and  the  propagation  of  microwave  energy  through  the 
atmosphere.  We  use  simulations  of  the  brightness  temperatures  to  understand  how  they 
vary  with  changing  atmospheric  properties  and  to  study  these  effects  in  isolation.  The 
selection  of  the  models  described  below  was  made  for  a  variety  of  reasons  including  the 
use  of  a  models)  in  similar  studies  done  previously,  the  most  current  model  available,  the 
computational  efficiency  of  the  model,  and  its  simplicity. 
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Figure  3.2:  Clear  sky  weigliiing  functions  for  the  SSM/T-1  channels  1-4  for  a  tropical 
atmosphere,  reflectivity  =  0.03,  and  a  column  water  vapcr  (CWV)  =  41.2  kg/m^. 
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3.3.1  Atmospheric  Profiles 

McClatchey  et  al.  (1972)  standard  atmosphere  profiles  of  mid-latitude  winter,  I 

mid-latitude  summer,  and  tropical  atmospheres  were  used  as  the  reference  for  height  (km), 
pressure  (mb  or  hPa),  temperature  (°K  or  K),  density  (gm~®)  and  water  vapor  density 
(gm~^).  The  radiative  transfer  code  also  included  a  featme  which  would  allow  the  water  ^ 

vapor  density  to  be  increased  or  decreased  by  a  common  factor  throughout  the  atmosphere. 

This  is  used  to  increase  the  ntunber  of  different  atmospheres  for  the  simulations  discussed 
in  Chapter  5. 

► 

3.3.2  The  Propagation  Model 

The  Liebe  millimeter  wave  propagation  model  (MPM92)  (which  is  an  updated  version 
of  MPM89  (Liebe,  1989)  and  Liebe,  personal  communication)  allows  the  specification  of  an  ^ 

input  atmospheric  profile.  The  atmospheric  profile  provides  the  necessary  meteorological 
parameters  to  allow  the  computation  of  the  volume  absorption  coefficient,  7.  Knowledge 
of  the  voliuue  absorption  coefficient  then  allows  the  computation  of  the  optical  depth,  r.  ^ 

From  the  optical  depth  the  calculation  of  the  transmittance  is  straightforward  and  the 
transmission-based  weighting  function  follows  as  described  above. 

The  MPM92  model  accoimts  for  dry  air  attenuation,  oxj'gen  and  water  vapor  line 

► 

absorption,  water  vapor  continuum,  emd  Rayleigh  cloud  droplet  absorption.  The  Zeeman 

effect  for  oxygen  for  pressures  less  than  0.7  kPa  and  Doppler  broadening  for  water  vapor 

are  also  included.  One  distinct  advantage  of  MPM92  is  the  small  number  of  lines  used 

in  the  model.  The  model  uses  44  O2  and  30  H2O  local  line  contributions  versus  the  over  ^ 

500  contributing  spectral  lines  for  fi-equencies  below  1  THz.  The  small  munber  of  lines 

in  MPM92  require  considerably  less  computation  time  than  do  other  codes  (i.e.  RAD- 

TRAN).  Versions  of  MPM  have  been  used  by  various  researchers  (e.g.  W€uig  et  al.  (1992)  > 

and  Fleming  et  al.  (1991))  to  study  atmospheric  propagation  at  millimeter/microwave 

wavelengths. 
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3.3.3  Cloud  Models 

To  simulate  cloud  effects  within  the  field  of  view  (FOV)  of  the  radiometer,  a  set  of  * 

cloud  models  was  chosen  to  be  incorporated  within  the  atmospheric  attenuation  profiles. 

The  cloud  models  chosen  (after  Isaacs  and  Deblonde,  1987)  are  firom  the  Air  Force  Geo¬ 
physics  Laboratory  FASCODE  model  (Falcone  et  aL,  1979).  Cloud  paxcuneters  for  the  five  I 

cloud  models  are  ^ven  in  Table  3.1. 


Table  3.1:  Cloud  Type  Characteristics  (after  Isaacs  and  Deblonde,  1987) 


Cloud 

Mode! 

Type 

Mode 

Radius, 

fun 

Liqmd 

Water 

Content, 

gmrZ 

Vertical 

Extent, 

km 

Stratus 

2.7 

0.15 

0.5-2.0 

Cumulus 

6.0 

1.00 

1.0-3.5 

Altostratus 

4.5 

0.40 

2.5-3.0 

Stratocumulus 

6.25 

0.55 

0.5-1.0 

5 

Nimbostrattis 

3.0 

0.61 

0.5-2.5 

Cloud  attenuation  was  calculated  using  the  Rayleigh  absorption  approximation  of 
Mie’s  forward  scattering  function  (Liebe,  1989).  Cloud  liqmd  water  content  {w,  gm~^)  and 
fi-equeiicy  are  the  required  variables.  Cloud  attenuation  (calculated  by  Liebe’s  MPM92 
model)  was  added  to  the  clear  sky  absorption  for  each  layer  of  the  atmosphere  for  the 
entire  cloud  layer  tis  specified  by  the  cloud  model.  The  MPM  model  assumes  that  the 
cloud  fills  the  entire  FOV  and  that  there  is  only  one  cloud  layer.  These  assumptions  will 
tend  to  overestimate  and  imderestimate  the  effects  of  clouds,  respectively. 


I 


Chapter  4 

SENSITIVITY  EFFECTS  FOR  THE  DMSP  SSM/T-2  ' 

This  chapter  is  devoted  primarily  to  the  discussion  of  forward  radiative  transfer  sim¬ 
ulations  of  the  SSM/T-2  brightness  temperatures  and  how  these  temperatures  respond  to  t 

a  variety  of  different  atmospheric  parameters.  The  approach  is  to  build  an  understanding 
of  this  radiative  transfer  problem  before  tackling  the  inverse  problem  —  that  of  retrieving 
moisture  information  from  brightness  temperature  measurements.  There  are  a  number  ^ 

of  factors  which  influence  brightness  temperatures  and  hence,  the  retrieval  of  moisture 
profiles  using  3  channels  near  183  GHz  and  a  window  channel  near  92  and  a  line  wing 
channel  at  150  GHz.  These  factors  will  be  referred  to  as  sensitivity  factors  and  include: 
the  thermometric  temperature  and  its  vertical  profile,  the  surface  reflectivity  (emissivity), 
the  hxuuidity  (moisture)  profile,  the  characteristics  of  clouds  (both  water  clouds  and  per¬ 
haps  ice  clouds),  the  sinface  elevation,  the  channel  noise,  and  scan  angle.  The  influence 
of  some  of  these  factors  on  T-2  simulated  brightness  temperatures  is  now  examined.  ^ 

4.1  DEPENDENCE  ON  MOISTURE  PROFILES 

For  the  simulation  studies  described  below,  three  basic  atmospheric  profiles  are  used:  > 

a  tropical  (T)  profile,  a  midlatitude  (MLS)  summer  profile,  and  a  midlatitude  winter 
(MLW)  profile  (McClatchey  et  al.,  1972).  The  water  vapor  density  data  for  each  profile 
will  also  be  multiplied  by  a  constant  factor  at  each  level  to  provide  a  method  of  varying  ^ 

the  coliunn  water  vapor  (CWV)  amour*  ce  each  moisture  profile  is  clmracterized 
by  a  different  column  water  vapor  amov  a.  since  each  atmosphere  has  a  different 
temperature  profile,  the  simulated  chai  cel .  .  '  T-2  will  respond  differently  to  each 

► 

profile. 


> 
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The  clear  sky  weighting  function  profiles  for  an  ocean  background  with  an  emissivily 
of  0,7  are  shown  in  Figures  4.1-4.3  with  altitude  as  the  ordinate  and  the  weighting  function 
^  (or  ^  as  designated  in  the  figures)  as  the  abscissa.  These  profiles  are  plotted  for  the 
nadir  angle  and  for  values  the  column  water  vapor  CWV  (kg  m~*)  as  given. 

Comparing  these  figures  it  can  easily  be  seen  how  the  amoimt  of  moisture  in  the 
atmosphere  affects  the  location  of  the  peaks  of  the  weighting  functions.  For  example,  for 
the  tropical  and  midlatitude  summer  atmospheres  the  182.31  GHz  clumnel  peaks  between 
7  and  8  km  whereas  for  the  midlatitude  winter  atmosphere  the  channel  peaks  between  5 
and  6  km.  Similar  effects  are  seen  for  the  other  rhannpla  where  the  more  moishire  in  the 
atmosphere  the  higher  in  the  atmosphere  the  channel  peaks.  Even  the  ‘window’  channel 
(91.655  GHz)  sees  the  effects  of  moisture  in  the  atmosphere  and  this  serves  as  a  reminder 
that  the  ‘window’  is  not  entirely  transparent. 

It  can  also  be  seen  from  Figures  4.1-4.3  the  broadness  of  the  weighting  functions  as 
depicted  for  these  cases.  This  broadness  is  an  obvious  limitation  of  the  system  in  that  the 
resolution  for  retrievals  is  limited  by  the  lack  of  sharpness  of  these  weighting  functions. 
The  resolution  of  the  T-2  was  expected  to  be  on  the  order  of  2-3  km  based  on  the  weighting 
functions  of  this  type  (e.g.  Schaerer  and  Wilheit,  1979). 

Additionally,  from  Figures  4. 1—4.3  it  is  relatively  easy  to  develop  some  intuition  as 
to  the  piupose  of  each  channel,  at  least  for  clear-sky  conditions.  The  two  cliannels  clos¬ 
est  to  the  183  GHz  line  relate  to  the  measurement  upper  tropospheric  moisture;  the  line 
wing  channel  (150  GHz)  to  low  level  moistiue;  and  the  window  channel  senses  low  level 
moisture  and  surface  characteristics,  including  the  emissivity  and  perhaps  the  thermomet¬ 
ric  (‘skin’)  temperatme.  The  remaining  channel  is  relevant  to  midtropospheric  moisture 
determinations. 

As  a  final  examination  of  how  moisture  changes  the  weighting  functions  of  the  SSM/T- 
2,  the  midlatitude  summer  profile  will  again  be  used.  In  this  case,  the  specific  hmnidity 
is  increased  by  a  factor  of  two  throughout  the  entire  atmosphere  (20  km).  Figure  4.4, 
compared  to  Figure  4.2,  shows  how  the  increased  moisture  at  all  levels  moisture  changes 
the  weighting  functions  significantly. 
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Figure  4.2:  The  same  as  Figure  4.1  except  for  a  midlatitude  sxunmer  atmosphere. 
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Figure  4.4:  The  same  as  Figure  4.2  except  with  2  X  water  vapor  profile. 
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These  fotir  figures  grt^hically  show  how  the  moisture  profiles  affect  the  weighting 
functions  when  it  is  plotted  against  height  (km).  Certainly  there  is  a  comprehensive 
change  in  the  weighting  functions  when  the  atmospheric  water  vapor  is  doubled.  This 
goes  back  to  the  discussion  earlier  about  the  nonlinearities  in  the  process  caused  by  the 
absorption  factor. 

Another  way  of  grt^hicaUy  showing  where  the  atmospheric  constituoits  diminish  the  * 

emission  of  the  atmosphere  and  siuface  is  to  plot  transmission  versus  altitude.  Clear  sky 
atmosphere  plots  of  the  transmission  corresponding  to  Figures  4.1-4.3  are  shown  in  Figures 
4.5-4.7.  Noteworthy  is  the  generally  similar  appearance  of  each  of  the  transmission  curves  > 

and  how  the  levels  of  the  most  rapid  decrease  in  transmission  with  decreasing  altitude 
correspond  to  the  peaks  of  the  weighting  functions. 

4.2  DEPENDENCE  ON  BACKGROUND  ► 

The  surface  reflectivity  R  (or  emissivify  £  =  1-R)  depends  on  the  characteristics  of 
the  siirface  of  the  earth  as  viewed  by  the  satellite  radiometer.  Land  and  coastal  regions 
have  large  variations  in  surface  emissivity  and  hence,  the  retrieval  of  moisture  profiles  can 
be  expected  to  suffer  degradation  for  these  regions  as  compared  with  ocean  surfaces. 

4.2.1  Ocean  Background 

The  emissivity  of  a  smooth  ocean  surface  can  be  calculated  from  the  dielectric  con- 
steint  of  water.  Calculations  using  the  dielectric  data  of  Ray  (1972)  result  in  a  value  of 
~0.7  for  the  emissivity  of  calm  ocean  water  at  183  GHz.  Lutz  et  al.  (1990)  have  devel- 
oped  an  algorithm  (from  Wilheit,  1990)  with  an  iterative  procedme  for  determining  the 
emissivity  for  wind-roughened  ocean  surfaces  or  land  smdaces.  The  effect  of  wind  rough¬ 
ening  procedure  degrades  the  ability  to  retrieve  moisriire  profiles  (versus  retrievals  over  a 
constant  background).  * 

The  low  emissivity  of  ocean  smfaces  is  one  important  aspect  which  helps  with  the 
retrieval  of  low  level  moisture.  This  is  because  the  rsuiiometer  senses  emission  from  the 
warm  atmosphere  in  contrast  to  a  cold  background  which  results  in  a  large  brightness  ► 

temperature  contrast.  Additionally,  the  emissivity  of  a  water  siuface  is  thoiight  to  be 
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Figure  4.5:  Clear  sky  transmission  for  the  tropical  atmosphere. 
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Figure  4.6:  Clear  sky  transmission  for  the  midlatitude  atmosphere. 
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Figure  4.7:  Clear  sky  transmission  for  the  midlatitude  winter  atmosphere. 


rather  stable,  especially  when  compared  to  land  surfaces.  One  problem  with  the  emissivity 
of  the  ocean  is  that  it  is  a  weak  function  of  the  smdace  wind  speed  above  it.  Various 
models  (e.g.  WUheit,  1979  and  Pandey  and  Kakar,  1982)  were  developed  for  determining 
the  emissivity  of  a  wind-roughened  ocean  surface  for  earlier  instruments  which  operated 
at  lower  frequencies  (e.g.  SMMR  and  SSM/I).  Currently,  there  is  a  dearth  of  models  which 
are  applicable  for  the  region  near  183  GELz. 

To  examine  the  effect  of  the  emissivity  of  the  ocean’s  surface,  the  emissivity  of  an 
ocean  surface  underlying  a  clear  sky  midlatitude  winter  atmosphere  was  varied  between 
0.82  and  0.64.  These  results  are  shown  in  Table  4.1. 


Table  4.1:  Brightness  temperatures  (K)  for  emissivities  firom  0.82  to  0.66  (reflectivities 
from  0.34  to  0.18)  for  a  midlatitude  winter  atmosphere. 


Midlatitude  Winter  Brightness  Temperatures 

FREQ. 

(GHz) 

Reflectivity 

0.34 

0.32 

0.30 

0.28 

0.26 

0.24 

0.22 

0.20 

0.18 

180.31 

254.99 

255.06 

255.13 

255.20 

255.27 

255.34 

255.41 

255.48 

255.55 

182.31 

245.77 

245.77 

245.77 

245.77 

245.77 

245.77 

245.77 

245.77 

245.77 

176.31 

243.01 

244.40 

245.80 

247.19 

248.59 

249.98 

251.37 

252.77 

254.16 

150.00 

203.81 

208.08 

212.35 

216.62 

220.90 

225.17 

229.44 

233.72 

237.99 

91.655 

193.83 

198.78 

203.74 

208.69 

213.64 

218.60 

223.55 

228.50 

233.46 

The  channel  closest  to  the  absorption  line  (182.31  GHz)  has  the  same  brightness 
temperatme  for  all  emissivities  since  the  atmosphere  is  so  opaque  that  microwave  thermal 
energy  emitted  and  reflected  by  the  sxuface  does  not  reach  the  satellite  (i.e.  the  surface 
term  has  no  contribution).  The  180.31  GHz  channel  is  also  basically  unaffected  by  the 
surface  reflectivity.  The  line  wing  channel  (150  GHz)  and  the  window  channel  (91.655 
GEbs)  brightness  temperatvues  change  appreciably  as  the  emissivity  changes  according  to 
the  values  depicted  in  the  table.  This  result  will  be  discussed  later  with  respect  to  the 
operational  algorithm  in  use. 

The  root  mean  square  (RMS)  error  calculated  pre-launch  performance  of  the  T-2  with 
the  a  priori  regression  coefRcients  using  the  multiple  linear  regression  technique  employed 
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operationally  showed  the  best  retrieval  statistics  occurrs  over  ocean  backgrounds  with 
their  relatively  uniform  surface  emissivity  (Boucher  et  aL,  1993). 

4.2.2  Land  Background 

The  emissivity  of  dry  land  backgrounds  is  thousht  to  vary  from  0.8  to  1.0  over  the 
range  of  frequencies  of  140  to  200  GHz  (Schaerer  and  Wilheit,  1979).  These  early  values  ) 

should  be  compared  with  those  of  Wang  et  al.  (1992)  below.  As  mentioned  above,  the 
Wilheit  (1990)  algorithm,  as  amended  by  Lutz  et  al.,  retrieves  the  land  emissivity  by  an 
iterettive  process  as  part  of  its  computations.  ^ 

The  presence  of  vegetation  cover,  however,  can  change  the  reflectivity  of  land  surfaces 
as  vegetation  generally  is  an  absorptive  medium  and  therefore  decreases  reflectivity  values. 

Frozen  soils  also  generally  have  low  reflectivity  values  (~  0.1).  Soils  with  moisture  contents 
of  less  than  50%  have  dielectric  constants  where  the  real  part  is  less  than  6  which  yields 
a  value  of  R  of  less  than  0.15  (Wang  and  Schmugge,  1990). 

Snow  cover  presents  a  different  challenge.  Wang  et  al.  (1992),  using  the  AMMS 
described  in  Chapter  1,  observed  emissivities  as  low  as  0.58  when  making  flights  over 
Alaska.  They  siirmised  that  the  reason  for  such  low  emissivities  is  due  to  scattering  by 
snow  particles  on  the  groimd.  Wang  et  oL  found  the  emissivities  on  these  flights  generally 
ranged  from  0.6-0.9  which  differ  from  the  estimates  of  Schaerer  and  Wilheit  of  more  than  ^ 

a  decade  ago. 

Conceptually,  low  level  moistiu«  retrievals  would  be  expected  to  be  more  difflcult  over 
land  because  of  the  high  and  variable  emissivity  of  this  siirface.  This  presents  the  possible  I 

situation  where  emission  from  the  lower  level(s)  of  the  atmosphere  and  the  emission  from 
earth’s  surface  may  be  nearly  at 

the  same  brightness  temperature.  Isaacs  and  Deblonde  (1985)  noted  that  surface  ^ 

emissivity  affects  the  retrieval  accuracy  at  all  levels  and  not  just  at  the  layers  near  the 
surface  based  on  the  statistical  retrieval  method  they  xised. 

Brightness  temperatmes  for  an  ocean  sxirface  with  an  emissivity  of  0.7  and  for  a  land 

I 

surface  with  an  emissivity  of  1.0  are  compared  in  Table  4.2. 
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The  three  channels  nearest  the  absorption  line  (182.31  GHz  channel  not  shown)  are 
virtually  unaffected  by  the  surface  emissivity.  However,  this  is  certainly  not  the  case 
for  the  150  and  91.655  GHz  channels.  Radically  different  brightness  temperatures  are 
computed  for  cases  with  emissivities  of  1.0  and  0.7  for  these  channels.  This  is  simply  the 
factor  of  a  warm  background  and  a  cold  backgroimd,  respectively. 

As  mentioned  above,  for  the  ocean  background  case  in  Table  4.1  brightness  tempera¬ 
tures  for  the  150  and  91.655  GHz  changed  appreciably  for  small  emissitrity  changes.  The 
operational  algorithm  retrieves  the  surface  emissivity  using  the  91.655  GHz  channel  emd 
a  physical  retrieval  model  (Grifhn  et  cd.,  1993). 

4.3  DEPENDENCE  ON  CLOVJDS 

Water  clouds  affect  microwave  radiation  primarily  through  the  processes  of  absorption 
and  emission.  As  discussed  in  the  previous  chapter,  Rayleigh  absorption  is  assumed  by 
the  cloud  model  (MPM92)  being  used.  The  process  of  scattering  is  usually  considered 
negligible  since  the  single  scatter  albedo,  w,  at  the  microwave/millimeter  wavelengths  is 
very  small. 

To  depict  the  effects  of  clouds  on  weighting  functions.  Figures  4.8-4.12  show  the  five 
cloud  models  for  a  tropical  atmosphere  over  an  ocean  (R  =  0.3  or  e  =  "  7)  backgrotmd.  For 
the  cloud  models  rised,  and  from  the  study  of  the  figures  shown,  the  following  conclusions 
are  deduced: 

1.  The  channel  closest  to  the  absorption  line  does  not  ‘see’  any  of  the  low  level  water 
clouds. 

2.  The  180.31  GHz  channel  sees  onlj'  those  clouds  with  cloud  tops  above  2  km. 

3.  The  weighting  functions  for  the  other  three  channels  are  radically  affected  by  clouds 
in  the  FOV. 

4.  The  weighting  functions  peak  sharply  at  the  clovd  top. 


Height  (km) 


Weight  (dTr/dz) 
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Figure  4.8:  Weighting  functions  as  in  Figure  4.1  except  for  stratus  cloud  {w=  O.lSgro"*'’) 
from  0.5-2.0  km  (cloud  model  1). 
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Figure  4.9:  The  same  as  Figiu^  4,8  except  for  cumulus  cloud  (w  =  1.0  gm  ’)  from  1.0-3.5 
km  (cloud  model  2}. 
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Figure  4.11:  The  same  as  Figure  4.8  except  for  stratocumulus  cloud  {w  =  0.55  gm  from 
0.5-1.0  km  (cloud  model  4). 
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Figure  4.12:  The  same  as  Figure  4.8  except  for  nimbostratus  cloud  {w  —  0.61  gm~®)  from 
0.5-2.5  km  (cloud  model  5). 
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4.3.1  Dependence  on  Layer  Thickness 

The  microwave  properties  of  clouds  in  the  Rayleigh  absoi^ptlon  re^me  are  defined  by  t 

several  factors.  One  factor  is  the  thickness  of  the  cloud  layer.  This  is  because  for  a  very 
thin  cloud,  much  of  the  microwave  energy  from  below  is  transmitted  through  the  cloud 
layer.  To  determine  the  effect  of  varying  cloud  thicknesses  on  the  brightness  temperatures 
at  satellite  altitude,  clouds  were  inserted  into  the  tropical  atmosphere  with  a  constant 
cloud  base  of  1  km.  Cloud  thicknesses  were  varied  from  1  to  5  km  by  1  km  increments. 

Table  4.3  shows  simulated  brightness  temperatures  for  clouds  of  variable  thickness  and 

liquid  water  content.  According  to  the  simulations  (svimmarised  in  this  table),  the  15C 

GHz  channel  and  to  a  lesser  extent  the  92  GHz  are  the  most  responsive  to  the  presence  of 

clouds.  The  responsiveness  in  the  150  GHz  channel  is  expected  since  it  is  designed  to  peak 

near  the  surface  as  depicted  in  Figure  4.1.  Note  that  brightness  temperatmes  decrease  for  * 

increasing  LWC  and  increasing  cloud  thickness  {A  z)  with  a  land  backgroimd.  The  reason 

for  this  is  that  the  clouds  diminish  the  effects  of  the  surface  term. 

4.3.2  Dependence  on  Layer  Location  ^ 

It  has  been  seen  that  the  weighting  functions  for  atmospheres  containing  clouds  tend 
to  peak  near  the  cloud  top.  This  gives  a  good  indication  cf  the  effects  of  where  a  cloud  is 
located  in  the  atmosphere  and  Its  influence  on  brightness  temperatures.  * 

To  determine  the  influence  of  the  layer  location  on  brightness  temperature.  Table  4.4 
shows  the  simulated  brightness  temperatures  for  a  tropical  atmosphere  with  a  constant 
cloud  thickness,  reflectivity,  and  a  varying  doud  base  and  liqmd  water  content.  The  data  > 

in  this  table  indicates  how  the  brightness  temperatures  decrease  as  a  doud  is  moved  higher 
into  the  atmosphere.  This  is  expected  since  the  temperature  profile  decreases  with  height 
and  the  doud  limits  the  microwave  emissions  from  below  it.  , 

f 

4.3.3  Dependence  on  Cloud  Liquid  Water  Content 

Figures  4.13-4  16  depict  brightness  temperatures  for  the  five  T-2  channels  for  several 
doud  models  with  the  LWC  varied  from  0.0  to  G.50  gm~®.  The  brightness  temperatuies  ^ 

for  the  3  channels  near  183  GHz  are  basically  constant  whereas  for  the  150  and  91.655 


> 
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Table  4.2:  Comparison  of  brightness  temperatures  for  land  (e  =  1.0)  and  ocean(e  =  0.7) 
smTaces. 


Brightness  Temperature 


Atmospheric 

case 

180.31  GHz 

176.3] 

GHz 

150.00  GHz 

91.655  GHz 

Land 

Ocean 

Land 

Ocean 

Land 

Ocean 

Land 

Ocean 

Clear  tropical 

264.71 

264.71 

277.58 

277.36 

291.19 

265.12 

295.95 

237.72 

Clear  ML  sum. 

263.54 

263.54 

276.32 

275.08 

290.91 

252.62 

295.84 

229.72 

Clear  ML  winter 

256.44 

255.31 

249.34 

291.92 

223.19 

295.93 

216.26 

Cloud  model  1 

256.89 

256.89 

wmm 

269.81 

284.18 

272.76 

292.04 

248.02 

Cloud  model  2 

264.67 

264.67 

276.89 

276.87 

289.13 

277.67 

294.28 

253.08 

Cloud  model  3 

262.32 

262.32 

271.23 

271.23 

278.64 

278.51 

283.51 

279.32 

Cloud  model  4 

263.92 

263.92 

275.58 

275.52 

286.91 

273.42 

292.59 

250.31 

Cloud  model  5 

264.71 

264.71 

277.22 

277.18 

289.55 

278.92 

294.77 

255.77 

Table  4.3:  Brightness  temperatures  for  variable  cloud  thickness  (A  z  (km))  and  LWC 
(gm~®)  for  a  fixed  cloud  base  (1  km),  emissivity  (0.97),  and  a  tropical  atmosphere. 


br; 

[GHTNESS  TEMPERATURES  (K) 

Cloud 

A  z 

FREQUENCY  (GHz) 

LWC 

(km) 

180.31 

176.31 

1.5P.00 

91.655 

Cloud 

X 

264.67 

276.92 

287.85 

289.91 

Model 

2 

264.08 

286.20 

288.91 

1 

3 

263.11 

274.26 

283.54 

287.14 

LWC  = 

4 

261.99 

271.69 

280.13 

284.67 

0.15 

268.58 

276.24 

281.69 

Cloud 

1 

264.62 

275.24 

284.30 

288.59 

2 

263.42 

272.68 

280.90'^ 

285.34 

2 

3 

27.5.91 

280.37 

LWC  = 

4 

258.78 

265.44 

269.86 

274.15 

1.00 

5 

mmitm 

267  92 

Cloud 

1 

264.65 

^6.50 

286.97 

289.59 

Model 

2 

263.90 

275.16 

284.43 

287.69 

3 

0 

262.63 

272.95 

280.49 

284,44 

LWC  = 

4 

261.17 

269.61 

275.41 

0.41) 

5 

259.13 

265.53 

269.92  1  274.90 
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Table  4.4:  Brightness  temperatures  for  variable  cloud  base  and  LWC  for  a  fixed  emissivity 
(e  =  1.0)  and  atmosphere  (tropical). 


BRIGHI 

fNESS  TEMPERATURE  (K) 

Cloud 
Base  (km) 

LWC 

(gm'®) 

FREQUENCY  (GHz) 

182.31 

180.31 

176.31 

150.00 

91.655 

1 

0.15 

251.24 

264.67 

276.92 

289.34 

294.49 

1.5 

0.15 

251.24 

264.50 

276.57 

288.69 

294.02 

2 

0.15 

251.23 

264.08 

276.14 

287.99 

293.55 

3 

0.15 

251.11 

263.13 

274.65 

286.24 

294.43 

5 

0.15 

249.51 

260.88 

271.16 

283.08 

290.25 

1 

^21311 

251.24 

264.62 

275.24 

284.54 

290.43 

1.5 

1.00 

251.24 

264.26 

273.99 

282.95 

288.81 

2 

1.00 

251.23 

263.42 

272.73 

281.44 

287.32 

3 

1.00 

251.04 

261.25 

269.71 

276.77 

283.11 

5 

1.00 

247.74 

255.82 

261.20 

265.86 

273.27 

GHz  channels  brightness  temperatures  vary  with  liquid  water  content,  cloud  location,  and 
cloud  thickness. 

Figures  4.15  and  4.16  are  identical  except  for  the  emissivity  (reflectivity).  The  latter  > 

figure  shows  the  effects  of  the  clouds  by  decreasing  the  91.655  and  150  GHz  brightness 
temperatmes  over  a  300  K  surface  as  the  cloud  liquid  water  content  is  increased.  Contrast 
this  to  Figure  4.15  where  the  brightness  temperature  for  a  low  emissivity  surface  (cold  ^ 

background  at  210  K,  as  discussed  above).  In  this  case,  the  brightness  temperatures  of 
the  window  channeb  systematically  increases  as  the  liquid  water  path  (LWP)  is  increased. 

This  yields  the  observation  that  clouds  generally  increase  brightness  temperatures 
ove’  oceans  and  decrease  brightness  temperatures  over  land.  These  effects,  however,  ap¬ 
pear  only  in  the  wrindow  channels  and  to  a  lesser  extent  in  the  183  ±  7  GHZ  channel, 
depending  on  cloud  altitude.  This  response  to  clouds  can  also  be  seen  in  Table  4.2.  The 
near  linear  response  of  the  window  (especially)  and  line  wing  channel  (€is  shown  in  Figure  ^ 

4.13)  to  changes  in  LWC  offers  a  simple  way  of  retrieving  this  cloud  property  by  using 
these  channels.  Several  authors  have  considered  this  (including  Isaacs  et  al.,  1985). 


BRIGHTNLrfS  TEMPERATURE  (K) 


69 


0.0  0.05  0.1  0.15  0.2  0.25  0.3  0.35  0.4  0.45  0.5 

Cloud  LWC  (g/m^) 


Figure  4.13:  Brightness  temperatures  for  Cloud  Model  1  and  a  tropical  atmosphere  with 
LWC  from  0.01  to  0.50  gm“*. 
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Figure  4.14;  Brightness  temperatures  as  in  Figure  4.13  except  for  Cloud  Model  2. 
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Figure  4.15:  Brightness  temperatures  as  in  Figure  4.13  except  for  Cloud  Model  3. 
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Figure  4.16:  Brightness  temperatures  as  in  Figure  4.15  except  for  R  =  0.0. 


4.4  SCAN  ANGLE  EFFECTS 


The  effects  of  scan  angle  on  brightness  temperatures  can  be  considered  through  the  Tise 
of  MPM92.  Scan  angle  affects  brightness  temperatures  since  it  changes  the  optical  path. 
As  shown  in  Figure  3.1,  the  scan  angle  (or  scene  station  view  angle  as  in  Figure  3.1)  on  the 
SSM/T-2  varies  between  1.5°  and  40.5°  as  the  sensor  scans  across  a  line.  Unfortunately, 
there  is  no  easy  way  to  correct  for  this  change  in  brightness  temperatiu*e  for  the  change 
in  the  atmospheric  path  which  is  viewed.  Figures  4.17  and  4.18  show  that  the  brightness 
temperature  can  both  increase  or  decrease  as  the  sensor  moves  away  from  nadir.  There 
are  competing  effects  here;  limb  darkening,  occurs  when  the  sensor  moves  away  from  nadir 
and  the  path  length  increases  and  the  weighting  functions  move  upwards  vertically  in  the 
atmosphere.  Limb  brightening  iccui^  for  the  wing  and  window  channels  which  sense  near 
the  surface.  The  temperature  of  a  cold  surface  is  replaced  by  higher  temperatures  as  the 
path  is  extended  by  viewing  away  from  nadir.  This  complex  problem  would  likely  be 
handled  operationally  through  the  use  of  regression  equations.  For  example,  regression 
equations  are  used  to  solve  the  scan  angle  <hfferences  for  the  MSU  (Grody,  1983). 

4.5  INTERPRETATION  OF  RESULTS 

The  SSM/T-2  responds  to  a  variety  of  atmospheric  phenomenon  in  a  manner  much 
different  than  those  instruments  which  preceded  it  in  the  remote  sensing  of  the  atmosphere. 
In  this  chapter  it  has  been  shown  that: 

1.  Weighting  functions  and  brightness  temperatures  are  dependent  on  the  moistxue  and 
temperature  profile  of  the  atmosphere. 

2.  The  backgroimd  and  its  emissivity  have  a  significant  effect  on  the  brightness  tem¬ 
peratures  of  the  window  channels. 

3.  Backgrounds  with  higher  emissivities  result  in  higher  brightness  temperatures. 

4.  Liquid  we*^**!*  clouds  that  form  low  in  the  atmosphere  have  a  profound  influence  on 
brightness  temperatures  in  the  window  chaimels. 
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Figure  4.17;  Brightness  temperatures  as  a  function  of  observation  angle  for  a  tropical 
atmosphere. 
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Figure  4.18;  The  same  as  Figure  4.17  except  for  a  midlatitude  summer  atmosphere  and 
T,/c  =  285  K. 
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(a)  Increasing  cloud  LWC  decreases  brightness  temperatxires. 

(b)  Higher  clouds  correspond  to  lower  brightness  temperatm-es. 

(c)  Liquid  water  content  affects  the  91  and  150  GHZ  channels  brightness  temper¬ 
atures  greatly. 

(d)  Clouds  generally  increase  brightness  temperatmes  over  ocean  stirfaces  and  de¬ 
crease  brightness  temperatures  ov»  land. 

Certainly  water  clouds  affect  brightness  temperatures  in  a  systematic  manner.  Clouds 
become  more  difficult  to  quantify  when  one  considers  the  possibility  of  mixed  phase  clouds, 
ice  clouds,  and  partial  cloud  cover. 


Chapter  5 
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WATER  VAPOR  BURDEN 

In  this  cluster  a  simple  method  is  described  for  deriving  the  column  vapor  bur¬ 
den  fix>m  brightness  temperatures.  It  will  be  demonstrated  how  the  measured  brightness  ^ 

temperatures,  especially  these  of  183  ±  1  GHz  (Channel  2)  and  183  ±  3  GHz  (Channel 
1),  correspond  to  a  level  of  constant  water  vapor  overburden.  This  result  together  with 
temperature-pressure  data  are  used  to  establish  the  burden  at  a  constant  pressure.  The  I 

method  is  demonstrated  using  only  Channel  2  brightness  temperatures. 

5.1  WATER  VAPOR  BURDEN  AND  THE  SSM/T-2  FREQUENCIES 

The  approach  followed  is  from  the  study  of  Rosenkranz  et  al.  (1982)  who  describe  the 
temperature  weighting  functions  for  both  the  oxygen  band  (e.g.  SSM/T-1  frequencies)  and 
for  the  water  vapor  absorption  line  at  183.31  GHz.  They  define  the  integrated  column 
water  vapor  above  a  given  level  as  the  water  vapor  burden  (U)  which  is  proportional  ^ 

to  r,  the  optical  depth.  Rosenkranz  et  al,  also  noted  that  when  weighting  functions 
are  expressed  as  functions  of  water  vapor  burden  rather  than  as  functions  of  pressure  or 
temperature,  these  wdghting  functions  are  nearly  symmetric  in  shape  and  their  peaks  t 

are  invariant  with  respect  to  changes  in  the  water  vapor  profile.  These  frmetions  are  also 
displaced  downward  as  the  frequency  A^es  away  from  the  absorption  line  at  183  GHz. 

Another  way  of  stating  this  principle  is  that  the  water  vapor  bixrden  may  be  expressed  j 

as: 

U{h)=  r  p{h')dh\  (5.1) 

where  p  (h)  is  the  water  vapor  density  as  a  function  of  height,  h.  This  means  that  the  I 

satelhte  brightness  temperature  is  a  weighted  average  of  the  temperature  profile.  The 
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weighting  is  defined  by  the  weighting  function,  which,  when  expressed  in  terms  of  U, 
is  c^proximately  invariant  with  respect  to  changes  in  p(h).  Therefore,  for  a  temperature 
profile  with  (say)  a  linear  lapse  rate  near  the  peak  of  the  weighting  function,  the  brightness 
temperatiu%  and  thermometric  temperature  of  the  atmosphere  will  be  equal  at  the  peak 
of  the  weighting  function  (Wilheit  and  al  Khalaf,  1993).  This  leads  to  the  hypothesis  that 
the  water  vzq>or  burden  at  the  level  where  the  atmospheric  temperattu^  (Tt)  equals  the 
brightness  temperature  (Tb)  will  be  constant  for  a  given  frequency  independent  of  the 
water  vapor  profile. 

To  test  this  hypothesis,  the  3  atmosph^c  profiles  described  in  Chapter  3  were  used  for 
clear  sky  temperature/pressure  profiles  and  various  water  vapor  profiles  (using  the  water 
vapor  multiplicative  factor)  were  employed  to  yield  47  cases.  The  water  vapor  burden  was 
calculated  using  (5.1)  where  the  height  h  is  the  altitude  where  Tt  =  Tt.  Table  5.1  shows 
the  results  for  these  water  vapor  burdens  for  the  brightness  temperatures  at  183.31  ±  1 
GHz  for  these  47  different  atmospheric  cases. 

Nine  cases  result  m  water  vapor  burdens  of  less  than  0.35  kg/m^  for  the  47  cases  in  , 

Table  5.1.  Certainly,  one  can  argue  whether  case  7  is  a  physically  realistic  atmosphere. 

The  other  cases  with  low  water  vapor  burdens  are  for  small  water  vapor  multiplicative 
factors.  Nonetheless,  the  hypothesis  appears  realistic  for  the  183.31  ±  1  GHz  channel.  For  ^ 

the  cases  considered  in  Table  5.1  a  mean  water  vapor  bmrden  of  0.3735  kg/m^  and  a  root 
mean  square  error  (RMSE)  value  of  0.0256  kg/m^. 

Water  vapor  burdens  were  also  calculated  for  the  same  47  cases  for  183.31  ±  3  GHz. 

I 

These  results  are  listed  in  Table  5.2.  Agmn  the  values  of  water  vapor  burden  are  reasonably 
constant  with  a  mean  and  RMS  values  of  1.2739  kg/m*  and  0.3709  kg/m*,  respectively. 

The  cases  with  very  low  water  vapor  burdens  are  characterized  by  optical  depths  of  less 

than  4,  meaning  the  surface  enussivities  are  influencing  the  brightness  temperatmres.  Com-  * 

puting  these  statistics  for  only  the  40  highest  cases  yields  a  mean  and  RMSE  of  1.4097 

and  0.1486  kg/m*,  respectively. 

In  Table  5.2,  case  8  has  a  sinface  relative  humidity  of  about  14%  with  the  surface  k 

temperature  near  20°C.  Cases  16,  17,  and  18  have  relative  hiunidities  of  25-47%  with  a 
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Ikble  5.1:  Water  Vapor  Overburdens  (U)  at  183  ±  1  Ghz  for  Tb  =  Tf 


HjO  MULT. 
FACTOR 


0.50 


0.76 


TEMPERATURE 

PROFILE 

WATER  VAPOR 
PROFILE 

U 

(kg/m^) 

CWV 

(kg/m») 

Topical 

TVopical 

.4083 

20.61 

TYopical 

IVopical 

.3954 

30.91 

Tropical 

Tropical 

.3910 

41.21 

IVopicai 

Tropical 

.3806 

51.51 

'IVopieal 

Tropical 

.3633 

61.82 

Tropical 

Ttopica! 

IKSlaH 

72.12 

I^pical 

Tropical 

.3389 

82.42 

ML  summer 

ML  summer 

7.32 

ML  summer 

ML  siunmer 

mssM 

14.84 

ML  summer 

ML  summer 

.3848 

21.95 

ML  summer 

ML  summer 

.3917 

2927 

ML  summer 

ML  summer 

.3947 

36.59 

ML  summer 

ML  summer 

.3886 

43.91 

ML  summer 

ML  summer 

.4057 

51.23 

ML  summer 

ML  summer 

.3774 

58.54 

ML  winter 

ML  winter 

mssm 

3.42 

ML  winter 

ML  winter 

.3584 

4.27 

ML  winter 

ML  vrinter 

.3767 

6.41 

ML  winter 

ML  winter 

.3671 

8.54 

ML  winter 

ML  winter 

.3518 

10.68 

ML  winter 

ML  winter 

.3366 

12.81 

ML  vrinter 

ML  winter 

wEsm 

14.95 

ML  winter 

ML  winter 

17.08 

ML  winter 

ML  summer 

.3695 

11.71 

ML  winter 

ML  summer 

.3697 

14.64 

ML  winter 

ML  summer 

.3722 

21.95 

ML  winter 

ML  summer 

.3762 

29.27 

ML  winter 

ML  summer 

M'j.l  liM 

36.59 

ML  winter 

ML  summer 

.3886 

43.91 

ML  winter 

ml  summer 

.3955 

51.23 

ML  winter 

ML  summer 

.4012 

58.54 

ML  summer 

ML  winter 

.3208 

3.42 

ML  summer 

ML  winter 

M'i.l'J  ■ 

4.27 

ML  summer 

ML  winter 

.4090 

6.41 

ML  summer 

ML  winter 

.3888 

8.54 

ML  summer 

ML  winter 

.3645 

10.88 

ML  summer 

ML  winter 

.3396 

12.81 

ML  summer 

ML  winter 

.3219 

11.95 

ML  summer 

ML  winter 

.3086 

17.08 

Tropical 

ML  summer 

.3980 

11.71 

IVopical 

ML  summer 

.4006 

14.64 

Tropical 

ML  summer 

.4007 

21.95 

Tropical 

ML  summer 

.4013 

29.27 

Tropical 

ML  summer 

.3995 

36.59 

Tropical 

ML  summer 

.3886 

43.91 

Tropical 

ml  summer 

.3784 

51.23 

Ttopical 


ML  summer 
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<rurface  temperature  of  -1°C  and  cases  32-34  have  relative  humidities  of  5-11%.  Certsinly, 
cases  32-34  would  be  extremely  rme,  a  surface  temperature  of  about  20°C  and  a  relative 
hiimidity  near  10%  (althou^  these  surface  conditions  are  occasionally  observed  in  Col¬ 
orado,  but  the  surface  pressiire  certainly  isn’t  1000  mb).  The  question  of  whether  these 
atmospheres  are  realistic  or  not  is  debatable  but  if  they  happen,  it  will  be  very  rarely. 

There  is  also  mo-e  scatter  of  the  data  for  the  ±  3  GHz  channel  compared  to  the  ± 
1  GHz  channel.  The  rerults  in  Tables  5.1  and  5.2  correspond  well  with  results  of  similar 
studies  conducted  by  Wilheit  and  al  Khalaf  (1993). 

Figures  5.1  and  5.2  are  plots  of  the  water  vapor  burden  for  the  47  cases  for  183.31 
±  1  and  3  GHz,  respectively  (these  are  the  data  given  in  Tables  5.1  and  5.2).  Figure  5.1 
shows  how  the  vapor  burdens  corresponding  to  the  level  at  Tj,  =  Tt  only  slightly  \'aries 
from  atmosphere  to  atmosphere  whereas  Figure  5.2  shows  that  the  water  vaoor  burden  is 
much  more  spread  although  the  outliers  are  of  a  debatable,  physically,  as  stated  above. 

Figures  5,3~5.6  show  weighting  functions  plotted  in  terms  of  the  bmden  U  for  the 
183.31  ±  1,  3,  7  GHz  and  150.00  GHz  channels,  respectively.  These  figures  grsqihicaliy 
portray  the  essence  of  the  hypothesis  stated  above.  While  the  magnitudes  of  the  weighting 
function  change  for  each  atmospheric  profile,  the  peaks  of  these  timetiens  and  their  shape 
are  remarkably  similar  for  these  different  atmospheric  profiles.  This  relative  independence 
of  the  shape  of  W’(U)  on  both  the  temperature  and  moisture  profiles  provides  a  useful  way 
of  interpreting  the  SSM/T-2  brightness  temperature  distributioiis  such  as  those  shown  in 
Figures  2.2-2.11. 

5.2  CHANNEL  2  WATER  VAPOR  BURDEN  MAX* 

Based  on  the  discussion  and  results  of  the  previous  section,  this  secuou  presents  a 
retrieval  of  the  water  vapor  burden  for  Channel  2  (183.31  ±  1  GHz)  on  a  constant  pres¬ 
sure  surface.  The  data  used  in  this  retrieval  consisted  of  the  T-2  brightness  temperatures 
described  above  and  European  Centre  for  Medium  Range  Weather  Forecasts  (ECMWF) 
monthly  mean  temperature  and  pressure  analyses  for  March  1992  (used  to  simulate  the 
temperature  and  pressure  data  which  will  normally  be  available  from  the  SSM/T-1  Tem¬ 
perature  Sounder).  ECMWF  monthly  mean  data  were  used  because  temperature  data 
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Table  5,2:  Water  Vi^or  Overburden  (U)  at  183.31  ±  3  GHz  for  =  Tt 


CASE 

HaO  MULT. 
FACTOR 

U 

(kg/m^) 

CWV 

(kg/m’) 

1 

0.50 

Tropi'jfil 

mm\  mm 

1.49.36 

20.61 

2 

0.75 

Tropical 

Tropical 

1.5346 

30.91 

3 

1.0 

Tropical 

i.5220 

41.21 

4 

1.25 

Tropical 

Tropical 

1.5389 

51.51 

5 

1.50 

Tropical 

Tropical 

1.5669 

61.82 

6 

1.75 

Tropica! 

Tropical 

1.5923 

72.12 

7 

2.00 

Tropica! 

Tropical 

1.6123 

82.42 

8 

0.25 

ML  simuner 

ml  summer 

0.9340 

7.32 

9 

0.50 

ML  surmner 

ML  siunmer 

1.1354 

14.64 

10 

0.75 

ML  suzuoer 

ML  cummer 

1,4097 

21.95 

11 

1.0 

ML  summer 

ML  summer 

1.4589 

29.27 

12 

1.25 

ML  summer 

ML  summer 

1.4913 

36.59 

13 

1.50 

ML  summer 

ML  summer 

1.5153 

43.91 

14 

1.75 

ML  summer 

ML  summer 

1.5346 

51.23 

15 

2.00 

ML  summer 

ML  summer 

0.2077 

58.54 

16 

0.40 

ML  vrinter 

ML  winter 

0.3S47 

3.42 

17 

0.50 

ML  winter 

ML  winter 

0.8060 

4  27 

18 

0.75 

ML  winter 

ML  winter 

1.0605 

6.41 

19 

1.00 

ML  winter 

ML  winter 

1.1879 

8.54 

20 

1.25 

ML  winter 

ML  winter 

1.2347 

10.68 

21 

1.50 

ML  winter 

ML  winter 

1.2901 

12.81 

22 

MW 

ML  winter 

ML  winter 

1.3126 

14.95 

23 

2.00 

ML  winter 

ML  winter 

1.3126 

17.08 

24 

0.40 

ML  winter 

ML  summer 

1.133! 

11.71 

25 

0.50 

ML  winter 

ML  swnmer 

1.2005 

14.34 

26 

0.75 

ML  winter 

ML  sunomer 

1.3006 

21.95 

27 

1.00 

ML  winter 

ML  summer 

1.3700 

29.27 

28 

1.25 

ML  winter 

ML  summer 

1.4i65 

36.59 

29 

1.50 

ML  winter 

ML  summer 

1.4555 

43.91 

20 

1.75 

ML  winter 

ML  summer 

1.4848 

51.23 

31 

2.00 

ML  winter 

ML  siimmt  - 

1.5121 

58.54 

32 

0.40 

ML  summer 

ML  winter 

0.1346 

3.42 

33 

0.50 

Ml-  summer 

ML  winter 

0.2674 

4.27 

34 

0.75 

ML  summer 

ML  winter 

0.7516 

alB 

35 

1.00 

ML  summer 

ML  winter 

1.1167 

SB 

36 

1.25 

ML  summer 

ML  winter 

1.3064 

10.68 

37 

1.50 

ML  summer 

ML  winter 

1.3914 

12.81 

38 

1.75 

ML  summer 

ML  winter 

1.4259 

39 

2.00 

ML  summer 

ML  winter 

1.4423 

17.08 

40 

0.40 

Tropical 

ML  summer 

1.2619 

11.71 

41 

0.50 

Tropical 

ML  summer 

1.3509 

14.64 

42 

0,76 

Tropical 

ML  summer 

1.4658 

21.95 

43 

1.00 

Tropical 

ML  summer 

1.4916 

29.27 

44 

1.25 

Tropical 

ML  summer 

1.5028 

36.59 

45 

1.50 

Tropical 

ML  summer 

1.5376 

43.91 

46 

1.75 

liopicai 

ML  summer 

1.5730 

51.23 

47 

2.00 

Tropical 

ML  summer 

1.6051 

lliEi£3il 

Water  Vapor  Overburden  (kg/rn 


0  5  10  15  20  25  30  35  40  45 


CASE 


Figure  5.1:  Water  Vapor  Overburdens  at  183.31  ±  1  GEb  for  47  Atmospheric  Cases. 
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183  ±1  GHz 


Figure  5.3:  Water  vapor  burden  for  183  ±  1  GHz.  Atmospheres  are  abbreviated  as  in 
text. 
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Figure  5.4:  Same  as  in  Figure  5.4  except  for  183  ±  3  GHz. 
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Figure  5.6:  The  same  as  Figure  5.3  except  for  150  GHz. 
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coincident  with  the  time  of  the  SSM/T-2  observations  was  unavulable  at  the  time  of  this 
research.  The  results  presented  here,  especially  in  the  tropics,  where  synoptic  temperattire 
variations  are  small,  will  not  greatly  alter  the  results  presented. 

The  retrieval  starts  from  the  assumption  that  the  brightness  temperatmre  and  the 
thermometric  temperature  of  the  atmosphere  are  equal  when  the  vapor  burden  is  0.3735 
kg/m^  (the  mean  value  from  Table  5.1).  Using  the  temperature  and  pressure  data  of 
ECMWF  from  March  1992,  the  pressure  p*  corresponding  to  the  level  where  Tt  =  Ts  is 
derived  by  interpolation.  Maps  of  p*  represent  the  contours  of  the  constant  U  smface. 
High  values  of  p*  correspond  to  dry  conditions  whereas  low  values  of  p*  correspond  to  a 
moist  atmosphere. 

A  plot  of  the  interpolated  p*  derived  in  this  way  is  presented  in  Figure  5.7.  The 
pressmre  varies  from  high  values  in  the  mid-to-high  latitudes  (varying  from  609-500  mb) 
to  re^ons  associated  with  deep  convection  in  the  equatorial  remans  that  have  p*  <  300 
mb.  These  areas  are  found  in  North-central  South  America,  the  West  Pacific  and  the 
South  Pacific  Convergence  Zone  (SPCZ),  Indonesia  extending  west  into  the  Indian  Ocean, 
and  portions  of  Central  Africa.  These  regions  may  be  considered  to  be  the  ascent  region 
of  the  Hadley  circulation  which  pumps  water  vapor  from  low  in  the  atmosphere  to  near 
the  tropopairse.  The  subtropical  high  regions  are  not  as  distinct  on  this  map  as  compared 
to  the  brightness  temperature  maps  although  some  regions  do  show  distinct  maxima  in  p* 
such  as  the  Bay  of  Bengal,  and  in  regions  across  the  northern  and  southern  Pacific.  This 
also  is  not  smprising  since  the  upper  level  water  vapor  is  carried  from  the  equatorial  region 
at  upper  levels  towards  the  poles  so  this  water  vapor  at  upper  levels  would  be  expected 
to  only  drop  off  slowly.  Pressures  increase  when  moving  away  from  the  equator.  This 
pressure  difference  between  the  tropics  and  the  arctic  regions  is  very  large  in  the  Southern 
Hemisphere  whereas  this  difference  is  not  nearly  as  great  in  the  Northern  Hemisphere. 

After  calculating  this  pressme  at  each  grid  point,  the  mean  pressine,  p*,  is  computed. 
A  simple  method  for  interpolating  the  water  vapor  information  contsdned  in  Figure  5.7 
to  a  constant  pressure  surface  follows.  This  surface  is  selected  {is  the  mean  pressme  p* 
between  40°  North  and  South  latitude  which  for  this  case  was  393  mb. 
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of  pressure,  p*,  where  Tf 
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The  water  vapor  burden  for  Channel  2,  U2  has  a  value  of  0.37  kg  m  ^  and  may  be 
e:q>res8ed  as  the  integral: 

U2  =  -  f  rdp,  (5.2) 

9  Jp' 

where  g  is  the  gravitational  constant  and  r  is  the  miying  ratio.  Assiuning  the  Tnivitig  ratio 
can  be  expressed  by:  r  =  malcing  the  substitution  p  =  ^  where  po  is  the 

surface  pressure  and  r.  the  saturation  mixing  ratio,  (5.2)  can  be  rewritten  as: 


(5.3) 


Finally,  making  the  substitution  p*  =  fr  and  integrating  results  in  an  expression  for  U2: 

po 


The  difference  between  U  (p*)  and  U  (p*)  at  each  grid  point  is  now  determined  and 
used  to  adjust  the  column  overburden.  This  difference,  A  U  can  be  written  as  (from  (5.2)): 

p^dp,  (5.5) 

9  Jp> 

which,  on  integration,  becomes: 


On  rearrangement  using  (5.4)  for  U2,  (5.6)  can  be  written: 

=  (T^[(P  )"+'  -  (^)^+^].  (5.7) 

Finally,  recalling  the  definitions  for  p  and  p  and  rearranging  results  in  an  equation  for  A  U 
which  is  added  to  the  value  of  U2  at  each  grid  point: 

(5.8) 

Since  p*,  the  pressure  onto  which  U  is  to  be  mapped,  is  ^ven  (393  mb),  p*  is  derived 
from  the  T-2  brightness  temperature,  smd  U2  is  known  (0.37  kg/m^),  then  A  U  can  be 
readily  determined  and  the  burden  U(p*)  =  U  +  A  U  then  follows.  The  water  vapor 
overburden  U  at  p*  were  derived  in  this  way  and  mapped  (Figure  5.8).  The  water  vapor 
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biirden  on  the  constant  presstue  surface  is  very  high  (~  0.7-0.8  kg/m*)  in  portions  of  the 
tropics  and  extremely  low  (below  0.2  kg/m*)  in  most  areas  north  of  Antarctica.  These 
results  correlate  well  with  the  previous  figure.  The  water  vapor  burden  at  'n^400  mb  would 
be  expected  to  be  larger  at  the  tropics  than  towards  the  poles.  Surprising  is  the  magnitude 
of  the  differences  between  these  remans. 

This  process  is  instructive  when  one  considers  that  part  of  the  operational  output 
from  the  T-2  is  layer  water  VE^or,  as  stated  in  Chapter  2.  Performing  calculations  which 
yield  the  water  vapor  burden  on  several  levels  means  the  layer  water  vapor  mass  can  also 
be  carried  out  in  this  way.  Certainly  the  simplified  approach  shown  above  would  not  be 
\ised  operationally  but  it  does  provide  a  framework  to  understand  how,  in  a  general  sense, 
the  layer  vary  in  the  atmosphere  and  may  provide  a  valuable  way  of  initializing  more 
sophisticated  retrieval  approaches. 


i  r 


s 


P  .. 


¥ 


1  ' 


\ 


> 


IS 
•Ns 

C35 

5  CSJ  CO 

^  w  Q) 

-  C7) 

KQ  t.  »-H 

C  0 
O  c  TJ 

a  CD 
U2  S 

V  cj  jQ  I 

C^si  (D  ^ 

I  > 

ffi  H  O  X 

^  ^  L  I 
$5  0  o 
K  CO  g-s: 


er  vapor  bixrden  on  a  constant  pressure  (393  mb)  surface. 


Chapter  6 

CONCLUSIONS  AND  SUMMARY 

6.1  SUMMARY 

The  opportunity  to  be  among  the  first  researchers  to  work  with  data  from  the  SSM/T- 
2  Microwave  Water  Vapor  Sounder  presented  a  unique  opporhmity  to  examine  some  of  the 
capabilities  of  this  instrument.  Since  the  SSM/T-2  is  a  new  instrument,  the  characteristics 
of  this  instnunent  and  the  other  instrmnents  flown  on  the  DMSP  spacecraft  F-11  were 
reviewed.  The  DMSP  microwave  suite,  (SSM/I,  SSM/T-1,  and  SSM/T-2  instruments), 
and  the  visible  and  infrared  scanners  result  in  a  multi-spectral  platform  from  which  obser¬ 
vations  of  the  earth  may  be  taken,  including  the  first  satellite  platform  with  firequencies 
near  the  183  GH2  absorption  line. 

Remote  sensing  in  the  microwave  re^me  has  occiured  since  the  1960’s.  Techniques 
for  remote  sounding  in  the  microwave  re^on  have  been  proposed  since  the  1970’s  but 
only  since  the  laimch  of  the  SSM/T-2  has  it  become  possible  to  profile  the  distribution 
of  water  vapor.  Since  this  soimding  approach  is  solved  using  weighting  fimctions  the 
properties  of  these  functions  for  the  SSM/T-2  frequendes  are  examined  for  a  number  of 
different  circumstances. 

Understanding  the  response  of  the  measmed  SSM/T-2  brightness  temperatm^s  to  a 
variety  of  factors  (called  sensitivity  effects,  herein)  was  achieved  using  a  number  of  models, 
such  as  a  model  of  the  propagation  of  microwave  radiation  through  the  atmosphere,  models 
of  the  meteorological  state  variables  of  the  atmosphere,  and  simple  models  of  the  charac¬ 
teristics  of  douds  in  the  atmosphere.  Ebcamining  how  brightness  temperatures  respond  to 
different  atmospheric  parameters  provides  conceptual  framework  necessary  to  understand 
the  properties  of  the  SSM/T-2  brightness  temperatures  and  the  weighting  functions  that 
characterize  these  temperatures 
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The  T-2  data  were  analyzed  by  averagiiig  these  data  in  I'’  latitude  x  lon^tude  bins 
over  two  time  periods.  Becatise  of  the  small  swath  width,  the  restilt  of  discarding  data  from 
the  4  scan  scenes  furthest  from  nadir  on  each  side,  some  regions  are  not  sampled  during 
the  two  time  periods  in  March  1993  (10-14  and  16-18).  The  global  maps  of  SSM/T-2 
brightness  temperatures  derived  for  these  two  time  periods  were  examined  and  related  to 
the  gross  characteristics  of  the  general  circulation  of  the  atmosphere.  Additionally,  some 
information  about  the  surface  characteristics  of  earth  and  the  moisture  profiles  above  the 
land  stirfaces  are  discernible  from  these  brightness  temperatures. 

Comparisons  of  brightness  temperatures  from  two  remans  of  the  tropical  Pacific  were  ^ 

presented.  One  of  these  re^ems  was  located  vdthin  the  TOGA  COARE  donudn  and  is 
characterized  by  deep  convection  and  moist  conditions.  The  second  region  was  over  the 
Eastern  Pacific  under  prevailing  subsidence  and  thus  represents  conditions  of  drying.  ^ 

The  concept  of  the  invariance  of  the  wrighting  functions  invariant  with  respect  to 
water  vapor  overburden  is  tested.  It  is  demonstrated  that  the  level  where  Channel  2  and 
Channel  1  brightness  temperatures  match  the  thermodynamic  temperatm-e  are  levels  of  ^ 

^proximately  equal  overburden  irrespective  of  the  atmospheric  profiles  of  temperature 
and  moistme. 

The  pressure  p*,  which  corresponds  to  the  level  where  the  brightness  temperature 

> 

for  Channel  2  equals  the  thermodynamic  temperature,  was  derived  by  interpolation  from 
temperature-pressure  data.  The  water  vapor  burden,  for  Channel  2,  U2,  was  then  mapped 
on  a  constant  pressure  surface  p*  =  393  mb. 

> 

6.2  CONCLUSIONS 

From  the  analyses  of  both  model  simulations  and  observations  performed  in  this  study, 
the  following  conclusions  were  made:  I 

1.  Weighting  functions  at  T-2  frequencies  (and  therefore  brightness  temperatures)  are 


strongly  dependent  on  the  profiles  of  water  vapor.  The  dependence  on  temperature 
is  weaker. 
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2.  Background  (surface)  emissivities  stron^y  affect  brightness  temperatures  at  91  GHz, 
ISO  GHz,  and  at  183  ±  7  GHz.  Scenes  characterized  by  backgrounds  with  high 
emissivities  are  scenes  of  high  brightness  temperatures. 

3.  Liquid  water  clouds  profoundly  affect  bri^tness  temperatures: 

(a)  Increased  doud  liquid  water  content  (LWC)  decreases  brightness  temperatures 
over  high  emissivity  backgrounds. 

(b)  Clouds  generally  increase  brightness  temperatiires  over  low  emissivity  (i.e. 
ocean)  stu^aces  and  decrease  brightness  temperatures  over  high  (i.e.  land)  emis* 
sivity  surfaces. 

4.  An  observational  study  of  the  T-2  brightness  temperature  maps  results  in  the  abil¬ 
ity  to  see  some  characteristics  of  the  known  general  circulation  of  the  atmosphere. 
The  subtropical  remans  with  thdr  assodated  subsidence  and  hence,  high  bright¬ 
ness  temperatures,  are  striking  features  on  the  T-2  brightness  temperature  maps 
for  all  chaimels  on  the  183  GHz  absorption  line.  The  South  Pacific  Convergence 
Zone  (SPCZ)  and  the  Intertropical  Convergence  Zone  (ITCZ)  appear  on  the  maps 
of  Channels  1  and  2  which  are  indicative  of  high  mid-to-upper  tropospheric  water 
vapor. 

5.  Surface  characteristics  can  also  be  observed  in  the  T-2  brightness  temperatiu'e  maps 
for  the  window  channels  (92  and  150  GHz).  Arid  land  areas  possess  large  emissivities 
and  high  brightness  temperatmes.  All  the  land  masses  in  the  Southern  Hemisphere 
^pear  warm  in  contrast  to  the  relatively  cold  adjacent  ocean  re^ons. 

6.  Water  vapor  burdens  are  relatively  constant  for  the  two  SSM/T-2  frequendes  nearest 
the  absorption  line.  The  water  vapor  burden  is  0.37  ±  0.02  kg/m*  for  Channel  2 
and  1.4  ±  0.15  kg/m*  for  Channel  1.  The  presstire  of  the  surface  of  constant  vapor 
overburden,  introduced  here  as  p*,  varies  from  around  600  mb  in  high  latitudes  to 
pressures  less  than  300  mb  in  regions  assodated  with  moist  convection. 
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7.  Brightness  temperatxires  and  temperature-pressure  profiles  obtained  from  ECMWF 
analyses,  combined  with  the  above  result,  were  then  mapped  onto  a  constant  pres8\ire 
surface.  This  provides  a  simple  method  for  retrieving  vapor  overburden.  The  method 
was  illustrated  for  a  constant  393  mb  pressure  surface. 

6.3  DISCUSSION  OF  FUTURE  WORJC 

Obviously,  since  the  SSM/T-2  represents  a  new  microwave  radiometer  on  a  satellite 
there  are  a  number  of  sensitivity  issues  which  still  need  to  be  examined.  The  first  issue 
that  comes  to  mind  is  that  both  mixed  phase  and  ice  clouds  are  not  discussed  very  much 
in  the  literature,  except  that  the  influence  of  ice  clouds  is  ne^gible  at  the  firequencies  of 
the  SSM/T-2.  Mixed-phase  clouds  also  represent  a  significant  hurdle  in  the  interpretation 
of  brightness  temperatures.  Certjunly,  the  ability  to  model  this  problem  and  quemtify  the 
effects  of  ice  and  mixed  phase  clouds  wotild  increase  the  ability  to  retrieve  water  vapor 
information. 

A  validation  test  of  the  retrieval  scheme  introduced  hers  should  also  be  undertaken. 
This  could  be  accomplished  using  data  from  special  intensive  observation  periods  (e.g. 
TOGA  COARE)  where  water  vapor  information  derived  from  radiosonde  observations 
could  be  compared  with  coinddent  T-2  data. 

Improved  models  of  the  surface  emissivity  at  the  183  GHz  region  are  also  needed. 
These  models  should  address  the  emissivities  ocean  surfaces  and  sea  ice  and  snow  surfaces. 
The  possibility  of  using  the  T-2  for  mow  cover  mapping  should  also  be  investigated. 

The  183  GHz  region  is  under  study  by  NASA  in  a  high  priority  program  to  place  a  Mi¬ 
crowave  Imager  and  Sotmder  on  a  future  geostationary  satellite  -  the  GEO-PLATFORM. 
A  microwave  sounder  of  this  type  would  certmnly  improved  the  sampling  capabilities  for 
mesoscale  forecasting. 

I  feel  the  immediate  future  (less  than  five  years)  will  see  an  enormous  amount  of 
research  along  the  lines  of  a  syner:^stic  approach/imified  retrieval.  The  scanners  and 
radiometers  on  DMSP  spacecraft  F-11  provide  a  data  set  which  is  as  comprehensive  as 
any.  Certainly  the  SSM/I  could  be  tised  for  the  ocean  surface  wind  speeds  and  this 
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information  could  be  used  to  help  in  the  detciiuination  of  the  background  enoissivity  over 
ocean  surfaces  for  the  183  GHz  channel?.  Ad^Iitionally,  the  SSM/I  and  SSM/T-2  can  be 
used  together  for  column  water  v«^or  and  cloud  liquid  water  content.  Here  it  may  be 
possible  to  explore  column  water  v!q>or  over  both  hmd  and  ocean  surfaces.  Present  SSM/I 
analysis  of  this  colmnn  quantity  is  restricted  to  ocean  areas  only  The  visible  and  infrared 
data  from  the  OLS  can  be  used  for  cloud  type  determumtion  (and  cloud  cover)  and  the 
responses  of  the  SSM/T-2  channels  can  be  evaluated  to  improve  the  imderstanding  of  the 
eftects  of  clouds  in  the  183  GHz  re^on. 

The  SSM/T-2  will  serve  to  only  expand  the  research  opportuidces  to  those  with  a 
keen  interest  in  water  vapor  and  it  will  undoubtedly  improve  th*;  ability  of  NWP  models  to 
forecast  the  state  of  the  atmosphere.  While  a  great  deal  of  work  remains  to  be  ione  on  the 
use  of  the  183  GHz  absorption  line  for  water  vapor  profiling,  it  has  been  an  enlightening 
opportunity  to  use  the  T-2  data. 
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